 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2819 | 回復(fù): 13 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||
743832049新蟲(chóng) (初入文壇)
|
[求助]
我用matlab編寫了一個(gè)風(fēng)力機(jī)葉片的優(yōu)化設(shè)計(jì)程序，有點(diǎn)問(wèn)題 已有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; %當(dāng)前只進(jìn)行0.1R處翼型截面的氣動(dòng)力計(jì)算試一試 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); 目標(biāo)函數(shù) 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; 條件函數(shù) 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中運(yùn)行的時(shí)候總是出問(wèn)題。我是菜鳥(niǎo)，該程序是邊學(xué)邊做的，可能有出錯(cuò)的地方，望大俠們幫幫吧。。。 |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 學(xué)。 |
» 猜你喜歡
334分 一志愿武理-080500 材料求調(diào)劑
已經(jīng)有10人回復(fù)
-
324分 085600材料與化工
已經(jīng)有7人回復(fù)
-
各位老師好，我的一志愿為北京科技大學(xué)085601材料專碩
已經(jīng)有11人回復(fù)
-
309求調(diào)劑
已經(jīng)有6人回復(fù)
-
調(diào)劑推薦
已經(jīng)有7人回復(fù)
-
一志愿西交大080500材料學(xué)碩349
已經(jīng)有7人回復(fù)
-
一志愿北京科技大學(xué)085601材料工程英一數(shù)二初試總分335求調(diào)劑
已經(jīng)有4人回復(fù)
-
352分-085602-一志愿985
已經(jīng)有6人回復(fù)
-
273求調(diào)劑
已經(jīng)有5人回復(fù)
-
349求調(diào)劑
已經(jīng)有4人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 跪求永磁同步風(fēng)力發(fā)電機(jī)組在ＭＡＴＬＡＢ搭建的仿真模型
已經(jīng)有5人回復(fù)
- [已完結(jié)]利用Isight與matlab結(jié)合優(yōu)化風(fēng)力機(jī)葉片氣動(dòng)外形的例子
已經(jīng)有12人回復(fù)
- 優(yōu)化 matlab編程
已經(jīng)有4人回復(fù)
- 付費(fèi)求助用matlab進(jìn)行多目標(biāo)優(yōu)化設(shè)計(jì)
已經(jīng)有4人回復(fù)
- 求助：遺傳算法的matlab編程求解資源優(yōu)化問(wèn)題
已經(jīng)有7人回復(fù)
- matlab優(yōu)化設(shè)計(jì)方向的課件
已經(jīng)有106人回復(fù)
- 求助使用MATLAB計(jì)算綠色葉片的覆蓋度的問(wèn)題
已經(jīng)有24人回復(fù)
743832049
新蟲(chóng) (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 112.8
- 帖子: 33
- 在線: 20.3小時(shí)
- 蟲(chóng)號(hào): 2370476
- 注冊(cè): 2013-03-23
- 專業(yè): 工程熱物理相關(guān)交叉領(lǐng)域
5樓2013-12-11 15:32:04
2樓2013-12-07 21:00:05
|
本帖內(nèi)容被屏蔽 |
3樓2013-12-08 10:55:54
743832049
新蟲(chóng) (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 112.8
- 帖子: 33
- 在線: 20.3小時(shí)
- 蟲(chóng)號(hào): 2370476
- 注冊(cè): 2013-03-23
- 專業(yè): 工程熱物理相關(guān)交叉領(lǐng)域
4樓2013-12-11 15:30:58
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 311（085601）求調(diào)劑 +9 | liziyeyeye 2026-03-28 | 10/500 |
|
|---|---|---|---|---|
|
|
[考研] 安徽大學(xué)專碩生物與醫(yī)藥專業(yè)(086000)324分，英語(yǔ)已過(guò)四六級(jí)，六級(jí)521，求調(diào)劑 +10 | 美味可樂(lè)雞翅 2026-03-26 | 11/550 |
|
|
[考研] 085601英二數(shù)二求調(diào)劑 總分325 +4 | 余航航 2026-03-31 | 4/200 |
|
|
|
[考研] 311求調(diào)劑 +9 | 藍(lán)月亮亮 2026-03-30 | 9/450 |
|
|
|
[考研]
 一志愿武漢理工，總分321，英一數(shù)二，求老師收留。
+12
|
nnnnnnn5 2026-03-25 | 12/600 |
|
|
|
[考研] 343求調(diào)劑 +8 | 愛(ài)羈絆 2026-03-28 | 8/400 |
|
|
|
[考研] 367求調(diào)劑 +7 | 芋泥啵�！� 2026-03-28 | 7/350 |
|
|
|
[考研] 一志愿北化085600材料專碩275|有文章專利｜求調(diào)劑 +15 | Micky11223 2026-03-25 | 16/800 |
|
|
|
[考研] 313求調(diào)劑 +6 | 賣個(gè)關(guān)子吧 2026-03-31 | 6/300 |
|
|
|
[考研] 284求調(diào)劑 +14 | junqihahaha 2026-03-26 | 15/750 |
|
|
|
[考研] 295求調(diào)劑 +5 | wei-5 2026-03-26 | 5/250 |
|
|
|
[考研] 275求調(diào)劑 +15 | Micky11223 2026-03-25 | 20/1000 |
|
|
|
[考研] 329求調(diào)劑 +10 | 鈕恩雪 2026-03-25 | 10/500 |
|
|
|
[考研] 349求調(diào)劑 +6 | 李木子啊哈哈 2026-03-25 | 6/300 |
|
|
|
[考研] 復(fù)試調(diào)劑 +3 | raojunqi0129 2026-03-28 | 3/150 |
|
|
|
[考研] 292求調(diào)劑 +14 | 鵝鵝鵝額額額額?/a> 2026-03-25 | 15/750 |
|
|
|
[考研] 086502化學(xué)工程342求調(diào)劑 +6 | 阿姨復(fù)古不過(guò) 2026-03-27 | 6/300 |
|
|
|
[考研]
 一志愿 西北大學(xué) 總分282 英語(yǔ)一62 求調(diào)劑
+7
|
18419759900 2026-03-25 | 8/400 |
|
|
|
[論文投稿] Journal of Mechanical Science and Technology +3 | Russ_ss 2026-03-25 | 5/250 |
|
|
|
[考研] 303求調(diào)劑 +6 | 藍(lán)山月 2026-03-25 | 6/300 |
|