声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1512|回复: 1

[编程技巧] 这个程序出问题了..那位大侠帮忙修改一下

[复制链接]
发表于 2007-5-23 16:27 | 显示全部楼层 |阅读模式

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

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

x
clear;
m=sqrt(-1);
delta=0.101043;
a1=-0.850848;
sample=32;      %number of sample spot
p=10;              %number of sample spot in coef method;
f1=0.05; f2=0.40; f3=0.42;
fstep=0.01;
fstart=-0.5;
fend=0.5;
f=fstart:fstep:fend;
nfft=(fend-fstart)/fstep+1;
%un=urn+juin
urn= normrnd(0,delta/2,1,sample);
uin= normrnd(0,delta/2,1,sample);
un=urn+m*uin;
%计算 zn
for n=1:sample-1
    zn(1)=un(1);
    zn(n+1)=-a1*zn(n)+un(n+1);
end
   
%计算 xn
    for n=1:sample
        xn(n)=2*cos(2*pi*f1*(n-1))+2*cos(2*pi*f2*(n-1))+2*cos(2*pi*f3*(n-1))+sqrt(2)*real(zn(n));
    end
   
%*********************************************************
x=xn';
for k=0:1:sample-1
    s=0;
    for n=1:sample-k,
       s=s+conj(x(1,n))*x(1,n+k); %calculate the value of rxx  
    end
   rxx(1,k+1)=(1/sample)*s;
end
Rx=zeros(sample,sample);
Rx=toeplitz(rxx(1,1:32));
[U,S,V]=svd(Rx);
Pmusicf=zeros(1,1/fstep+1);
ei=zeros(1,sample);
     for i=1:length(f)
            for j=1:sample
                 ei(1,j)=exp(-2*pi*(j-1)*f(i)*m);
            end;
            sum=0;
            for k=7:sample
                 sum=sum+abs(ei*V(:,k))^2;
            end
     Pmusicf(1,i)=10*log10(1/sum);
     end
   figure
plot(f,Pmusicf);
回复
分享到:

使用道具 举报

发表于 2007-5-23 16:30 | 显示全部楼层
原帖由 chenhainiany 于 2007-5-23 16:27 发表
clear;
m=sqrt(-1);
delta=0.101043;
a1=-0.850848;
sample=32;      %number of sample spot
p=10;              %number of sample spot in coef method;
f1=0.05; f2=0.40; f3=0.42;
fstep=0.01;
fs ...



请先看看置顶贴:聚宝盆
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-12-12 16:36 , Processed in 0.076649 second(s), 18 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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