数学模型实验课(一),求解非线性方程,问题:已知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),