多径传播matlab仿真·信号分析处理

多径传播matlab仿真

writer： M宝可梦

输入信号

      s(t)   =   Acos(wc*t)

经过多条路径时变的衰落和时延得到接收信号

接收信号：

  r(t) = Σ ai(t)*cos{wc*[t - Ti(t)]}

% 根据模拟通道数为100，所以求和次数为100次

  = Σ ai(t)*cos[wc*t + PHi(t)]      % phase 相位; PHi(t) = -wc * Ti(t)% 仿真起始点   = Σ ai(t)*cos[PHi(t)]*cos(wc*t) - Σai(t)*sin[PHi(t)]*sin(wc*t)= X(t)*cos(wc*t) - Y(t)*sin(wc*t)= V(t)cos[wc*t + PH(t)]

其中：

         X(t) = Σai(t)*cos[PHi(t)]  Y(t) = Σai(t)*sin[PHi(t)]  V(t) = sqrt[X^2 + Y^2]PH(t) = arctan(Y/X)

仿真思路：
1.衰减ai(t),以及时延PHi(t) 的模拟：
·幅值衰减由随机数产生，符合[0,1]上的均匀分布
·PHi(t)时延产生方式同幅值衰减
·为了保持与原理上的一致性，引入变量t影响两个参数
2.wc 为超参数##########未使用
3.输出V(T),PH(t)：

代码说明

·为方便比较最终接收端的信号分布，做出了瑞利分布图象
·sum_signal函数为接受所有通道的信号之叠加

仿真结果如图

% 参考：
https://blog.csdn.net/weixin_44606638/article/details/106762784
https://blog.csdn.net/qq_42691315/article/details/115015058

clear;
clc;
% 参数说明：
% signal 第一列为幅值，第二列为相位
% f_width，x_width 为幅值的概率和数值
% f_phase,x_phase  为相位的概率和数值for i = 1:10000 %仿真次数[signal(i,1),signal(i,2)] = sum_signal(100);
end[f_width,x_width]=ksdensity(signal(:,1));
[f_phase,x_phase]=ksdensity(signal(:,2));% 将相位划归到[-pi, pi]区间上
for i = 1: 100delta = max(x_phase(1,:)) - min(x_phase(1,:));x_phase(1,i) = (x_phase(1,i) -min(x_phase(1,:)))/delta;x_phase(1,i) = x_phase(1,i)*2*pi-pi;
end% rayleigh分布的pdf
cs2=var(signal(:,1))*2;  %方差为实部或虚部的方差
r=linspace(0,50,1000);
fx2=r./(cs2).*exp(-r.^2/2/(cs2));% 作图
figure(1);
subplot(121);
plot(x_width,f_width);
hold on;plot(r,fx2,'r:');hold off;
ylabel('幅度分布');xlabel('幅度');title('幅度分布');axis([-2 20 0 0.2]);
subplot(122);
plot(x_phase,f_phase);
hold on;plot([-pi pi],[1/2/pi 1/2/pi],'r:');hold off;
ylabel('相位分布');xlabel('相位/rad');title('相位分布');axis([-pi-0.5 pi+0.5 0 0.2]);function [width,phase] = sum_signal(channel_num)% 通过叠加的形式求得最终接收信号,此为一次仿真结果
% 参数信息、
% input：
%       channel_num    ： 路径数目
% output:
%       phase :  接收到信号相位移动
%       width ： 接收到信号的幅值X = 0;
Y = 0;
for i = 1:channel_numa(i) = rand(1,1);           % 幅值衰减随机赋值为[0,1]    ph(i) = 2*pi*rand(1,1) - pi;% 将时延相位划归到 [-pi,pi]X = X + a(i)*cos(ph(i));Y = Y + a(i)*sin(ph(i));
end
width = sqrt(X^2+Y^2);
phase = atan2(Y,X);
end

多径传播matlab仿真·信号分析处理相关推荐

通信原理实验多径传播 MATLAB仿真
通信原理实验多径传播 MATLAB仿真一.实验目的通过实验掌握多径传播.信道的频率选择性.相干带宽等概念,理解多径信道对信号传输的影响. 二.实验原理多径信道指信号传输的路径不止一条,接收端同 ...
《MATLAB语音信号分析与合成（第二版）》：第4章语音信号的线性预测分析
<MATLAB语音信号分析与合成(第二版)>:第4章语音信号的线性预测分析前言 1. 数据与函数路径设置 2. MATLAB仿真一:LPC的频谱分析 3. MATLAB仿真二:普通预测 ...
《MATLAB语音信号分析与合成（第二版）》：第8章基音周期的估算方法
<MATLAB语音信号分析与合成(第二版)>:第8章基音周期的估算方法前言 1. 数据与函数路径设置 2. MATLAB仿真一:基音周期提取的预处理 3. MATLAB仿真二:倒谱法的 ...
《MATLAB语音信号分析与合成（第二版）》：第7章语音信号的减噪
<MATLAB语音信号分析与合成(第二版)>:第7章语音信号的减噪前言 1. 数据与函数路径设置 2. MATLAB仿真一:自适应滤波器LMS滤波减噪一 3. MATLAB仿真二:自适 ...
《MATLAB语音信号分析与合成（第二版）》：第10章语音信号的合成算法
<MATLAB语音信号分析与合成(第二版)>:第10章语音信号的合成算法前言 1. 数据与函数路径设置 2. MATLAB仿真一:重叠相加法语音合成 3. MATLAB仿真二:重叠存储 ...
matlab对声场仿真,基于声线模型的水声传播MATLAB仿真.pdf
基于声线模型的水声传播MATLAB仿真.pdf 科技广场 2007.9 基于声线模型的水声传播MATLAB仿真熊光耀杨琴 Ciong GuangyaoYang Qin (江西中医学院计算机系,江 ...
《MATLAB语音信号分析与合成（第二版）》：第3章语音信号在其他变换域中的分析技术和特性
<MATLAB语音信号分析与合成(第二版)>:第3章语音信号在其他变换域中的分析技术和特性前言 1. 数据与函数路径设置 2. MATLAB仿真一:信号倒谱图.声道冲激响频谱.声门激励 ...
降压斩波占空比matlab,升降压斩波电路的MATLAB 仿真及分析.docx
升降压斩波电路的 MATLAB仿真及分析院系: 班级: 姓名: 学号: 时间: 地点: 一．问题介绍设输入电压Uo=50V 纹波电压小于0.02% 脉冲周期T=1e-4s 负载R=450Ω 电感 ...
三相变压器matlab仿真与分析,三相变压器励磁涌流的MATLAB仿真与分析.doc
三相变压器励磁涌流的MATLAB仿真与分析三相变压器励磁涌流的MATLAB仿真与分析摘要:在Simulink仿真平台上建立三相变压器励磁涌流仿真模型,对变压器空载合闸投入和外部故障切除后两种情况进 ...