已阅读5页,还剩18页未读, 继续免费阅读
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
MATLAB 语言与数学运算实验课程任务书一、实验教学目标与基本要求上机实验是本课程重要的实践教学环节。实验的目的不仅仅是验证理论知识,更重要的是通过上机加强学生的实验手段与实践技能,掌握应用MATLAB 语言求解问题的方法,培养学生分析问题、解决问题、应用知识的能力和创新精神,全面提高学生的综合素质。上机学时为4学时,主要内容是基于上课内容对课后典型习题应用MATLAB进行求解,基本掌握常见数学问题的求解方法与命令调用,更深入地认识和了解MATLAB语言的强大的计算功能。上机实验最终以书面报告的形式提交,作为期末成绩的考核内容。二、实验内容(4 学时)第一部分MATLAB 语言编程、科学绘图与基本数学问题求解(4 学时)主要内容:掌握MATLAB 语言编程、科学绘图、微积分、线性代数问题等基本数学问题的求解与应用。练习题:1、安装MATLAB 软件,应用demo 命令了解主要功能,熟悉基本功能,会用help 命令。2、用MATLAB 语句输入矩阵A 和B 矩阵前面给出的是4 4 矩阵,如果给出A(5,6) = 5 命令将得出什么结果?答案:3、假设已知矩阵A,试给出相应的MATLAB 命令,将其全部偶数行提取出来,赋给B 矩阵,用A =magic(8) 命令生成A 矩阵,用上述的命令检验一下结果是不是正确。答案:4、用数值方法可以求出,试不采用循环的形式求出和式的数值解。由于数值方法采用double 形式进行计算的,难以保证有效位数字,所以结果不一定精确。试采用符号运算的方法求该和式的精确值。5、答案;答案:答案:1、2、3、答案:答案:ans =1编写一个vander()函数如下:function A=vander(v)n=length(v);v=v(:);A=sym(ones(n);for j=n-1:-1:1,A(:,j)=v.*A(:,j+1);end将此函数存放在 c:MATLAB6p5work下,然后开始运行工作命令即可。需要编写funm()与diagm()函数如下:function F=funm(A,fun,x) V,J=jordan(A); v1=0,diag(J,1); v2=find(v1=0), length(v1)+1; for i=1:length(v2)-1 v_lambda(i)=J(v2(i),v2(i); v_n(i)=v2(i+1)-v2(i); end m=length(v_lambda); F=sym(); for i=1:m J1=J(v2(i):v2(i)+v_n(i)-1,v2(i):v2(i)+v_n(i)-1); fJ=funJ(J1,fun,x); F=diagm(F,fJ); end F=V*F*inv(V);function A=diagm(A1,A2)A=A1; A(size(A1,1)+1:size(A1,1)+size(A2,1),size(A1,2)+1:size(A1,2)+size(A2,2)=A2;function fJ=funJ(J,fun,x)lam=J(1,1); f1=fun;fJ=subs(fun,x,lam)*eye(size(J);H=diag(diag(J,1),1); H1=H;for i=2:length(J) f1=diff(f1,x); a1=subs(f1,x,lam); fJ=fJ+a1*H1; H1=H1*H/i;end第二部分数学问题求解与数据处理(4 学时)主要问题:掌握代数方程与最优化问题、微分方程问题、数据处理问题的MATLAB 求解方法。练习题:5、用数值求解函数求解下述的方程的根,并对得出的结果进行检验。OPT=optimset;OPT.Gradobj=on; 注:求解精度le-10,le-20 不好用可以不要。写OPT语句即可,得到相同结果:OPT=optimset写目标函数: function y=exc6fun6(x) y=exp(x(1)*(4*x(1)2+2*x(2)2+4*x(1)*x(2)+2*x(2)+1);写约束条件:function c,ce=exc6fun6a(x)ce=;c=x(1)+x(2); x(1)*x(2)-x(1)-x(2)+1.5; -10-x(1)*x(2);A=; B=; Aeq=; Beq=; xm=-10; -10; xM=10; 10;x0=(xm+xM)/2;ff=optimset; ff.TolX=1e-10; ff.TolFun=1e-20; x=fmincon(exc6fun6,x0,A,B,Aeq,Beq,xm,xM,exc6fun6a,ff)Warning: Large-scale (trust region) method does not currently solve this type of problem,switching to medium-scale (line search). In C:MATLAB6p5toolboxoptimfmincon.m at line 213Maximum number of function evaluations exceeded;increase OPTIONS.MaxFunEvalsx = 0.41950.4195编写一个M文件如下:(不要加上function语句)i=1; x=x0;while (1)x,a,b=fmincon(exc6fun6,x,A,B,Aeq,Beq,xm,xM,exc6fun6a,ff);if b0, break; endi=i+1;end给M函数起个名字如:kk,回到工作页调用。 kk x,ix = 1.1825 -1.7398i = 533f=-592 381 273 55 48 37 23; A=3534 2356 1767 589 528 451 304;B=119567; intlist=1;1;1;1;1;1;1;ctype=-1; xm=zeros(7,1);xM=inf*ones(7,1); res,b=ipslv_mex(f,A,B,intlist,xM,xm
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 机电调试施工施工方法及工艺要求
- 观光车辆项目节能分析报告
- 母婴产业园项目人力资源管理分析报告
- 轻量化铸件生产项目运营管理
- 传统文化实践活动设计与策划 课件 活动10 悟国画意蕴习书法技艺
- 消防安全培训的标准与目标
- 企业市场营销与工作总结
- 春节习俗英语作文
- 安徽省2024九年级英语全册专项十六词语运用课件新版人教新目标版
- 企业工作总结及团队建设
- 35kV复测报告模板(1)
- 叉车安全管理操作规程
- 运营铁路隧道的病害整治
- 《年度考核证明模板 》
- 10以内加减法练习题(4套200道)
- 数控加工中心进给系统设计 毕业设计
- 广播体操 强身健体ppt课件
- 东阳工业用地项目准入评审办法
- 低压断路器的选用和整定原则及方法
- 大学生早餐调查报告
- GB∕T13912-2002金属覆盖层钢铁制件热浸镀锌层技术要求及试验方法
评论
0/150
提交评论