基於matlab的tcm編解碼
⑴ pcm 的matlab實現未編碼、編碼、解碼波形和不同幅度的量化信噪比,我的程序有問題,請大神幫幫忙
我修改了一下,能夠運行了,不知道是不是你想要的結果?
clear all;
close all;
t=0:0.01:10; %定義時間抽樣點
%vm1=-70:1:0; %輸入的信號幅度的db值
vm1=linspace(-70,0,1001);
vm=10.^(vm1/20); %輸入信號幅度
figure(1);
for k=1:length(vm)
for m=1:2
%x=vm*sin(2*pi*t+2*pi*rand(t)); %輸入語音信號
x=vm.*sin(2*pi*t); %輸入語音信號
v=1;
xx=x/v;
sxx=floor(xx*4096);
y=pcm_encode(sxx); %PCM編碼
yy=pcm_decode(y,v); %PCM解碼後信號幅值
nq(m)=sum((x-yy).*(x-yy))/length(x);%噪音功率
sq(m)=mean(yy.^2); %信號均值
snr(m)=(sq(m)/nq(m)); %信噪比
drawnow
subplot(211)
plot(t,x);
title('采樣序列'); %畫出采樣序列的圖形
subplot(212)
plot(t,yy);
title('解碼序列'); %畫出PCM解碼後的序列圖
end
snrq(k)=10*log10(mean(snr)); %量化信噪比
end
figure(2)
plot(vm1,snrq);
axis([-60 0 0 60]); %X軸范圍是(-60,0)Y軸范圍是(0,60)
grid;
function out=pcm_decode(in,v);
n=length(in);
in=reshape(in',8,n/8)'; %將in值變換成8行
slot(1)=0;
slot(2)=16;
slot(3)=32;
slot(4)=64;
slot(5)=128;
slot(6)=256;
slot(7)=512;
slot(8)=1024;
step(1)=1;
step(2)=1;
step(3)=2;
step(4)=4;
step(5)=8;
step(6)=16;
step(7)=32;
step(8)=64;
for i=1:n/8
ss=2*in(i,1)-1;
tmp=in(i,2)*4+in(i,3)*2+in(i,4)+1;
st=slot(tmp);
dt=(in(i,5)*8+in(i,6)*4+in(i,7)*2+in(i,8))*step(tmp)+0.5*step(tmp);
out(i)=ss*(st+dt)/4096*v; %量化輸出值
end
function [out]=pcm_encode(x); %定義A率13折線壓縮特性
n=length(x);
for i=1:n
if x(i)>0
out(i,1)=1; %代表正值
else
out(i,1)=0; %代表負值
end
if abs(x(i))>=0&&abs(x(i))<16
out(i,2)=0;out(i,3)=0;out(i,4)=0;step=1;st=0;
elseif 16<=abs(x(i))&&abs(x(i))<32
out(i,2)=0;out(i,3)=0;out(i,4)=1;step=1;st=16;
elseif 32<=abs(x(i))&&abs(x(i))<64
out(i,2)=0;out(i,3)=1;out(i,4)=0;step=2;st=32;
elseif 64<=abs(x(i))&&abs(x(i))<128
out(i,2)=0;out(i,3)=1;out(i,4)=1;step=4;st=64;
elseif 128<=abs(x(i))&&abs(x(i))<256
out(i,2)=1;out(i,3)=0;out(i,4)=0;step=8;st=128;
elseif 256<=abs(x(i))&&abs(x(i))<512
out(i,2)=1;out(i,3)=0;out(i,4)=1;step=16;st=256;
elseif 512<=abs(x(i))&&abs(x(i))<1024
out(i,2)=1;out(i,3)=1;out(i,4)=0;step=32;st=512;
elseif 1024<=abs(x(i))&&abs(x(i))<2048
out(i,2)=1;out(i,3)=1;out(i,4)=1;step=64;st=1024; %由抽樣值定義段落編碼
end
if(abs(x(i))>=2048)
out(i,2:8)=[1 1 1 1 1 1 1];step=128;st=2048;
else
tmp=floor((abs(x(i))-st)/step);
if tmp<0
a=1;
end
t=dec2bin(tmp,4)-48; %dec2bin函數表示輸出ASCII字元串值,48表示0
out(i,5:8)=t(1:4); %輸出段內碼
end
end
out=reshape(out',1,8*n); %將out值變換成1行8n列
⑵ MATLAB: 數字通信系統信道編碼 AMI 編譯碼
程序如下,現在原始序列長度20的隨機0,1串,要變自己改。
clc;
clear;
source = randint(1,20);
%%%%%%%%%%%% Encode %%%%%%%%%%%%%
perbit = -1;
for i=1:length(source);
if source(i)==1
encoded(i) = (-1)*perbit;
perbit = encoded(i);
else
encoded(i) = source(i);
end
end
%%%%%%%%%%%% Decode %%%%%%%%%%%%
for i=1:length(source);
if encoded(i)~=0
decoded(i) = 1;
else
encoded(i) = 0;
end
end
source
encoded
decoded
⑶ 求TCM網格編碼的matlab程序代碼
不懂,關注。
⑷ 用matlab viterbi演算法怎麼求傳輸
[mlse_Eq.zip] - 含有已知信道下BPSK,QPSK,8PSK,16QAM的mlse均衡,MMSE(頻域)均衡SER,BER性能的對比。
[wanzhengchengxu.rar] - 信號在信道進行傳輸時,經常會出現串擾,採用基於最大似然序列估計(mlse)的電子色散均衡器的方法,克服光纖通信中的由各種色散引起的碼間干擾。 研究基於最大似然序列估計(MLSE)的均衡器,採用維特比演算法實現,用matlab模擬出採用MLSE後的性能指標(眼圖和誤碼率)的提升。
[TCMchengxu.rar] - matlab實現8PSK TCM網格編碼調制,包含維特比解碼和與未編碼的QPSK誤碼性能比較
[OFDM_code.rar] - 本代碼是關於OFDM通信系統的模擬源程序,其中信道編碼,QPSK調制,QPSK解調,FFT,IFFT,viterbi編解碼,誤比特統計,濾波等模塊。
⑸ 想問下關於循環碼解碼的程序用matlab怎麼編寫啊急求啊
推薦本是你看下吧,看完自己基本就能搞定了
<matlab 通信模擬開發手冊>
⑹ 如何使用matlab通訊工具箱的函數實現RS編解碼
rsenc RS編碼器
rsdecof 將RS編碼的ASCII文件解碼
rsencof 對一個ASCII文件進行RS編碼
主要使用這幾個函數,你仔細看看
祝你愉快,滿意請採納。
⑺ 求信道極化碼polar碼的MATLAB編解碼程序
親,這個沒法給,因為採用不同解碼方式代碼都有不同的
⑻ MATLAB編程完成PCM編、解碼(交源代碼及說明)
出售與定製Matlab畢業設計和課程設計或者文獻翻譯,成品部分羅列http://975941553.qzone.qq.com/infocenter?ADUIN=975941553&ADSESSION=1241335843&ADTAG=CLIENT.QQ.1855_QQUrlReportBlankZone.0
你可以在需要的那個帖子下留下聯系方式或者直接加該QQ為好友祥談定製您需要的畢業設計.請期待更多的作品,希望得到您的認可,滿足您的畢業要求.