求出方波的频谱
	引用:
	
	
		
			
				 
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; 
			
		 | 
	
	
 网上流传的一段代码
都看得懂,就是第一行看不懂
明明不是方波,乘了t明明
	引用:
	
	
		
			
				ft=sym('Heaviside(t+4)-Heaviside(t)+Heaviside(t-4)');  
 
Fw=simplify(fourier(ft));  
			
		 | 
	
	
 但是这样写却报错
本来想用fft,可是我写不出一个离散点形式的方波y,望高手指教!