请教matlab求解非线性方程组的问题

shootboy712

我有fsolve求解非线性方程组，方程组中有一个系数是变化的，c是变量，程序如下：
function f = myfun(p,c)
a=40;
angle1=85*pi/180;
angle5=25*pi/180;
f=[(a+p(1))^2+(p(1)-a)^2-2*(p(1)*p(1)-a*a)*cos(angle5)-(2*c*sin(angle1*0.5))^2;...
   (p(1)-a)^2+(2*c*sin(angle1/2))^2-2*(p(1)-a)*2*c*sin(angle1/2)*cos(p(2)+(pi-angle1)/2)-(a+p(1))^2;...
   -p(3)^2+(p(1)-a)^2+c^2-2*(p(1)-a)*c*cos(p(2))]
%%求解过程
p0=[170;1.3;111];
options=optimset('Display','off')
c = 65:10:105;
for i=1:1:length(c);
cc=c(i);
p = fsolve(@ (p) myfun(p,cc),p0,options);
p1(i)=p(1);
p2(i)=p(2);
p3(i)=p(3);
end
可是运行时提示我??? p = fsolve(@ (p) myfun(p,cc),p0,options);          |
Error: "identifier" expected, "(" found.
这是怎么回事呢，谢谢各位了

ch_j1985

原帖由 shootboy712 于 2008-5-23 16:59 发表
我有fsolve求解非线性方程组，方程组中有一个系数是变化的，c是变量，程序如下：
function f = myfun(p,c)
a=40;
angle1=85*pi/180;
angle5=25*pi/180;
f=[(a+p(1))^2+(p(1)-a)^2-2*(p(1)*p(1)-a*a)*cos(angle5 ...

试试这段代码：

1. 
   
2. %函数
   
3. function f = myfun(p)
   
4. global c
   
5. a=40;
   
6. angle1=85*pi/180;
   
7. angle5=25*pi/180;
   
8. f=[(a+p(1))^2+(p(1)-a)^2-2*(p(1)*p(1)-a*a)*cos(angle5)-(2*c*sin(angle1*0.5))^2;...
   
9.    (p(1)-a)^2+(2*c*sin(angle1/2))^2-2*(p(1)-a)*2*c*sin(angle1/2)*cos(p(2)+(pi-angle1)/2)-(a+p(1))^2;...
   
10.    -p(3)^2+(p(1)-a)^2+c^2-2*(p(1)-a)*c*cos(p(2))];
    
11. end

复制代码

1. 
   
2. %主程序
   
3. clc
   
4. clear
   
5. global c
   
6. p0=[170;1.3;111];
   
7. options=optimset('Display','off');
   
8. cc = 65:10:105;
   
9. for i=1:1:length(cc);
   
10. c=cc(i);
    
11. p = fsolve(@myfun,p0,options);
    
12. p1(i)=p(1);
    
13. p2(i)=p(2);
    
14. p3(i)=p(3);
    
15. end

复制代码