作者issac10383 (幻，直覺影像)
看板Electronics
標題[問題] matlab程式請各位大師幫幫我找錯誤
時間Tue Oct 8 21:51:35 2013
雖然是簡單的程式但找了一天的我，還是不知道程式哪裡了錯誤
是假設值還是架構上出問題，還請各位老手幫幫忙
感激不盡。
%<風力機模型>
%CP值參照論文的12個測試點
P = 1.225;
Awt = 12;
Vw = [ 0.3 0.25 0.22 0.23 0.19 0.17 0.15 0.13 0.12 0.08 0.05 0.025 ];
CP = 3;
for i = 1:12
Pwt(i) = 1/2*P*CP(i)*Vw(i)*Awt;
Pwt = Pwt
%<太陽能模型>
%Esun是假設值
%η=N
Apv = 100;
Npv = 0.15;
NINV = 0.95;
NWIRE = 0.9;
Esun = [ 1 2 3 4 5 6 7 8 9 10 11 12];
Epv(i) = Apv*Esun(i)*Npv*NINV*NWIRE;
Epv = Epv
%<電池模型>
%系統 ED=[ET(i)-ER(i)]
%ET和ER是假設值
NN = 8760;
ET = [5 6 7 8 9];
ER = [1 2 3 4 5];
EL = [ 0 1 2 3 4 5 6 7 8 9 10 11 12];
Ndischar = 1;
Nchar = 1;
ED(i) = [ET(i)-ER(i)];
if (ED(i)<0)
EBat(i) = EBat(i-1)*NINV*NWIRE*Ndischar-EL(i)+ED(i);
end
if (ED(i)>0)
EBat(i) = EBat(i-1)*Nchar+ED(i);
end
if (ED(i)==0)
EBat(i) = EBat(t-1)+ED(i);
end
ED = ED
end
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 220.136.234.200
1F：噓 LaLuLa:第一，這裡應該不是code debug板 10/08 22:48
2F：→ LaLuLa:第二，你連你打算做甚麼都不說，誰知道你在幹嗎？ 10/08 22:48
3F：噓 LaLuLa:第三，你自己都覺得簡單，錯了還不肯自己找，無解 10/08 22:50
4F：推 bxxl:設一些中斷點,看跑到那個點的時候,變數值跟你想的是否一樣. 10/08 22:54
5F：→ bxxl:這程式這麼短,就算你每一行的結果都印出來看也能找到吧 10/08 22:55
6F：→ bxxl:還是說你自己也不曉得中間過程的變數應該是什麼值? 10/08 22:56
7F：推 bxxl:還有個很大的問題,你沒有說明這程式的"錯誤"是怎樣的型態? 10/08 23:00
8F：→ bxxl:是會出現錯誤訊息?還是可以執行完但是某些值結果不如預期? 10/08 23:00
9F：→ s4300026:借標題問一下，debug功能是不是有什麼限制? 10/09 09:27
10F：→ s4300026:之前我寫小程式時可以用，但我把它寫成一堆函式後 10/09 09:27
11F：→ s4300026:就沒辦法再使用debug功能了... 就只能用沒加";"的方法 10/09 09:28
12F：→ s4300026:和跳出error的方式看變數的變化與錯誤的位置 10/09 09:29