MATLAB编程应用实例.docx

MATLAB编程应用实例.docx

《MATLAB编程应用实例.docx》由会员分享，可在线阅读，更多相关《MATLAB编程应用实例.docx（14页珍藏版）》请在冰豆网上搜索。

MATLAB编程应用实例

1.求解线性方程组AX=B

程序：

A=[11.5297;03.60.5-44;710-32233;378.5216;38090-20]

B=[3-420516]'

x=inv（A）*B

结果：

ans=

3.5653

-0.9255

-0.2695

0.1435

0.0101

2.计算函数值

程序：

x=input（'请输入x'）

ifx<1

x

elseifx<10

2*x-1

else

3*x-11

End

结果：

第二题

输入x39

x=

39

ans=

106

3.求一元二次方程

程序：

a=input（''）

b=input（''）

c=input（''）

d=b*b-4*a*c

ifd<0

disp（'方程无解'）

else

（-b+d）/（2*a）

（-b-d）/（2*a）

end

结果：

输入a4

a=4

输入b8

b=8

输入c2

c=2

d=32

ans=3

ans=-5

4.输入三角形三条边，求面积

程序：

a=input（''）

b=input（''）

c=input（''）

ifa+b>c&abs（a-b）

p=（a+b+c）/2

s=sqrt（p*（p-a）*（p-b）*（p-c））

end

结果：

输入三角形三边3

a=3

输入三角形三边4

b=4

输入三角形三边5

c=5

p=6

s=6

5.绘图：

程序：

x=-100:

pi/100:

100;

y=100./（1+x.*x）;

plot（x,y）

6.

程序：

x=linspace（-10,10,100）;

y=exp（-（x.*x）/2）/（2*pi）;

plot（x,y）

7.

程序：

t=linspace（-10,10,1000）;

x=t.^2;

y=5.*t.^3;

plot（x,y）

8.李萨如图形

程序：

t=0:

0.01:

2*pi;

k=1

forn=1:

5

x=cos（k*t+（（n-1）*（pi/4）））;

y=cos（k*t）;

subplot（2,3,n）

plot（x,y）

ifn==1

title（'相位差=0'）

elseifn==2

title（'相位差=pi/4'）

elseifn==3

title（'相位差=3pi/2'）

elseifn==4

title（'相位差=3*pi/4'）

elseifn==5

title（'相位差=pi'）

end

xlabel（'x'）

ylabel（'y'）

subplot（2,3,6）

plot（x,y）

xlabel（'x'）

ylabel（'y'）

holdon

end

9.绘制下列极坐标图

（1）

（2）

（3）

程序：

t=0:

pi/60:

2*pi;

polar（t,5*cos（t）+4）,text（5,5,'5cos（\theta）+4'）

holdon

polar（t,12./（sqrt（t）））;text（-5,-5,'12/sqrt（\theta）'）

polar（t,pi/3.*（t.^2））;text（-5,3,'\pi\theta^2/3'）

10.某工厂2005年度各季度生产值（单位：

万元）分别为：

450.6，395.9,410.2,450.9，绘制饼图

程序：

pie（[450.6,395.9,410.2,450.9]）

11.绘图：

程序：

t=0:

pi/30:

8*pi

x=cos（t）

y=sin（t）

z=t

plot3（x,y,z）

12.

程序：

u=0:

pi/30:

2*pi

v=0:

pi/30:

2*pi

x=（1+sin（u））.*cos（v）

y=（1+cos（v））.*sin（u）

z=sin（u）

plot3（x,y,z）

13.半径为10的球面

程序：

[x,y,z]=sphere（10）

surf（10*x,10*y,10*z）

axisequal

14.动画：

单摆

程序：

fort1=0:

pi/60:

4*pi;

a=sin（t1）;

y1=-10*cos（a）,x1=-10*sin（a）;

t=0:

pi/30:

2*pi

x=cos（t）+x1,y=sin（t）+y1

plot（x,y）,fill（x,y,'b'）,axis（[-20,20,-20,20]）;

line（[x1,0],[y1,20]）

m=getframe

end

movie（m）

