作者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