 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 792 | 回復(fù): 2 | ||||
dreamjie618新蟲 (初入文壇)
|
[求助]
基于Wilson理論的風(fēng)力機(jī)葉片設(shè)計(jì) 已有2人參與
|
|
各位大神，下邊是基于Wilson理論編寫的風(fēng)力機(jī)葉片設(shè)計(jì)程序，分別是主函數(shù)，目標(biāo)函數(shù)和約束函數(shù)，但是調(diào)試有問題，運(yùn)行出錯(cuò)，懇請(qǐng)各位大神幫忙看下原因，不勝感激！ 報(bào)錯(cuò)顯示如下 “錯(cuò)誤使用 barrier (line 22) Objective function is undefined at initial point. Fmincon cannot continue. 出錯(cuò) fmincon (line 799) [X,FVAL,EXITFLAG,OUTPUT,LAMBDA,GRAD,HESSIAN] = barrier(funfcn,X,A,B,Aeq,Beq,l,u,confcn,options.HessFcn, ... 出錯(cuò) 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; |
回復(fù)此樓» 猜你喜歡
317求調(diào)劑
已經(jīng)有12人回復(fù)
-
一志愿武漢理工材料工程專碩調(diào)劑
已經(jīng)有9人回復(fù)
-
【考研調(diào)劑】化學(xué)專業(yè) 281分，一志愿四川大學(xué)，誠(chéng)心求調(diào)劑
已經(jīng)有7人回復(fù)
-
324求調(diào)劑
已經(jīng)有5人回復(fù)
-
288求調(diào)劑
已經(jīng)有16人回復(fù)
-
317求調(diào)劑
已經(jīng)有9人回復(fù)
-
287求調(diào)劑
已經(jīng)有8人回復(fù)
-
一志愿 西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。
已經(jīng)有4人回復(fù)
-
求調(diào)劑，一志愿:南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分
已經(jīng)有4人回復(fù)
-
085600材料與化工
已經(jīng)有8人回復(fù)
chendequan
鐵蟲 (小有名氣)
- 應(yīng)助: 152 (高中生)
- 金幣: 918.6
- 紅花: 17
- 帖子: 262
- 在線: 208.6小時(shí)
- 蟲號(hào): 994474
- 注冊(cè): 2010-04-11
- 專業(yè): 結(jié)構(gòu)工程
2樓2016-12-13 16:35:43
pdl9527
專家顧問 (小有名氣)
-
專家經(jīng)驗(yàn): +8 - 應(yīng)助: 100 (初中生)
- 金幣: 2109.3
- 紅花: 23
- 帖子: 282
- 在線: 112小時(shí)
- 蟲號(hào): 1227333
- 注冊(cè): 2011-03-09
- 性別: GG
- 專業(yè): 分離過程
- 管轄: 計(jì)算模擬
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
dreamjie618: 金幣+20, ★★★很有幫助 2016-12-14 15:18:30
感謝參與，應(yīng)助指數(shù) +1
dreamjie618: 金幣+20, ★★★很有幫助 2016-12-14 15:18:30
|
請(qǐng)檢查你的目標(biāo)函數(shù)和初始值，你可以將你目前的初值帶入目標(biāo)函數(shù)中，在對(duì)x0賦值的下一行輸入objfun(x0)，運(yùn)行后可以看到得到的復(fù)數(shù)，改變x0的值，發(fā)現(xiàn)還是得到的復(fù)數(shù)，因此fmincon無法繼續(xù)進(jìn)行了。 另外不要用global了，直接用內(nèi)嵌函數(shù)就行。 |
3樓2016-12-13 20:55:47
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 dreamjie618 的主題更新
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 288求調(diào)劑 +16 | 于海海海海 2026-03-19 | 16/800 |
|
|---|---|---|---|---|
|
[考研] 一志愿蘇州大學(xué)材料求調(diào)劑，總分315（英一） +5 | sbdksD 2026-03-19 | 5/250 |
|
|
|
[考研] 350求調(diào)劑 +5 | weudhdk 2026-03-19 | 5/250 |
|
|
|
[考研] 求調(diào)劑一志愿南京航空航天大學(xué)289分 +3 | @taotao 2026-03-19 | 3/150 |
|
|
|
[考研] 316求調(diào)劑 +5 | 梁茜雯 2026-03-19 | 5/250 |
|
|
|
[考研] 265求調(diào)劑 +12 | 梁梁校校 2026-03-19 | 13/650 |
|
|
|
[考研] 本人考085602 化學(xué)工程 專碩 +19 | 不知道叫什么！ 2026-03-15 | 21/1050 |
|
|
|
[考研] 260求調(diào)劑 +3 | 朱芷琳 2026-03-20 | 3/150 |
|
|
|
[考研] 299求調(diào)劑 +3 | shxchem 2026-03-20 | 4/200 |
|
|
|
[考研] 288求調(diào)劑，一志愿華南理工大學(xué)071005 +5 | ioodiiij 2026-03-17 | 5/250 |
|
|
|
[考研] 一志愿福大288有機(jī)化學(xué)，求調(diào)劑 +3 | 小木蟲200408204 2026-03-18 | 3/150 |
|
|
|
[考研] 304求調(diào)劑 +6 | 司空. 2026-03-18 | 6/300 |
|
|
|
[考研] 311求調(diào)劑 +4 | 冬十三 2026-03-18 | 4/200 |
|
|
|
[考研] 299求調(diào)劑 +5 | △小透明* 2026-03-17 | 5/250 |
|
|
|
[考研] 材料，紡織，生物（0856、0710），化學(xué)招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|
|
|
[考研] 考研調(diào)劑 +3 | 淇ya_~ 2026-03-17 | 5/250 |
|
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國(guó)防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|
|
|
[考研] 085601材料工程315分求調(diào)劑 +3 | yang_0104 2026-03-15 | 3/150 |
|
|
|
[考研] 297一志愿上交085600求調(diào)劑 +5 | 指尖八千里 2026-03-14 | 5/250 |
|
|
|
[考研] 一志愿哈工大材料324分求調(diào)劑 +5 | 閆旭東 2026-03-14 | 5/250 |
|