资源描述:
《数学建模matlab练习4》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库。
1、MATLAB练习四1.绘制两个一元函数的曲线图(例如y二sin(x)、y二cos(x)),加上标注说明,包括坐标、抬头、曲线说明等。%lianxi41.mx=0:pi/30:2*pi;y=sin(x);plot(xzy,*-ob1z'markersize1,10);holdonyl=cos(x);plot(x,yl,'-.dm','markersize1,10);legend(1sinx1,'cosx1);title(1sin和cos函数的曲线图')xlabel(1x*);ylabel(1y1);2.寻找一组数据拟合多项式,写出M文件,包括拟
2、合的次数,数值的计算,数学表达式的呈现,图形的呈现等。%lianxi42x=0:pi/20:pi;y=sin(x);n=input(*n=');p=polyfit(x,y,n);xl=0:pi/20:pi*2;yl=sin(xl);y2=polyval(p,xl)plot(xlzyl,'-ob1,xl,y2z':dm*);legend(1原曲线J】拟合曲线「);xlabel(1x1);ylabel(*y1);title(1sin的拟合曲线和原曲线');axis([-0.2,6.5,-1.2,1.2])1.求解一个代数方程,写出M文件,应该包括
3、函数式M文件和命令式M文件。%fc3.mfunctionf=fc3(x)f=13*x-l=01;%lianxi43clc;clear;close;sprintf('原方程和Jt解析解的结果为:,)s=solve(fc3)sprintf(•保留10位有效数字的结果为:Jdigits(10)sprintf('数值解的结果为:Tx=vpa(s)2.求解一个代数方程组,写出M文件,应该包括函数式M文件和命令式M文件。%fcl.mfunctiony=fcl(x)y(1)=x(l)-0•7*sin(x(1))-0.2*cos(x(2));y(2)=x(2
4、)-0.7*cos(x(l))+0.2*sin(x(2));y=[y(1)y(2)];%lianxi44clcclear;close;sprintf('输入初始值:')x0=input('x0=')%[0.50.5];sprintf('原方程的解为:,)x=fsolve(*fcl',xO)1.求解一个常微分方程,写出M文件,应该包括函数式M文件和命令式M文件。functiondy=fc(x)dy=1Dy=-5*x1;closeclearclcs=dsolve(fc,1x1)2.求解一个常微分方程组,写出M文件,应该包括函数式M文件和命令式M文
5、件。%fc.mfunctiondy=fc(t,y)dy=[y(2);3*(l-y(l)x2)*y(2)-y(1)];%lianxi.mclc;clear;close;[T/Y]=ode45('fc,,[020],[2;3]);plot(T,Y(:,1),心*,T,Y(:,2)宀」)title(9DE方程的解,)xlabel('T1)ylabel(1Y1)legend('Y1'z'Y21)3.输入一个矩阵,求出其行列式的值、逆矩阵、特征值、特征向量。%lianxi47clearclca=input「矩阵A—)%矩阵的输入必须带[]%a=[111
6、1;1234;13610;141020];det=det(a)%对应行列式的值inv=inv(a)%矩阵的逆[vzd]=eig(a)%矩阵的特征值和特征向量trace=trace(a)%矢巨卩车白勺迹rank=rank(a)%矢冃卩车的秩poly=poly(a)%特征多项式cond=cond(a)%条件数1.绘制一个二元函数的曲面图(例如y=sin(x+y)),加上标注说明,包扌舌坐标、抬头、曲面说明、网格、着色等。%lianxi48clear;close;clc;x=0:pi/30:pi*2;y=0:pi/30:pi*2;[X,Y]=mes
7、hgrid(x,y);%将向量x,y指定的区域转化为矩阵X,Yz=sin(X+Y);colormap(cool);mesh(X,Y,z)%生成网格图xlabel(1x轴1);ylabel('y轴*);zlabel('z轴');title('sin(x+y)的曲线图T;legend('sin(x+y)曲面');1.给出一个优化问题,写出M文件,进行求解。%lianxi491clear;close;clc;sprintf(,输入目标函数系数:T;f=input(1f=');%f=[-5,4,2];sprintf(•输入系数矩阵:»);A=inpu
8、t(1A=');%[6-11;124];sprintf(,输入右侧常数:
9、;b=input(1b=1);%[8z10];vlb=[-1,0,0];vub=[3,2]