Fs=1000; T=1/Fs;
Tp=0.06; N=Tp*Fs; w=99*pi;
n=1:N;
Xn=sin(w*n*T+10*pi/180)+0.5*sin(3*w*n*T+20*pi/180)+0.5*sin(5*w*n*T+40*pi/180)+0.4*sin(7*w*n*T+60*pi/180)+0.3*sin(9*w*n*T+80*pi/180)+0.2*sin(11*w*n*T+90*pi/180);
Xn=Xn/max(abs(Xn));%请问这句话什么意思,干什么的
wn=boxcar(N);
Xn1=Xn.*wn';
Xk=fft(Xn1,4096);
fk=Fs*(0:4095);
plot(fk,abs(Xk)/max(abs(Xk)));
xlabel('Hz');
ylabel('幅值');
title('矩形窗,Tp=0.06,Fs=1000')
请问注释处是干什么的,为了求什么?
Tp=0.06; N=Tp*Fs; w=99*pi;
n=1:N;
Xn=sin(w*n*T+10*pi/180)+0.5*sin(3*w*n*T+20*pi/180)+0.5*sin(5*w*n*T+40*pi/180)+0.4*sin(7*w*n*T+60*pi/180)+0.3*sin(9*w*n*T+80*pi/180)+0.2*sin(11*w*n*T+90*pi/180);
Xn=Xn/max(abs(Xn));%请问这句话什么意思,干什么的
wn=boxcar(N);
Xn1=Xn.*wn';
Xk=fft(Xn1,4096);
fk=Fs*(0:4095);
plot(fk,abs(Xk)/max(abs(Xk)));
xlabel('Hz');
ylabel('幅值');
title('矩形窗,Tp=0.06,Fs=1000')
请问注释处是干什么的,为了求什么?