| Labfans是一个针对大学生、工程师和科研工作者的技术社区。 | 论坛首页 | 联系我们(Contact Us) | 
![]()  | 
	
| 		
			
			 | 
		#1 | 
| 
			
			 初级会员 
			
			
			
			注册日期: 2008-06-16 
				
				年龄: 45 
				
					帖子: 11
				 
				
				
				声望力: 18 ![]()  | 
	
	
	
		
		
			
			 
			
			这是我设计中的一个关键例子:(望各位高手帮帮忙) 
		
		
		
		
		
		
		
	设信号为:x(t)=sin(2π20t),在信号中加入方差为0.1的白噪声序列,取采样频率fs=200Hz,采样点数N=200,用Burg算法分析信号。(应该出现谱分裂现象,可是我弄的程序不出现,调了很多次,也没出来,请各信号处理高手指点,小女子将万分感激) 我编的程序如下: clear all fs=200; nfft=200; order=20; t=0:1/fs:1; f1=20; x=sin(2*pi*f1*t)+sqrt(0.1)*randn(size(t)); [p,f]=pburg(x,order,nfft,fs); p1=10*log10(p); plot(f,p1) grid on; xlabel('Frequency(Hz)'); ylabel('相对功率谱密度(dB/Hz)'); title('Burg PSD Estimate');  | 
| 
		 | 
	
	
	
		
		
		
		
			 
		
		
		
		
		
		
		
			
		
		
		
	 | 
| 		
			
			 | 
		#2 | 
| 
			
			 初级会员 
			
			
			
			注册日期: 2008-09-04 
				
				年龄: 41 
				
					帖子: 13
				 
				
				
				声望力: 18 ![]()  | 
	
	
	
		
		
			
			 
			
			看来你真是初学者, 
		
		
		
		
		
		
		
	clear all fs=200; nfft=200; order=20; t=0:1/fs:1; f1=20; x=sin(2*pi*f1*t)+sqrt(0.1)*randn(size(t)); [p,f]=pburg(x,order,nfft,fs); p1=10*log10(p); figure, plot(f,p1) grid on; xlabel('Frequency(Hz)'); ylabel('相对功率谱密度(dB/Hz)'); title('Burg PSD Estimate');  | 
| 
		 | 
	
	
	
		
		
		
		
			 
		
		
		
		
		
		
		
			
		
		
		
	 | 
| 		
			
			 | 
		#3 | 
| 
			
			 初级会员 
			
			
			
			注册日期: 2008-09-04 
				
				年龄: 41 
				
					帖子: 13
				 
				
				
				声望力: 18 ![]()  | 
	
	
	
		
		
			
			 
			
			你必须加个figure语句在plot的上面!!!
		 
		
		
		
		
		
		
		
	 | 
| 
		 | 
	
	
	
		
		
		
		
			 
		
		
		
		
		
		
		
			
		
		
		
	 | 
| 		
			
			 | 
		#4 | 
| 
			
			 初级会员 
			
			
			
			注册日期: 2008-06-16 
				
				年龄: 45 
				
					帖子: 11
				 
				
				
				声望力: 18 ![]()  | 
	
	
	
		
		
			
			 
			
			加与不加没有区别,主要是要出现谱分裂。
		 
		
		
		
		
		
		
		
	 | 
| 
		 | 
	
	
	
		
		
		
		
			 
		
		
		
		
		
		
		
			
		
		
		
	 |