📄 b.m
字号:
clear all;
% A=zeros(1,800);
% A(1,1:50)=1;A(351:400)=1;
% plot(abs(fftshift((fft(x)))))
% t=0:0.002:4;
% x=zeros(1,length(t));
% for m=1:52
% x=x+sin(2*pi*(m*312500+5000000000)*t);
% end
% plot(x);
% X=fft(x,64);
% P=X.*conj(X)/64;
% f=(64*312500+5000000000)*(1:64)/64;
% plot(f,abs(X));
t=0:0.001:0.6;
x=zeros(1,length(t));
for m=1:52
x=x+sin(2*pi*((m-1)*50+50)*t);
end
f_max=400;
m_s=1;
t_0=0;
T_s=0.2E-6;
[C1,F1,TH1,C2,F2,TH2,F01,F02,RHO,F_RHO,q_l,T]=F_S_K_p(9,'tu',f_max);
[y_t,T,t_0]=F_S_K(x,f_max,m_s,T,t_0,q_l,C1,F1,TH1,C2,F2,TH2,F01,F02,RHO,F_RHO,T_s);
y=y_t+1.5*randn(1,length(t));
Y=fft(y,512);
P=Y.*conj(Y)/512;
f=1000*(0:255)/512;
plot(f,P(1:256))
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -