24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

743832049

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 112.8
帖子: 33
在线: 20.3小时
虫号: 2370476
注册: 2013-03-23
专业: 工程热物理相关交叉领域

[求助] 我用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中运行的时候总是出问题。我是菜鸟，该程序是边学边做的，可能有出错的地方，望大侠们帮帮吧。。。

回复此楼

» 猜你喜欢

孩子确诊有中度注意力缺陷已经有12人回复
2025冷门绝学什么时候出结果已经有3人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复
AI论文写作工具：是科研加速器还是学术作弊器？已经有3人回复
2026博士申请-功能高分子，水凝胶方向已经有6人回复
论文投稿，期刊推荐已经有4人回复
硕士和导师闹得不愉快已经有13人回复
请问2026国家基金面上项目会启动申2停1吗已经有5人回复
同一篇文章，用不同账号投稿对编辑决定是否送审有没有影响？已经有3人回复

1楼2013-12-04 10:28:26

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

743832049

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 112.8
帖子: 33
在线: 20.3小时
虫号: 2370476
注册: 2013-03-23
专业: 工程热物理相关交叉领域

引用回帖:

9楼: Originally posted by yqsky610 at 2013-12-12 20:47:19
我初步测试发现你贴出来的程序没有input输入指令，fmincon中的子程序objfun和confun要用@，不是 ' '，fmincon要求你的两个子程序返回一个标量值，但是你把x0这个矩阵带进去了，显示错误。
暂时看到这么多，你看看 ...

想把金币给你，你再回复我一下，点击“确定回帖应助”。