菜鸟提问：方波幅频变换
 

求出方波的频谱[quote]

ft=sym('(t+4)/2*Heaviside(t+4)-t*Heaviside(t)+(t-4)/2*Heaviside(t-4)');
Fw=simplify(fourier(ft));
Ff=subs(Fw, '2*pi*f', 'w');
Ff_conj=conj(Ff);
GF=sqrt(Ff*Ff_conj);
ezplot(GF,[-0.5 0.5]);grid;

[/quote]

网上流传的一段代码



都看得懂，就是第一行看不懂



明明不是方波，乘了t明明

[quote]

ft=sym('Heaviside(t+4)-Heaviside(t)+Heaviside(t-4)');


Fw=simplify(fourier(ft));

回复: 菜鸟提问：方波幅频变换
 

求出方波的频谱[quote]

ft=sym('(t+4)/2*Heaviside(t+4)-t*Heaviside(t)+(t-4)/2*Heaviside(t-4)');
Fw=simplify(fourier(ft));
Ff=subs(Fw, '2*pi*f', 'w');
Ff_conj=conj(Ff);
GF=sqrt(Ff*Ff_conj);
ezplot(GF,[-0.5 0.5]);grid;

[/quote]
网上流传的一段代码

都看得懂，就是第一行看不懂

明明不是方波，乘了t明明
[quote]
ft=sym('Heaviside(t+4)-Heaviside(t)+Heaviside(t-4)');

Fw=simplify(fourier(ft));

[/quote]
但是这样写却报错

本来想用fft，可是我写不出一个离散点形式的方波y，望高手指教！

回复: 菜鸟提问：方波幅频变换
 

抱歉一开始不会编辑帖子

求出方波的频谱[quote]

ft=sym('(t+4)/2*Heaviside(t+4)-t*Heaviside(t)+(t-4)/2*Heaviside(t-4)');
Fw=simplify(fourier(ft));
Ff=subs(Fw, '2*pi*f', 'w');
Ff_conj=conj(Ff);
GF=sqrt(Ff*Ff_conj);
ezplot(GF,[-0.5 0.5]);grid;

[/quote]
网上流传的一段代码

都看得懂，就是第一行看不懂

明明不是方波，乘了t明明
[quote]
ft=sym('Heaviside(t+4)-Heaviside(t)+Heaviside(t-4)');

Fw=simplify(fourier(ft));

[/quote]
但是这样写却报错

本来想用fft，可是我写不出一个离散点形式的方波y，望高手指教！