 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 793 | 回復: 2 | |||
[求助]
基于Wilson理論的風力機葉片設計 已有2人參與
|
|
各位大神，下邊是基于Wilson理論編寫的風力機葉片設計程序，分別是主函數(shù)，目標函數(shù)和約束函數(shù)，但是調試有問題，運行出錯，懇請各位大神幫忙看下原因，不勝感激！ 報錯顯示如下 “錯誤使用 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; |
回復此樓» 猜你喜歡
279分求調劑 一志愿211
已經有12人回復
-
297求調劑
已經有9人回復
-
085700資源與環(huán)境308求調劑
已經有13人回復
-
材料 336 求調劑
已經有4人回復
-
華東師范大學-071000生物學-293分-求調劑
已經有3人回復
-
280求調劑
已經有8人回復
-
324分 085600材料化工求調劑
已經有4人回復
-
一志愿華中科技大學，080502，354分求調劑
已經有5人回復
-
一志愿武理材料305分求調劑
已經有7人回復
-
271材料工程求調劑
已經有8人回復
chendequan
鐵蟲 (小有名氣)
- 應助: 152 (高中生)
- 金幣: 918.6
- 紅花: 17
- 帖子: 262
- 在線: 208.6小時
- 蟲號: 994474
- 注冊: 2010-04-11
- 專業(yè): 結構工程
2樓2016-12-13 16:35:43
pdl9527
專家顧問 (小有名氣)
-
專家經驗: +8 - 應助: 100 (初中生)
- 金幣: 2109.3
- 紅花: 23
- 帖子: 282
- 在線: 112小時
- 蟲號: 1227333
- 注冊: 2011-03-09
- 性別: GG
- 專業(yè): 分離過程
- 管轄: 計算模擬
3樓2016-12-13 20:55:47
相關版塊跳轉
我要訂閱樓主 dreamjie618 的主題更新
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿武理材料305分求調劑 +6 | 想上岸的鯉魚 2026-03-18 | 7/350 |
|
|---|---|---|---|---|
|
|
[考研] 一志愿重慶大學085700資源與環(huán)境專碩，總分308求調劑 +3 | 墨墨漠 2026-03-18 | 3/150 |
|
|
|
[考研] 一志愿西南交大，求調劑 +5 | 材化逐夢人 2026-03-18 | 5/250 |
|
|
|
[考研] 288求調劑 +16 | 于海海海海 2026-03-19 | 16/800 |
|
|
[考研] 085600材料與化工 +8 | 安全上岸！ 2026-03-16 | 8/400 |
|
|
|
[考研]
 材料與化工 322求調劑
+4
|
然11 2026-03-19 | 4/200 |
|
|
|
[考研] 求調劑 +5 | Mqqqqqq 2026-03-19 | 5/250 |
|
|
|
[考研] 0817 化學工程 299分求調劑 有科研經歷 有二區(qū)文章 +22 | rare12345 2026-03-18 | 22/1100 |
|
|
|
[考研] 298-一志愿中國農業(yè)大學-求調劑 +9 | 手機用戶 2026-03-17 | 9/450 |
|
|
|
[考博] 招收博士1-2人 +3 | QGZDSYS 2026-03-18 | 3/150 |
|
|
|
[考研] 307求調劑 +9 | 冷笙123 2026-03-17 | 9/450 |
|
|
|
[考研] 【同濟軟件】軟件（085405）考研求調劑 +3 | 2026eternal 2026-03-18 | 3/150 |
|
|
|
[考研] 化學工程321分求調劑 +15 | 大米飯！ 2026-03-15 | 18/900 |
|
|
|
[考研] 0703化學調劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|
|
|
[考研] 生物學071000 329分求調劑 +3 | 我愛生物生物愛?/a> 2026-03-17 | 3/150 |
|
|
|
[考研] 301求調劑 +4 | A_JiXing 2026-03-16 | 4/200 |
|
|
|
[考研] 考研化學學碩調劑，一志愿985 +4 | 張vvvv 2026-03-15 | 6/300 |
|
|
|
[考研] 考研調劑 +3 | 淇ya_~ 2026-03-17 | 5/250 |
|
|
|
[考研] 277材料科學與工程080500求調劑 +3 | 自由煎餅果子 2026-03-16 | 3/150 |
|
|
|
[考研] 085601材料工程315分求調劑 +3 | yang_0104 2026-03-15 | 3/150 |
|