%BPSK信号的产生
figure (74);
clear
df=0.15;
fs=20;%20MHz
fc=2;
rb=0.2;
ts=1/fs;
m=100;%m是码元数,自己随便定义
n=80;%n是每个码元取点数,即对应的正弦波4个周期取多少个点
a=sj(m,n);
%sj(m,n)的解释参见sj.m,sj(m,20)的作用相当于上面的for循环
t=0:1:m*n-1;
Fc1=sin(8*pi*t/n);%码元0发射的载波
Fc2=sin(8*pi*t/n+pi);%码元1发射的载波
b=1-a;
u0=a.*Fc2+b.*Fc1;
% for k=1:20 %对应程序结尾的end
%产生白噪声
snr=13;%S/N
snr_lin=10^(snr/10);%线性信噪比
signal_power=(norm(u0)^2)/length(u0);
noise_power=signal_power/snr_lin;
noise_std=sqrt(noise_power);
noise=noise_std*randn(1,length(u0));
u1=u0+noise;%加白噪声
u=u0;
ls=length(u);
%信号的一步小波分解
[cA1,cD1]=dwt(u,'db9');
%由小波系数构造一步近似与细节
A1=upcoef('a',cA1,'db9',1,ls);
D1=upcoef('d',cD1,'db9',1,ls);
%显示一步近似与细节
subplot(2,1,1);
plot(A1);title('近似A1');
subplot(2,1,2);
stem(abs(D1));title('细节D1');
grid
v=abs(D1);
% %求瞬时平方包络
% v=blpf(u);
%均值
miu=mean(v);
%方差
delt=var(v);
%求R参数
r=delt/(miu^2)
%求AWGN因子
e1=mean(v);
nfft=length(v);
j=0;
for i=1:nfft
if v(i)>e1
j=j+1;
ee(j)=v(i);
end
end
e2=mean(ee);
j=0;
for i=1:nfft
if v(i)>e2
j=j+1;
eee(j)=v(i);
end
end
e3=mean(eee);
a1=(e2-e1)/e1
a2=(e3-e2)/e1
% %通过逆小波变换重新构造信号
% A0=idwt(cA1,cD1,'db4',ls);
% %信号的多步小波分解
% [C,L]=wavedec(u,5,'db4');
% %提取近似与细节系数
% cA5=appcoef(C,L,'db4',5);
% cD5=detcoef(C,L,5);
% cD4=detcoef(C,L,4);
% cD3=detcoef(C,L,3);
% cD2=detcoef(C,L,2);
% cD1=detcoef(C,L,1);
% %重构第五层近似和各层细节
% A5=wrcoef('a',C,L,'db4',5);
% D5=wrcoef('d',C,L,'db4',5);
% D4=wrcoef('d',C,L,'db4',4);
% D3=wrcoef('d',C,L,'db4',3);
% D2=wrcoef('d',C,L,'db4',2);
% D1=wrcoef('d',C,L,'db4',1);
% %显示信号多层分解结果
% figure (3);
% subplot(3,2,1);plot(A5);title('近似A5');
% subplot(3,2,2);plot(D5);title('近似D5');
% subplot(3,2,3);plot(D4);title('近似D4');
% subplot(3,2,4);plot(D3);title('近似D3');
% subplot(3,2,5);plot(D2);title('近似D2');
% subplot(3,2,6);plot(D1);title('近似D1');
% %由5层小波分解结构重构原始信号
% A0=waverec(C,L,'db4');
% %显示重构信号与原始信号
% figure (4);
% subplot(2,1,1);plot(A0);title('重构信号A0');
% subplot(2,1,2);plot(u);title('原始信号');
% %信号除噪
% [thr,sorh,keepapp]=ddencmp('den','wv',u);
% de_noise=wdencmp('gbl',C,L,'db4',3,thr,sorh,keepapp);
% figure (5);
% subplot(2,1,1);plot(u);title('原始信号');
% subplot(2,1,2);plot(de_noise);title('除噪声后信号');
林当时
- 粉丝: 113
- 资源: 1万+
最新资源
- SparkSQL进阶操作相关数据
- java制作的小游戏,作为巩固java知识之用.zip
- Java语言写的围棋小游戏 半成品A Go game written in golang(Semi-finished).zip
- 基于Java-swing的俄罗斯方块游戏:源码+答辩文档+PPT.zip
- florr map详细版
- shiahdifhiahfiqefiwhfi weifwijfiwqufiqweefijeq0jfe
- registry-2.8.3<arm/amd>二进制文件
- Kotlin接口与抽象类详解及其应用
- 51单片机加减乘除计算器系统设计(proteus8.17,keil5),复制粘贴就可以运行
- lv_0_20241114231223.mp4
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