WHT.zip_LFM WHT_LFM 检测_WHT_WHT MATLAB

% 用chirp函数产生LFM信号 % %第一路LFM信号 % t1=0:0.001:1; % N1 = length(t1); % y1 = chirp(t1,100,1,20); % y1 = hilbert(y1); % fft1=fftshift(fft(y1)); % % figure(1); % subplot(211); % plot(t1,real(y1)); % title('第一路LFM信号时域波形'); % xlabel('时间'); % ylabel('幅度'); % grid on; axis([0 1 -1 1]); % subplot(212); % plot(abs(fft1)); % title('第一路LFM信号频谱'); % xlabel('频率'); % ylabel('幅度'); % grid on; %第二路LFM t2=0:0.001:1; N2 = length(t2); y2 = chirp(t2,200,1,20); y2 = hilbert(y2); fft2=fftshift(fft(y2)); figure(2); subplot(211); plot(t2,real(y2)); % title('第二路LFM信号时域波形'); xlabel('时间'); ylabel('幅度'); grid on; subplot(212); plot(abs(fft2)); title('信号频谱'); xlabel('频率'); ylabel('幅度'); grid on; % % %两路信号分别加噪 % y1_n = awgn(y1,-3); % figure(3); % subplot(211); % plot(t1, real(y1_n)); % title('第一路信号加噪后时域波形'); % xlabel('时间');ylabel('幅度');grid on; % % fft3=fftshift(fft(y1_n)); % subplot(212); % plot(abs(fft3)); % title('第一路信号加噪后的频谱'); % xlabel('频率');ylabel('幅度');grid on;axis tight; % % y2_1n = awgn(y2,-10); % y2_2n=awgn(y2,-5); % y2_3n=awgn(y2,0); % y2_4n=awgn(y2,5); % % figure(4); % % subplot(211); % % plot(t2, real(y2_n)); % % title('第二路信号加噪后时域波形'); % % xlabel('时间');ylabel('幅度');grid on; % % % % fft4=fftshift(fft(y2_n)); % % subplot(212); % % plot(abs(fft4)); % % title('第二路信号加噪后的频谱'); % % xlabel('频率');ylabel('幅度');grid on;axis tight; % % % % % % % Wigner-Ville Distribution % % %第一路信号WVD % % [tfr1,t1,f1] = wv(y1); % % % f = f * (N1-1)/2/N1; % % t1 = t1 * 1/1000; % % [F1, T1] = meshgrid(f1, t1); % % figure(5); % % subplot(211); % % mesh(F1, T1, abs(tfr1)); % % title('第一路信号WVD变换三维图'); % % xlabel('频率');ylabel('时间');zlabel('幅度');grid on;axis tight; % % % %第二路信号WVD % % % % [tfr2,t2,f2] = wv(y2); % % % % % f = f * (N1-1)/2/N1; % % % % t2 = t2 * 1/1000; % % % % [F2, T2] = meshgrid(f2, t2); % % % % % % % % figure; % % % % subplot(121); % % % % mesh(F2, T2, abs(tfr2)); % % % % % title('第二路信号WVD变换三维图'); % % % % xlabel('频率');ylabel('时间');zlabel('幅度');grid on;axis tight; % % % % figure(6); % % contour(f1,t1,abs(tfr1)); % % title('第一路信号WVD'); % % xlabel('频率'); % % ylabel('时间');grid on; % % % % % figure(1); % % % subplot(122); % % % contour(f2,t2,abs(tfr2)); % % % % title('第二路信号WVD'); % % % xlabel('频率'); % % % ylabel('时间');grid on; % % % % figure(8); % % [tfr3,t3,f3] = wv(y1_n); % % t3 = t3 * 1/1000; % % [F3, T3] = meshgrid(f3, t3); % % mesh(F3, T3, abs(tfr3)); % % title('第一路信号加噪WVD变换三维图'); % % xlabel('频率');ylabel('时间');zlabel('幅度');grid on;axis tight; % % % figure(1); % [tfr1,t1,f1] = wv(y2_1n); % t1 = t1 * 1/1000; % [F1, T1] = meshgrid(f1, t1); % subplot(221); % mesh(F1, T1, abs(tfr1)); % title('snr=-10'); % xlabel('f');ylabel('t');zlabel('A');grid on;axis tight; % % [tfr2,t2,f2] = wv(y2_2n); % t2 = t2 * 1/1000; % [F2, T2] = meshgrid(f2, t2); % subplot(222); % mesh(F2, T2, abs(tfr2)); % title('snr=-5'); % xlabel('f');ylabel('t');zlabel('A');grid on;axis tight; % % [tfr3,t3,f3] = wv(y2_3n); % t3= t3 * 1/1000; % [F3, T3] = meshgrid(f3, t3); % subplot(223); % mesh(F3, T3, abs(tfr3)); % title('snr=0'); % xlabel('f');ylabel('t');zlabel('A');grid on;axis tight; % % % [tfr4,t4,f4] = wv(y2_4n); % t4= t4 * 1/1000; % [F4, T4] = meshgrid(f4, t4); % subplot(224); % mesh(F4, T4, abs(tfr4)); % title('snr=5'); % xlabel('f');ylabel('t');zlabel('A');grid on;axis tight; % % % % % % % % % figure(2); % subplot(221); % contour(f1,t1,abs(tfr1)); % title('snr=-10'); % xlabel('f'); % ylabel('t');grid on; % % subplot(222); % contour(f2,t2,abs(tfr2)); % title('snr=-5'); % xlabel('f'); % ylabel('t');grid on; % % subplot(223); % contour(f3,t3,abs(tfr3)); % title('snr=0'); % xlabel('f'); % ylabel('t');grid on; % % subplot(224); % contour(f4,t4,abs(tfr4)); % title('snr=5'); % xlabel('f'); % ylabel('t');grid on; % % % % figure(11); % % contour(f4,t4,abs(tfr4)); % % title('第二路信号加噪后WVD'); % % xlabel('频率'); % % ylabel('时间');grid on; % % % % %两个信号相加 % % % % y3=y1+y2; % % t=0:0.001:1; % % fft5=fftshift(fft(y3)); % % figure(12); % % subplot(211); % % plot(t,real(y3)); % % title('合成信号时域波形'); % % xlabel('时间'); % % ylabel('幅度'); % % grid on; axis([0 1 -1 1]); % % subplot(212); % % plot(abs(fft5)); % % title('合成信号频谱'); % % xlabel('频率'); % % ylabel('幅度'); % % grid on; % % % % %合成信号WVD % % % % figure(13); % % [tfr6,t6,f6] = wv(y3); % % t6 = t6 * 1/1000; % % [F6, T6] = meshgrid(f6, t6); % % mesh(F6, T6, abs(tfr6)); % % title('合成信号加噪WVD变换三维图'); % % xlabel('频率');ylabel('时间');zlabel('幅度');grid on;axis tight; % % % % figure(14); % % contour(f6,t6,abs(tfr6)); % % title('合成信号WVD'); % % xlabel('频率'); % % ylabel('时间');grid on; % % % % %合成信号加噪 % % % % y3_n = awgn(y3,-3); % % figure(15); % % subplot(211); % % plot(t, real(y3_n)); % % title('合成信号加噪后时域波形'); % % xlabel('时间');ylabel('幅度');grid on; % % % % fft6=fftshift(fft(y3_n)); % % subplot(212); % % plot(abs(fft6)); % % title('合成信号加噪后的频谱'); % % xlabel('频率');ylabel('幅度');grid on;axis tight; % % % % %合成信号与合成信号加噪后的WVD % % % % figure(16); % % [tfr5,t5,f5] = wv(y3_n); % % t5 = t5 * 1/1000; % % [F5, T5] = meshgrid(f5, t5); % % mesh(F5, T5, abs(tfr5)); % % title('合成信号加噪WVD变换三维图'); % % xlabel('频率');ylabel('时间');zlabel('幅度');grid on;axis tight; % % % % figure(17); % % contour(f5,t5,abs(tfr5)); % % title('合成信号加噪WVD'); % % xlabel('频率'); % % ylabel('时间');grid on; % % % % % % % % % Hough Transform % % % % %第一路LFM信号WHT % % % % [ht1, rho1, theta1] = hough(tfr1, f1, t1); % % figure(18); % % mesh(theta1*180/pi, rho1, abs(ht1)); % % title('第一路信号WHT'); % % xlabel('角度'); ylabel('极坐标半径'); zlabel('幅度');grid on; % % axis tight; % % % % [ht2, rho2, theta2] = hough(tfr2, f2, t2); % % figure(19); % % mesh(theta2*180/pi, rho2, abs(ht2)); % % title('第二路信号WHT'); % % xlabel('角度'); ylabel('极坐标半径'); zlabel('幅度');grid on; % % axis tight; % % % % [ht3, rho3, theta3] = hough(tfr3, f3, t3); % % figure(19); % % mesh(theta3*180/pi, rho3, abs(ht3)); % % title('第一路信号加噪WHT'); % % xlabel('角度'); ylabel('极坐标半径'); zlabel('幅度');grid on; % % axis tight; % % % % [ht4, rho4, theta4] = hough(tfr4, f4, t4); % % figure(20); % % mesh(theta4*180/pi, rho4, abs(ht4)); % % title('第二路信号加噪WHT'); % % xlabel('角度'); ylabel('极坐标半径'); zlabel('幅度');grid on; % % axis tight; % % % % [ht6, rho6, theta6] = hough(tfr6, f6, t6); % % figure(21); % % mesh(theta6*180/pi, rho6, abs(ht6)); % % title('合成信号WHT'); % % xlabel('角度'); ylabel('极坐标半径'); zlabel('幅度');grid on; % % axis tight; % % % % [ht5, rho5, theta5] = hough(tfr5, f5, t5); % % figure(22); % % mesh(theta5*180/pi, rho5, abs(ht5)); % % title('合成信号加噪WHT'); % % xlabel('角度'); ylabel('极坐标半径'); zlabel('幅度');grid on; % % axis tight; % % % % % figure; image(theta*(180/pi), rho, abs(ht)); % % % xlabel('theta'); ylabel('rho'); % % %