four_case_doppler.m

多径多普勒信道模型

clear;
an=[1,0.3,-0.8,0.5,-0.4,0.2];
tn=[0,1,2,3,4,5;0,1,2,3,4,5;0,0.1,0.2,0.3,0.4,0.5;0,0.1,0.2,0.3,0.4,0.5];
wn=[0,2,10,6,8,4;0,20,100,60,80,40;0,2,10,6,8,4;0,20,100,60,80,40];
f=-1:0.01:1;
w=2*pi*f;
for k=1:4;
for i=1:6;
h1(i,:)=an(i)*exp(-j*w*tn(k,i))*exp(j*wn(k,i)*0);
h2(i,:)=an(i)*exp(-j*w*tn(k,i))*exp(j*wn(k,i)*0.02);
end
h1_1(k,:)=sum(h1(:,1:end));
h1_2(k,:)=sum(h2(:,1:end));
end
subplot(4,2,1);plot(f,abs(h1_1(1,:)));ylabel('case 1');ylim([0 4]);
hold on
plot(f,abs(h1_2(1,:)) ,'r');title('amplitude');
%legend('t_0=0msec','t_0=20msec')
hold off
subplot(4,2,2);plot(f,angle(h1_1(1,:))); ylim([-4 4]);
hold on
plot(f,angle(h1_2(1,:)) ,'r'); title('phase'); %
hold off
subplot(4,2,3);plot(f,abs(h1_1(2,:)));ylabel('case 2');ylim([0 4]);
hold on
plot(f,abs(h1_2(2,:)) ,'r');
hold off
subplot(4,2,4);plot(f,angle(h1_1(2,:)));ylim([-4 4]);
hold on
plot(f,angle(h1_2(2,:)) ,'r');
hold off
subplot(4,2,5);plot(f,abs(h1_1(3,:)));;ylabel('case 3');ylim([0 4]);
hold on
plot(f,abs(h1_2(3,:)) ,'r');
hold off
subplot(4,2,6);plot(f,angle(h1_1(3,:))); ylim([-4 4]);
hold on
plot(f,angle(h1_2(3,:)) ,'r');
hold off
subplot(4,2,7);plot(f,abs(h1_1(4,:)));xlabel('Frequency(MHz)');
ylabel('case 4');ylim([0 4]);
hold on
plot(f,abs(h1_2(4,:)) ,'r');
hold off
subplot(4,2,8);plot(f,angle(h1_1(4,:)));xlabel('Frequency(MHz)');
ylim([-4 4]);
hold on
plot(f,angle(h1_2(4,:)) ,'r');
hold off;