Matlab程序设计(中高级).pdf

祁彬彬 2005 年年9 月月，全国大学生数学建模比赛成功参赛奖全国大学生数学建模比赛成功参赛奖 2006 年9月，全国大学生数学建模比赛北京市二等奖 2007 年2 月，MCM提名奖 2007 年9月，全国研究生数学建模国家三等奖 2008 年9月，全国研究生数学建模国家一等奖 20099 年年99 月，全国研究生数学学建模国家二等奖 个人能力个人能力 朋友圈子 思考问题的方式 感谢老师 界面设置界面设置 编辑和调试 解方程 画图 数据读写 拟合 积分 优化优化 Lingo小例子（圆的问题） DDesktktopDDesktktop LLayoutt‐DDeffaultlt Preference设置字体，属性等 Set Path添加路径信息 MM 文件文件：函数文件和脚本文件函数文件和脚本文件 多行注释：1.Ctrl+R,2.if 0 end,3.%{   %} 自动缩进：Ctrl+I 双击显示所有变量 做标记Ctrl+F2 更多右键信息息 设置断点，F5,F10,F11,Go until Cursor 分块执行分块执行%%%% 线性方程组线性方程组 非线性方程 隐函数方程 非线性方程组 带参数的情况 微分方程 更多可参考《高等应用数学问题的Matlab 求解》 AAx=BB e−3t sin(4t +2) +4e−0.5t cos(2t) 0.5 fsolve 基本用法 x = fsolve(fun,x0) fzero 基本用法 xx  ==  fzerofzero(fun(fun,x0)x0) 更多用法请看帮助文件，选项和输出 根据图像选取根据图像选取 ⎧⎧ 22 22 1717x −1616 || x || y +1717y 225225 ⎨ 2 2 2 x +y ⎩⎩xx cos(cos(xx ++yy )) ++yy ee 00 画图获得初值的范围画图获得初值的范围 用fsolve 三维以上三维以上，就得根据经验设置初值就得根据经验设置初值，可以通过循环多可以通过循环多 次计算，找到最优的解 cllear;cllc;cllose allll f{1}=17*x^2‐16*abs(x)*y+17*y^2‐225; f{2}=x^2*cos((x+y^2))+y^2*exp((x+y)); hold on; h=arrayfun(@(x)ezplot(f{x},[‐5 5 ‐5 5]),[1,2]); set((h((1),),color,,r)) （（‐33 ，，‐11）） function y=myfun_1(x) y=[[17*x((1)).^2‐16*abbs((x((1)))).*x((2))+17*x((2)).^2‐225; x(1).^2.*cos(x(1)+x(2).^2)+x(2).^2.*exp(sum(x))];  fsolve(@myfun_1,[‐3,‐1]) Equation solved. fsolve completed because the vector of function values is near zero as measured by the default value of the function tolerance, and the problem appears regular as measured by the gradient. stoppiing  criiteriia  ddetailils ans = ‐2.93894670984629         ‐1.16856745664698 ⎧⎧ 22 22 22 x +ax +66bb +33y 00