基于matlab关于多径效应的MATLAB仿真程序，共四个.rar

%this matlab code generates the frequency domain view of transmitted %signals and received signals for both cases. clear all; s=[ones(1,10),zeros(1,90)];%transmitted signal s_f=fft(s); x=s_f([1:50]); y=s_f([51:100]); signal_f=[y,x];%input spectrum dt=5/10;%each time interval is 0.01 micro sec ? df=1/(100*dt); f_f=df*([0:99]-50;%frequency vector an=[1 0.3 -0.8 0.5 -0.4 0.2];%amplitudes f=f_s; w=2*pi*f; tn_1=[0 1 2 3 4 5]l%arrival times for case 1 for i=1:6 h1(i,:)=an(i)*exp(-j*w*tn_1(i)); end h_1=sum(h1(:,1:end));%transfer function y_1=h_1.*signal_f;%output spectrum tn_2=[0 0.1 0.2 0.3 0.4 0.5];%arrival times for case2 for i=1:6 h2(i,:)=an(i)*exp(-j*w*tn_2(i)); end h_2=sum(h2(:,1:end));%transfer function y_2=h_2.*signal_f;%output spectrum figure(1); subplot(2,3,1); plot(f_s,abs(signal_f)); ylabel('magnitude'); title('I/P spectrum'); subplot(2,3,4); plot(f_s,angle(signal_f)); ylable('Phase'); xlable('Frequency(MHz)'); subplot(2,3,2); plot(f,abs(h_1); title('channel 1'); subplot(2,3,5); plot(f,angle(h_1)); xlable('Frequency(MHz)'); subplot(2,3,3); plot(f,abs(h_2); title('channel 2'); subplot(2,3,6); plot(f,angle(h_2)); xlable('Frequency(MHz)'); figure(2); subplot(2,3,1); plot(f_s,abs(signal_f)); ylabel('magnitude'); title('I/P spectrum'); subplot(2,3,4); plot(f_s,angle(signal_f)); ylable('Phase'); xlable('Frequency(MHz)'); subplot(2,3,2); plot(f,abs(y_1); title('O/P spectrum 1'); subplot(2,3,5); plot(f,angle(y_1)); xlable('Frequency(MHz)'); subplot(2,3,3); plot(f,abs(y_2); title('O/P spectrum 2'); subplot(2,3,6); plot(f,angle(y_2)); xlable('Frequency(MHz)');