Question

MATLAB程序 这个导入之后 运行老是出错 求大佬帮解答 万分感激

x0=[32*pi,4,20,2.92,3,1.2,4];
goal=[2,0.5,0.5,200];
w=[0.25,4,4,0.0001778];
lb=[20*pi,2,15,2.5,1.5,1,3];
ub=[50*pi,6,25,3.5,4.5,2,5];
options=optimset('Largescale','off');
[x,fn,attainfactor,exitflag]=fgoalattain(@eg9_7_mubiao,x0,goal,w,[],[],[],[],lb,ub,@eg9_7_yueshu,options)

function f=eg9_7_mubiao(y)
g=9.8;
km=0.2;
p=7.1*10^3;
q=7850;
h=0.3;
D=y(4);
fy=inline('sqrt((2*pi^2*x^2+cos(2*pi*x)-1)/(2*pi*x-sin(2*pi*x))+1)-D','x','a');
id=fzero(fy,0.8,[],D);
b=2*pi*id;
f(1)=y(1)*y(2)/1000/(2*pi*y(4))*(b-sin(b)-sqrt(y(4)^2-1)*(1-cos(b))+b*b/sqrt(2)*tan(y(3)*pi/180));
f(1)=1/f(1);
f(2)=g*sin(y(3)*pi/180)/y(1)*(sqrt(y(4)^2-1)*(1-cos(b))-b+sin(b));
f(2)=abs(f(2));
f(3)=0.14*(km*p*h*y(5)*y(6)+q*y(5)*y(6))*y(1)^2*y(2)^2/1000/1000/1000;
f(4)=54.95*y(1)^2*y(2)*y(5)*y(6)/y(7)^2/1000;

function[g,ceq]=eg9_7_yueshu(x)
gc=9.8;k=10;
g(1)=x(1)^2*x(2)/1000/gc-k;
mocajiao=20*pi/180;
x3hudu=x(3)*pi/180;
g(2)=sqrt(abs(gc*sin(x3hudu-mocajiao)/(x(2)/1000)))-x(1);
ceq=[];

Answer 1

该代码在matlab2010a可以运行。 运行结果

追问

我运行不了唉 有分成一个主函数和两个字函数吗？