close
matlab疊代程式問題
猜a,aa值，帶入式1~6，然後再由式7,8得到新的a,aa值，疊帶多次最後結果應該是收斂到約0.11010.33270.00620.01535.41*10^-30.0244程式碼：cl=0.989;V1=1;%m/sR=0.3;TSR=6;a=0;aa=0;a(1)=0.3;%初始猜值aa(1)=0;%初始猜值B=3;fori=1:100ifabs((a(i+1)-a(i))/a(i))<0.001&&abs((aa(i+1)-aa(i))/a(i))<0.001%收斂條件disp(theta(i));...顯示更多猜a,aa值，帶入式1~6，然後再由式7,8得到新的a,aa值，疊帶多次最後結果應該是收斂到約0.11010.33270.00620.01535.41*10^-30.0244程式碼：cl=0.989;V1=1;%m/sR=0.3;TSR=6;a=0;aa=0;a(1)=0.3;%初始猜值aa(1)=0;%初始猜值B=3;fori=1:100ifabs((a(i+1)-a(i))/a(i))<0.001&&abs((aa(i+1)-aa(i))/a(i))<0.001%收斂條件disp(theta(i));disp(a(i));disp(aa(i));disp(chord(i));disp(beta(i));disp(solidity(i));elsetheta(i)=(atan((1-a(i))/6/(1+aa(i))));%式1beta(i)=theta(i)-(6*pi/180);%式2Cn(i)=0.989*cos(theta(i));%式3ct(i)=0.989*sin(theta(i));%式4chord(i)=8*pi*0.3*aa(i)*sin(theta(i))*cos(theta(i))/(3*ct(i)*(1+aa(i)));%式5solidity(i)=3*chord(i)/(2*pi*0.3);%式6a(i+1)=1/(1+4*sin(theta(i))*sin(theta(i))/(solidity(i)*0.989*cos(theta(i))));%式7aa(i+1)=1/((4*cos(theta(i))/(solidity(i)*0.989))-1);%式8endend最後輸出的收斂結果就是0.11010.33270.00620.01535.41*10^-30.0244麻煩matlab高手幫我看一下謝謝~~~
我重寫了程式碼, 但最後出來的結果是錯的, 是否可以請您再檢查您的程式看看 clc cl = 0.989; V1 = 1; %m/s R = 0.3; TSR = 6; B = 3; % Initial Guess a = 0.3; aa = 0.1; % Iteration Setting k = 0; kmaxit = 100; Tola = 1.0; Tolaa = 1.0; while k < kmaxit && Tola > 0.001 || Tolaa > 0.001 theta = (atan((1 - a)/6/(1 + aa))); %式1 beta = theta - (6*pi/180); %式2 Cn = 0.989*cos(theta); %式3 ct = 0.989*sin(theta); %式4 chord = 8*pi*0.3*aa*sin(theta)*cos(theta)/(3*ct*(1 + aa)); %式5 solidity = 3*chord/(2*pi*0.3); %式6 a_new = 1/(1 + 4*sin(theta)*sin(theta)/(solidity*0.989*cos(theta))); %式7 aa_new = 1/((4*cos(theta)/(solidity*0.989))-1); %式8 Tola = abs((a_new - a)/a); Tolaa = abs((aa_new - aa)/aa); a = a_new; aa = aa_new; k = k + 1; end theta = (atan((1 - a)/6/(1 + aa))); %式1 chord = 8*pi*0.3*aa*sin(theta)*cos(theta)/(3*ct*(1 + aa)); %式5 beta = theta - (6*pi/180); %式2 solidity = 3*chord/(2*pi*0.3); %式6disp(theta) disp(a) disp(aa) disp(chord) disp(beta) disp(solidity) ----------------------------------結果如下 2.1182e-013 1.0000 0.1000 1.3726e-007 -0.1047 2.1846e-007
以上文章來自奇摩知識家，如有侵犯請留言告知
https://tw.answers.yahoo.com/question/index?qid=20110119000016KK07913