15.有一电子在电场中运动，画出轨迹图

（1）沿x轴以初速度为0匀速运动，电场为匀强电场，动画

（2）沿x轴以初速度为0匀速运动，电场为周期变化的余弦电场，动画

（3）画出

（1）

（2）的曲线

（4）在复合电磁场中，电场方向为x方向，磁场方向垂直于电场方向，都为匀强场。

程序：

t=0:

0.01:

6*pi;

v0x=1;

v0y=1;

a1=1;

y1=v0y*t+0.5*a1*t.*t;

x1=v0x*t;

a2=cos（t）;

y2=v0y*t+0.5*a2.*t.*t;

x2=v0x*t;

subplot（2,2,3）

plot（x1,y1,'g',x2,y2,'r'）

gridon

title（'叠加'）

legend（'匀强电场','周期变化电场',3）

x3=t;

y3=cos（t）;

z3=sin（t）;

fori=1:

length（t）

subplot（2,2,1）

plot（x1（1:

i）,y1（1:

i））

axis（[0,8,0,30]）

gridon

holdon

subplot（2,2,1）

plot（x1（i）,y1（i）,'r.','markersize',20）

title（'匀强电场'）

gridon

holdoff

subplot（2,2,2）

plot（x2（1:

i）,y2（1:

i））

axis（[0,8,-30,30]）

title（'周期变化电场'）

gridon

holdon

subplot（2,2,2）

plot（x2（i）,y2（i）,'r.','markersize',20）

gridon

holdoff

subplot（2,2,4）

plot3（x3（1:

i）,y3（1:

i）,z3（1:

i））

holdon

gridon

axis（[0,10,-1,1,-1,1]）

title（'匀强电场和磁场'）

subplot（2,2,4）

plot3（x3（i）,y3（i）,z3（i）,'r.','markersize',20）

view（-100,500）

holdoff

M（i）=getframe;

end

movie（M,1）

16.抛射体的射程与抛射角、初速度的关系

设某一抛射体的初速度为

时，抛射角为

，将其运动在X、Y轴上进行正交分解，则有

水平方向速度为

=

竖直方向速度为

=

gt

质点的坐标为（x,y）,所以x（t）=

y（t）=

g

，则质点的轨迹运动方程为y=x

抛射体能达到的最大高度为H=

其到达最大高度所需时间为T=

空中飞行时间为t=2T=

，抛射体的最大射程为X=

所以抛射体的射程跟初速度

和抛射角

有关，他们之间的关系由MATLAB的绘图功能可以更直观的体现出来。

其程序为：

x=linspace（0,pi/2,100）;%产生行向量发射角

g=10;%重力加速度

v1=10;%初速度取10

v2=15;%初速度取15

v3=20;%初速度取20

v4=25;%初速度取25

y1=v1^2*sin（2*x）/g;%初速度为10下的射程

y2=v2^2*sin（2*x）/g;%初速度为15下的射程

y3=v3^2*sin（2*x）/g;%初速度为20下的射程

y4=v4^2*sin（2*x）/g;%初速度为25下的射程

subplot（2,2,1）;%选择2*2个区的一号区

plot（x,y1）;%输出初速度为10下的射程曲线

title（'v0=10'）;%加图形标题

text（pi/4,10,'射程为10'）;%在最大射程处加图形说明

subplot（2,2,2）;%选择2*2个区的二号区

plot（x,y2）;%输出初速度为15下的射程曲线

title（'v0=15'）;%加图形标题

text（pi/4,22.5,'射程为22.5'）;%在最大射程处加图形说明

subplot（2,2,3）;%选择2*2个区的三号区

plot（x,y3）;%输出初速度为20下的射程曲线

title（'v0=20'）;%加图形标题

text（pi/4,40,'射程为40'）;%在最大射程处加图形说明

subplot（2,2,4）;%选择2*2个区的四号区

plot（x,y4）;%输出初速度为25下的射程曲线

title（'v0=25'）;%加图形标题

text（pi/4,62.5,'射程为62.5'）;%在最大射程处加图形说明

程序运行结果如下图所示: