| 查看: 2801 | 回复: 13 | |||
743832049新虫 (初入文坛)
|
[求助]
我用matlab编写了一个风力机叶片的优化设计程序，有点问题 已有1人参与
|
||
|
主程序 clear; global r; global R; global phi; global lamb; V1=7;P=400;lamb0=6;B=3; D=2.5; R=D/2; N=60*lamb0*V1/pi/D; omiga=V1*lamb0/R; r=0.1*R; %当前只进行0.1R处翼型截面的气动力计算试一试 a0=lamb0^2/2000+0.027*lamb0-0.038; b0=10/lamb0/exp(12*r/R); lamb=lamb0*r/R; x0=[a0;b0]; for i=1:10 lb=[0;0]; ub=[0.5;1]; options=optimset('largescale','off','Algorithm','active-set'); [x,fval]=fmincon('objfun',x0,[],[],[],[],lb,ub,'confun',options); a=x.(1); b=x.(2); xi=x; if((abs(a-a0)<=0.00001)&&(abs(b-b0)<=0.00001)),break,end x0=xi; end C0=0.21897; for i=0:10 Re=[100000 500000 1000000 2000000 3000000 4000000 5000000 9000000]; Cl=[1.2424 1.0853 1.001 0.9616 0.9663 0.9693 0.9016 1.0023]; alpha=[8 6 5 4.5 4.5 4.5 4 5]; Re(i)=V1*(1-a)*C0/1.48e-5/sin(phi); Cl(i)=interpi(Re,Cl,Re(i),'spline'); alpha(i)=interpi(Re,alpha,Re(i),'spline'); S=8*pi*sin(phi)*sin(phi)*(1-a*F)*a*F/cos(phi)/(1-a)^2; C(i+1)=S*r/B/Cl(i); C=C(i+1); if(abs(C-C0)<=0.00001),break,end C0=C(i+1); end theta=phi-alpha(i); printf(fid,'%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f\n',a,b); 目标函数 function fval=objfun(x0) global lamb0; global B; global R; global r; global phi; lamb=lamb0*r/R; phi=atan((1-x0.(1))/(1+x0.(2))/lamb); f=B*(R-r)/2/r/sin(phi); F=2*acos(exp(-f))/pi; fval=-8*x0.(2)*(1-x0.(1))*F*lamb^3/lamb0^2; 条件函数 function [c,ceq]=confun(x0) global lamb0; global B; global R; global r; global phi; lamb=lamb0*r/R; phi=atan((1-x0.(1))/(1+x0.(2))/lamb); f=B*(R-r)/2/r/sin(phi); F=2/pi*acos(exp(-f)); ceq=x0.(1)*(1-x0.(1)*F)-x0.(2)*(1+x0.(2))*lamb^2; c=[]; 程序采用的是Wilson算法，但是在MATLAB中运行的时候总是出问题。我是菜鸟，该程序是边学边做的，可能有出错的地方，望大侠们帮帮吧。。。 |
回复此楼» 收录本帖的淘帖专辑推荐
 学。 |
» 猜你喜欢
308求调剂
已经有4人回复
-
NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗
已经有14人回复
-
材料与化工一志愿南昌大学327求调剂推荐
已经有6人回复
-
化学调剂0703
已经有7人回复
-
327求调剂
已经有11人回复
-
调剂
已经有8人回复
-
梁成伟老师课题组欢迎你的加入
已经有7人回复
-
伙伴们，祝我生日快乐吧
已经有24人回复
-
中科院材料273求调剂
已经有3人回复
-
材料工程专硕274一志愿211求调剂
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 跪求永磁同步风力发电机组在ＭＡＴＬＡＢ搭建的仿真模型
已经有5人回复
- [已完结]利用Isight与matlab结合优化风力机叶片气动外形的例子
已经有12人回复
- 优化 matlab编程
已经有4人回复
- 付费求助用matlab进行多目标优化设计
已经有4人回复
- 求助：遗传算法的matlab编程求解资源优化问题
已经有7人回复
- matlab优化设计方向的课件
已经有106人回复
- 求助使用MATLAB计算绿色叶片的覆盖度的问题
已经有24人回复
|
本帖内容被屏蔽 |
9楼2013-12-12 20:47:19
2楼2013-12-07 21:00:05
|
本帖内容被屏蔽 |
3楼2013-12-08 10:55:54
743832049
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 112.8
- 帖子: 33
- 在线: 20.3小时
- 虫号: 2370476
- 注册: 2013-03-23
- 专业: 工程热物理相关交叉领域
4楼2013-12-11 15:30:58
743832049
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 112.8
- 帖子: 33
- 在线: 20.3小时
- 虫号: 2370476
- 注册: 2013-03-23
- 专业: 工程热物理相关交叉领域
5楼2013-12-11 15:32:04
|
本帖内容被屏蔽 |
6楼2013-12-11 22:14:27
743832049
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 112.8
- 帖子: 33
- 在线: 20.3小时
- 虫号: 2370476
- 注册: 2013-03-23
- 专业: 工程热物理相关交叉领域
7楼2013-12-12 16:46:36
|
本帖内容被屏蔽 |
8楼2013-12-12 20:34:58
743832049
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 112.8
- 帖子: 33
- 在线: 20.3小时
- 虫号: 2370476
- 注册: 2013-03-23
- 专业: 工程热物理相关交叉领域
10楼2013-12-13 21:52:44
