|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
<P><FONT face=宋体>clear all;<BR>close all;<BR>clc;<BR>%</FONT><FONT face=宋体>产生数据</FONT><BR><FONT face=宋体>N=128*8;<BR>fs=20;<BR>f=0:fs/N:fs/N*(N-1);<BR>n=0:N-1;<BR>f1=0.6381;<BR>fai1=2*pi*rand(1);<BR>%fai1=0.142;<BR>f2=0.8345;<BR>fai2=2*pi*rand(1); <BR>%fai2=1.253;<BR>f3=1.4726;<BR>fai3=2*pi*rand(1);<BR>%fai3=1.413;<BR>f4=0.4909;<BR>fai4=2*pi*rand(1);<BR>%fai4=2.365;<BR>f5=1.7671;<BR>fai5=2*pi*rand(1); <BR>%fai5=3.127;<BR>f6=2.2580;<BR>fai6=fai4+fai5;<BR>%fai6=5.492;<p></p></FONT></P>
<P><FONT face=宋体>%e=rand(1,N);e=e-mean(e);<BR>t=[0:1/fs:(N-1)/fs];<BR>s1=cos(2*pi*f1*t+fai1);<BR>s2=cos(2*pi*f2*t+fai2);<BR>s3=cos(2*pi*f3*t+fai3);<BR>s4=cos(2*pi*f4*t+fai4);<BR>s5=cos(2*pi*f5*t+fai5);<BR>s6=cos(2*pi*f6*t+fai6);<BR>%x=s1+s2+0.5*s3+s1.*s2+e;<BR>x=s1+s2+s3+s4+s5+s6;<BR>x=x-mean(x);<p></p></FONT></P>
<P><FONT face=宋体>%</FONT><FONT face=宋体>双谱估计平画图</FONT><BR><FONT face=宋体>figure(2);<BR>%x=s1+s2+0.5*s3+s1.*s2+e;<BR>[bspec,waxis]=bispecd(x);<BR>title('</FONT><FONT face=宋体>双谱估计平画图</FONT><FONT face=宋体>')<BR>figure(3)<BR>bspec_abs=abs(bspec);<BR>mesh(bspec_abs);<BR>title('双谱估计三维图')<p></p></FONT></P>
<P><FONT face=宋体>%</FONT><FONT face=宋体>累积量切片</FONT><BR><FONT face=宋体>t1=64;<BR>%figure(5)<BR>for k=-5:0.1:5<BR>y_cum = cumest(x,3,t1,128,50,'biased',k,k)<BR>end<BR>%plot(y_cum)<p></p></FONT></P>
<P><FONT face=宋体>%</FONT><FONT face=宋体>切片的fft变换</FONT><BR><FONT face=宋体>nfft=256;<BR>yy=abs(psd(y_cum,nfft));<BR>%yy=abs(fft(y_cum,nfft));<BR>figure(6)<BR>plot((0:nfft/2-1)/nfft*fs,yy(1:nfft/2));<BR>xlabel('</FONT><FONT face=宋体>频率</FONT><FONT face=宋体> f/Hz');<BR>ylabel('fft')<p></p></FONT></P>
<P>上面是我写的高阶谱分辨耦合和非耦合情况的程序,这个程序分辨不出来这两种情况,但是不知道错在哪里,希望高手们指点!</P> |
|