数学模型实践课.PPTVIP

数学模型实验课（一） 求解非线性方程 问题：已知L=30, C=0.5, M=51, 求W, 使得 精确到 0.001 MW - LC - W =0 是关于未知量W的非线性方程 难于解析求解 不必要解析求解 求问题的近似的数值解 一. M文件的编辑与建立 M文件具有普通的文本格式，可以用任何编辑程序建立和编辑。 可以使用Matlab提供的M文件窗口来进行。 在Matlab命令窗口的File菜单中选择New命令， 出现子菜单， 再选择m-file的命令，将得到M文件的编辑窗口 可以编辑M文件：可以是一段程序或一个函数。 编辑完成后，在此窗口的File菜单中选择save as的命令，得到对话框，输入文件名 ***.m 编辑M文件 pap.m function y=pap(W) M=51 L=30 C=0.5 Y=(M.*W-W.*sqrt(C.^2-W.^2)-L*C Save as pap.m 二. 求解非线性方程 图像求解 W=0:0.001:0.5;plot(W,pap(W),’r’,W,0,’—b’) 迭代求解 n=0.25; while pap(n)0, n=n+0.001; end; n 零点求解 Y=fzero(‘pap’，0.2） * *