《matlab实习》word版

《《matlab实习》word版》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、1 矩阵运算 矩阵建立，转置、加、减、乘、求逆等运算。计算A+B，AB、

2、A

3、、、AB.>>A+Bans=043569342>>det(A)ans=12>>A.^-1ans=1.00000.50000.33330.25000.20000.16671.00000.33331.0000>>ABans=1.1667-1.83331.50000.66670.66670-1.16670.8333-0.5000>>A*Bans=77913192146102.编写M文件的编写与运行输出闰年i=1;for(j=1000:3000)if(rem(j,100)==0)if(rem(j,400)==0)a(i)

4、=ji=i+1;endelseif(rem(j,4)==0)a(i)=ji=i+1;endendk=i-1end部分数据：Columns463through4692908291229162920292429282932Columns470through4762936294029442948295229562960Columns477through4832964296829722976298029842988Columns484through48529922996 3 数据文件的建立与读取 例：建立班级一门课成绩数据库，编程读取并排序、另存。学号123456789成绩85968587566596

5、8575S.sc=[859685875665968575],S.num=[123456789]S=sc:[859685875665968575]S=sc:[859685875665968575]num:[123456789]fori=1:9forj=i:9ifS.sc(i)

6、(2*t)),polar(t,r,'k')第二次命令t=(-1:0.01:1)*pi/4；r=sqrt(25*cos(2*t)),polar(t,r,'k')参数方程p=-2*pi:0.01:2*pi;r=sqrt(cos(2*p));x=r.*cos(p);y=r.*sin(p);plot(x,y)螺线a=2;t=0:.1:20;r=a*t;polar(t,r)球面X,Y]=meshgrid([-1:.05:1]);%生成网格点，从-1到1间隔0.05Z=sqrt(1-X.^2-Y.^2)+eps;%求出每个点上对应的Z值Z((X.^2+Y.^2)>1)=0;%在满足x^2+y^2>1的网格

7、点上的Z值变成0mesh(X,Y,abs(Z))%画图title('球面x.^2+y.^2+z.^2=1')%加标题axis([-12-1101])%设置X,Y,Z轴坐标范围>>5 多项式运算 例：多项式运算求值、求根。 ++x-3=0a=[111-3]a=111-3>>r=poly(a)r=10-68-3>>q1=poly2str(a,'x')q1=x^3+x^2+x–3f=inline('x^3+x^2+x-3','x')f=Inlinefunction:f(x)=x^3+x^2+x-3>>f(3)ans=366 曲线拟合 例：数据的多项式拟合X0.51.01.52.02.53.0y1.7

8、52.453.814.807.008.60使用了MATLAB中curvefitting，进行了该多形式的拟合。由图可知这些点是五阶多项式linearmodelPoly5:f(x)=p1*x^5+p2*x^4+p3*x^3+p4*x^2+p5*x+p6Coefficients:p1=-1.6p2=13.74p3=-44.07p4=65.67p5=-42.63p6=11.35 7 .线性规划 例：线性规划求解（如教材中）。Min-5x1-4x2-6x3s.tx1-x2+x3<=20;3x1+2x2+4x3<=42;3x1+2x2<=30;x1、x2、x3>=0A=[1-11;324;320]A=1

9、-11324320>>b=[20;42;30]b=204230>>[x,fval,exitflag,output,lambda]=linprog(f,A,b,[],[],lb)Optimizationterminated.x=0.000015.00003.0000fval=-78.0000exitflag=1output=iterations:6algorithm:'large-scale:int