MATLA编程应用实例.docVIP

1.求解线性方程组AX=B 程序： A=[1 1.5 2 9 7;0 3.6 0.5 -4 4; 7 10 -3 22 33;3 7 8.5 21 6;3 8 0 90 -20] B=[3 -4 20 5 16] x=inv(A)*B 结果： ans = 3.5653 -0.9255 -0.2695 0.1435 0.0101 2.计算函数值 程序： x=input(请输入x) if x1 x elseif x10 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 if d0 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() if a+bcabs(a-b)c 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 for n=1:5 x=cos(k*t+((n-1)*(pi/4))); y=cos(k*t); subplot(2,3,n) plot(x,y) if n==1 title(相位差=0) elseif n==2 title(相位差=pi/4) elseif n==3 title(相位差=3pi/2) elseif n==4 title(相位差=3*pi/4) elseif n==5 title(相位差=pi) end xlabel(x) ylabel(y) subplot(2,3,6) plot(x,y) xlabel(x) ylabel(y) hold on end 9.绘制下列极坐标图 （1） （2） （3） 程序： t=0:pi/60:2*pi; polar(t,5*cos(t)+4) ,text(5,5,5cos(\theta)+4) hold on 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) axis equal 14.动画：单摆 程序：for t1=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])