[問題] fft已回收
我想在command window裡面寫一個每週期取32點的fft
利用兩個迴圈
可是不管怎麼跑出來都怪怪的>"<
已經debug兩天了>"<
拜託各位板大幫忙看一下
感激不盡~~~
t是時間、v是電壓、i是電流
目前只做到電壓的部份就卡住了...
clc
load t_v_i_on.txt
t=(t_v_i_on(:,1))'
v=(t_v_i_on(:,2))'
i=(t_v_i_on(:,3))'
m=[0 sin(pi/16) sin(2*pi/16) sin(3*pi/16) sin(4*pi/16) sin(5*pi/16)
sin(6*pi/16) sin(7*pi/16) sin(8*pi/16) sin(9*pi/16) sin(10*pi/16)
sin(11*pi/16) ...
sin(12*pi/16) sin(13*pi/16) sin(14*pi/16) sin(15*pi/16) sin(16*pi/16)
sin(17*pi/16) sin(18*pi/16) sin(19*pi/16) sin(20*pi/16) sin(21*pi/16) ...
sin(22*pi/16) sin(23*pi/16) sin(24*pi/16) sin(25*pi/16) sin(26*pi/16)
sin(27*pi/16) sin(28*pi/16) sin(29*pi/16) sin(30*pi/16) sin(31*pi/16)]
n=[1 cos(pi/16) cos(2*pi/16) cos(3*pi/16) cos(4*pi/16) cos(5*pi/16)
cos(6*pi/16) cos(7*pi/16) cos(8*pi/16) cos(9*pi/16) cos(10*pi/16)
cos(11*pi/16) ...
cos(12*pi/16) cos(13*pi/16) cos(14*pi/16) cos(15*pi/16) cos(16*pi/16)
cos(17*pi/16) cos(18*pi/16) cos(19*pi/16) cos(20*pi/16) cos(21*pi/16) ...
cos(22*pi/16) cos(23*pi/16) cos(24*pi/16) cos(25*pi/16) cos(26*pi/16)
cos(27*pi/16) cos(28*pi/16) cos(29*pi/16) cos(30*pi/16) cos(31*pi/16)]
for k=1:131
a_1=0
b_1=v(1,k)
for j=1:32
v_j(j)=v(1,k+j)
a(k)=(a_1+v_j(j)*(m(j))')/16
b(k)=(b_1+v_j(j)*(n(j))')/16
end
end
t_k=t(1,1:131)
plot(t_k,a,t_k,b)
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 120.126.8.178
推
06/01 13:19, , 1F
06/01 13:19, 1F
討論串 (同標題文章)