声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2080|回复: 1

[其他] 请教高阶谱的二次相位耦合

[复制链接]
发表于 2006-6-22 09:45 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

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>
回复
分享到:

使用道具 举报

发表于 2013-8-1 10:56 | 显示全部楼层
看的眼花
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-5 23:30 , Processed in 0.065567 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表