| 查看: 883 | 回复: 2 | |||
[求助]
基于Wilson理论的风力机叶片设计 已有2人参与
|
|
各位大神，下边是基于Wilson理论编写的风力机叶片设计程序，分别是主函数，目标函数和约束函数，但是调试有问题，运行出错，恳请各位大神帮忙看下原因，不胜感激！ 报错显示如下 “错误使用 barrier (line 22) Objective function is undefined at initial point. Fmincon cannot continue. 出错 fmincon (line 799) [X,FVAL,EXITFLAG,OUTPUT,LAMBDA,GRAD,HESSIAN] = barrier(funfcn,X,A,B,Aeq,Beq,l,u,confcn,options.HessFcn, ... 出错 Wilson (line 19) [x,fval]=fmincon(@objfun,x0,[],[],[],[],[],[],@confun);” Wilson.m clear; global r; global R; global phi; global lamb; V1=61.1;P=20000;lamb0=7;Cp=0.4;eta1=0.95;eta2=0.92;B=2; Mu=17.9*10^-6;Cl=0.8361;Cd=0.0171; D=sqrt(P/0.49/V1^3/Cp/eta1/eta2); %D=round(D); d=roundn(D,2); R=d/2; N=60*lamb0*V1/pi/d; omiga=V1*lamb0/R; fid=fopen('Wilson.txt','w'); for i=0.072:0.072:0.72 r=i;a=0;b=0; lamb=lamb0*r/R;phi=atan((1-a)*V1/(1+b)/omiga/r); x0=[lamb0^2/2000+0.027*lamb0-0.038;10/lamb0/exp(12*r/R)]; [x,fval]=fmincon(@objfun,x0,[],[],[],[],[],[],@confun); a=x(1); b=x(2); fprintf(fid,'%f\t%f\n',a,b); end fclose(fid); p=load('Wilson.txt'); fid=fopen('wilson1.txt','w'); for i=1:10 r=i; a=p(i,1); b=p(i,2); phi=atan(1-a)*V1/(1+b)/omiga/r; F=2/pi*acos(exp(-B/2*(R-r)/r*sin(phi))); S=8*pi*a*F*(1-a*F)*sin(phi)*sin(phi)/(1-a)^2/cos(phi); C=S*r/B/Cl; Cp=2/lamb0^2*b*(1-a)*F*lamb^4; fprintf(fid,'%f\t%f\t%f\t%f\n',F,phi,C,Cp); end fclose(fid); objfun.m function f=objfun(x) global r; global R; global phi; global lamb; B=2;lamb0=7; c=B/2*(R-r)/r/sin(phi);F=2/pi*acos(exp(-c)); f=-8/lamb0^2*x(2)*(1-x(1))*F*lamb^3; confun.m function[g,ceq]=confun(x) global r; global R; global phi; global lamb; B=2; c=B/2*(R-r)/r/sin(phi);F=2/pi*acos(exp(-c)); g=[]; ceq=x(1)*(1-x(1)*F)-x(2)*(1+x(2))*lamb^2; |
回复此楼» 猜你喜欢
咱们一起用铁证分析2026国家社科基金中标与否
已经有17人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有47人回复
-
各位道友，我要去昆明玩几天，回来见。
已经有6人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有11人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有4人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有8人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有37人回复
-
好奇怪的filecode
已经有7人回复
-
filecode
已经有9人回复
-
请教兼职经验
已经有3人回复
chendequan
铁虫 (小有名气)
- 应助: 152 (高中生)
- 金币: 918.6
- 红花: 17
- 帖子: 262
- 在线: 208.7小时
- 虫号: 994474
- 注册: 2010-04-11
- 专业: 结构工程
2楼2016-12-13 16:35:43
pdl9527
专家顾问 (小有名气)
-
专家经验: +8 - 应助: 100 (初中生)
- 金币: 2110.3
- 红花: 23
- 帖子: 282
- 在线: 112小时
- 虫号: 1227333
- 注册: 2011-03-09
- 性别: GG
- 专业: 分离过程
- 管辖: 计算模拟
3楼2016-12-13 20:55:47
