matlab程序设计实验34109new

《matlab程序设计实验34109new》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、MATLAB程序设计实验解方程和最优化问题求解姓名：龙刚学号：1404110427班级：电信1104班一．实验目的1.掌握线性方程组的数值求解方法。2.掌握常微分方程的数值求解方法。3.掌握非线性方程以及最优化问题的求解方法。二．实验内容1.对下列线性病态方程组：a)求方程组的解。b)将方程组右边的向量元素b3（0.52）修改为0.53，再次求方程组的解。比较b3的变化和解的相对变化。c)求系数矩阵的条件数，并分析条件数对线性方程组求解的影响。%1.aA=[1/21/31/4;1/31/41/5;1/41/51/6];B=[0.950.6

2、70.52]';disp('1.a方程组的解为：')xa=AB%1.bC=[0.950.670.53]';disp('1.b方程组的解为：')xb=ACchazhi=xa-xb%1.ct1=cond(A)t2=cond(B)t3=cond(C)1.求方程的解。a)funx.m文件functionfx=funx(x)fx=x.^41+x.^3+1;调用z=fzero(@funx,-1)得b)fun1.m文件如下：functionm=fun1(x)m=x-sin(x)./x;a)，fun2.m文件如下：functionf=fun2(p)x

3、=p(1);y=p(2);z=p(3);f(1)=sin(x)+y.^2+log(z)-7;f(2)=3.*x+2.^y-z.^3+1;f(3)=x+y+z-5;2.求常微分方程的数值解。a)，fun3.mfunctionyp=fun3(t,y)yp=-(1.2+sin(10*t)*y);调用代码：结果：a)，fun4.mfunctionyp=fun4(t,y)yp=-y/(1+t^2)+cos(t);调用函数如下：结果：1.求函数在指定区间的最大值。，functionf=fmin(x)f=(1+x.^4)./(1+x.^2);则他的最大

4、值为：2.有400万资金，要求在4年内使用完，若在一年内使用资金x万元，则可得效益万元（效益不能再次使用），当年不用的资金可以存入银行，年利率为10%。制订资金的使用方案，使4年效益之和最大。效益y四年的使用金额分别为=程序代码：x0=[200;200;200];A=[1,0,0;2.1,1,0;2.21,2.1,1];b=[400;440;484];lb=[0;0;0];[x,f]=fmincon(@fun5,x0,A,b,[],[],lb)