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.26950.14350.01012.计算函数值程序：x=input('请输入x')ifx1xelseifx102*x-1else3*x-11End结果：第二题输入x39x=39ans=1063.求一元二次方程0a2cbxx程序：a=input('')b=input('')c=input('')d=b*b-4*a*cifd0disp('方程无解')else(-b+d)/(2*a)(-b-d)/(2*a)end结果：输入a4a=4输入b8b=8输入c2c=2d=32ans=3ans=-54.输入三角形三条边，求面积程序：a=input('')b=input('')c=input('')ifa+bc&abs(a-b)cp=(a+b+c)/2s=sqrt(p*(p-a)*(p-b)*(p-c))end结果：输入三角形三边3a=3输入三角形三边4b=4输入三角形三边5c=5p=6s=65.绘图：21100yx程序：x=-100:pi/100:100;y=100./(1+x.*x);plot(x,y)6.2221yxe程序：x=linspace(-10,10,100);y=exp(-(x.*x)/2)/(2*pi);plot(x,y)7.35xtyt程序：t=linspace(-10,10,1000);x=t.^2;y=5.*t.^3;plot(x,y)8.李萨如图形程序：t=0:0.01:2*pi;k=1forn=1:5x=cos(k*t+((n-1)*(pi/4)));y=cos(k*t);subplot(2,3,n)plot(x,y)ifn==1title('相位差=0')elseifn==2title('相位差=pi/4')elseifn==3title('相位差=3pi/2')elseifn==4title('相位差=3*pi/4')elseifn==5title('相位差=pi')endxlabel('x')ylabel('y')subplot(2,3,6)plot(x,y)xlabel('x')ylabel('y')holdonend9.绘制下列极坐标图（1）4c5os（2）12（3）23程序：t=0:pi/60:2*pi;polar(t,5*cos(t)+4),text(5,5,'5cos(\theta)+4')holdonpolar(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.绘图：tztytsincosx程序：t=0:pi/30:8*pix=cos(t)y=sin(t)z=tplot3(x,y,z)12.uzvuyvuxsinsin)cos1(cos)cos1(程序：u=0:pi/30:2*piv=0:pi/30:2*pix=(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)axisequal14.动画：单摆程序：fort1=0:pi/60:4*pi;a=sin(t1);y1=-10*cos(a),x1=-10*sin(a);t=0:pi/30:2*pix=cos(t)+x1,y=sin(t)+y1plot(x,y),fill(x,y,'b'),axis([-20,20,-20,20]);line([x1,0],[y1,20])m=getframeendmovie(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')gridontitle('叠加')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])gridonholdonsubplot(2,2,1)plot(x1(i),y1(i),'r.','markersize',20)title('匀强电场')gridonholdoffsubplot(2,2,2)plot(x2(1:i),y2(1:i))axis([0,8,-30,30])title('周期变化电场')gridonholdonsubplot(2,2,2)plot(x2(i),y2(i),'r.','markersize',20)gridonholdoffsubplot(2,2,4)plot3(x3(1:i),y3(1:i),z3(1:i))holdongridonaxis([0,10,-1,1,-1,1])title('匀强电场和磁场')subplot(2,2,4)plot3(x3(i),y3(i),z3(i),'r.','markersize',20)view(-100,500)holdoffM(i)=getframe;endmovie(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;%初速度取10v2=15;%初速度取15v3=20;%初速度取20v4=25;%初速度取25y1=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');%在最大射程处加图形说明程序运行结果如下图所示: