数学实验第三次讲稿new.pptVIP

具体方法如下: 建立函数m文件: function ff funpara x,a ff sin a*x -x; fsolve 语句的用法 建立相应的命令m文件: B zeros 100,1 ; for a 1:100 x0 pi/ 2*a +0.01; B a fsolve @funpara,x0,[],a ; end plot B * * * * 例：某物体的边缘呈圆形，通过测量边沿上是一个点的坐标，数据为（见程序）。使用fsolve计算物体边缘的曲线方程。 解：假定物体边缘的曲线方程为（x-a）^2+ y-b ^2 r^2 为了得到a，b和r的值。将上述十一组数据代人方程，得到 （xi-a）^2+ yi-b ^2 r^2（i 1,…,11）。次数方程个数是11个，未知量个数是3个。使用fsolve求解，其函数m文件为： function y funcircle x A [6.7630 5.1313 2.4713 -0.3435 -2.3887 -2.9927 -1.9572 0.3778 3.2455 5.7042 6.9465]; B [23.2879 25.6492 26.7268 26.1668 24.1531 21.3470 18.6699 17.0010 16.8883