資源簡介
包含線性調(diào)頻信號的參數(shù)估計,兩個:初始頻率和調(diào)制斜率。

代碼片段和文件信息
clc;
close?all;
clear;
%?T=10e-6;
%?B=30e6;
%?K=B/T;
%?fs=2*B;
%?ts=1/fs;
%?N=T/ts;
%?t=linspace(-T/2T/2N);
%?st=exp(j*pi*K*t.^2);
%?t=-1*ts/2:1/fs:ts/2;
%?S=sqrt(ts/fs);
%?dx=sqrt(ts*fs);
%?tt=-1*dx/2:1/dx:dx/2;%歸一化之后的數(shù)據(jù)
fs=256;T=2;
f0=10;k=10;
S=sqrt(T/fs);
dx=sqrt(T*fs);
ts=1/fs;
point=ts*fs+1;
t=-1:0.01:1;
tt=-1*dx/2:1/dx:dx/2;
f1=20;k1=20;
x=exp(1i*2*pi*f0*tt*S+1i*pi*k*tt.*tt*S*S)+exp(1i*2*pi*f1*tt*S+1i*pi*k1*tt.*tt*S*S);%exp(j*2*pi*((f0/fs)*t+0.5*k*t.*t));
p=0.9:0.01:1.8;
u=-dx/2:1/dx:dx/2;
Pp3=zeros(length(u)length(p));
for?k=1:1:length(p)
????p3=frft(xp(k));
????%vx3=frft(vxa);
????%vy3=frft(vya);
????Pp3(:k)=p3.*conj(p3)/point;
end
plot(pPp3);
figure;
surf(puPp3);
[CXCI]=max(max(Pp3));
p(CI);
Pp31=max(max(Pp3));
[x1?y1]?=?find(Pp3==Pp31);
u?=?u(x1);
p?=?p(y1);
[ab]=jisuan(pufsT)
?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件?????????888??2019-05-30?11:41??工程\ce.m
?????文件????????1638??2019-05-29?09:24??工程\frft.m
?????文件?????????268??2019-05-30?11:04??工程\jisuan.m
?????目錄???????????0??2019-10-11?08:57??工程\
- 上一篇:哈工大同義詞詞林.zip
- 下一篇:ACM比賽必帶模板
評論
共有 條評論