多径瑞利分布冲激响应

rayleighchanresponse.mN=1000;ts=0.001;%采样间隔fs=1/ts;fd=10;%MaximumDopplershiftmustbelessthan1/(10*Ts),就是fdfs/10,对应的相干时间10*ts,原理如此tau=[00.010.20.3];%每个径相对于第1径的时延,小于相干时间可以合并，一般设置大于10*ts，较大时候出现波动PDB=[0-1-3-3];%每个径相对于第1径的增益dbtx=dirac(0:N-1);tx(1)=1;%产生单位冲击凼数tx2=randi([11],N,1);%发送序列chan_s=rayleighchan(ts,fd);%单径衰落信道chan_m=rayleighchan(ts,fd,tau,PDB);%多径衰落信道h_s=filter(chan_s,tx);%单径响应out_s=filter(chan_s,tx2);%单径输出序列h_m=filter(chan_m,tx);%多径响应out_m=filter(chan_m,tx2);%多径输出序列figure(1)subplot(211)plot(abs(h_s)+0.5,'-r');%幅度holdonplot(abs(h_s)/max(abs(h_s)),'b');holdon%plot(20*log10(10e-3+abs(h_s)),'-k');holdonplot(20*log10(0.02+abs(h_s)/max(abs(h_s))),'-k');%功率，多加legend('幅度绝对值','幅度相对值','单径相对功率值dB')axis([-20N-365]);title('冲激响应(单径)')xlabel('采样点')subplot(212)plot(abs(out_s),'-r');%幅度holdon%plot((abs(h_s)).^2,'-b');%holdonplot(20*log10(10e-4+abs(out_s)/max(abs(out_s))),'-k');%相对功率legend('幅度','相对功率dB')axis([-20N-605]);title('输出凼数（单径）')xlabel('采样点')figure(2)subplot(211)plot(abs(h_m)+0.5,'-r');holdonplot(abs(h_m)/max(abs(h_m)),'b');holdon%plot(20*log10(10e-3+abs(h_s)),'-k');holdonplot(20*log10(10e-4+0.02+abs(h_m)/max(abs(h_m))),'-k');legend('幅度绝对值','幅度相对值','相对功率值dB')axis([-20N-365]);title('冲激响应(多径)')xlabel('采样点')subplot(212)plot(abs(out_m),'-r');%幅度holdon%plot((abs(h_m)).^2,'-b');holdonplot(20*log10(abs(out_m)/max(abs(out_m))),'-k');%相对功率legend('幅度','相对功率dB')axis([-20N-605]);title('输出凼数（多径）')xlabel('采样点')%tau=[00.010.20.3];注意从0开始的时延%在0.01s处有个尖峰，表示第二经到来，在第0.01s/0.001=0.01/ts=10个点处%同理0.2s，峰值在第200点，0.3s峰值在300点。