int main ()
{
char c
while ((c = getchar()) != '\n')
{
if(c==' ')
{
printf("*")
continue
}
c-=64
if (c!=-32)
printf ("%d ",c)
}
printf ("凳凳销.")
return 0
}
试试看符不枣游符粗缓合你的要求
% clc% clear all
% 产生信号
load handel
x=y(1:20000)%取前20000个采样点
sound(x,Fs)
%PCM编码
x1=x/0.8.*2048
yy=pcm_encode(x1)
figure
subplot(2,1,1)
stem(yy(1:80),'.')
title('PCM编码后的波形')
%加噪声
snr=10
sp=mean(yy.^2)
attn=sp./ 10^(snr/10)
attn = sqrt(attn)
noise=randn(1,length(yy)).*attn
np=mean(noise.^2)
snr1=10*log10(sp/np)
data=yy+noise
% data=yy%不加噪声
subplot(2,1,2)
stem(data(1:80),'羡厅和.r')
title('PCM加噪声后波形')
%译码
demodata=data>0.5
zz=pcm_decode(demodata,0.8)
figure
subplot(2,1,1)
plot(x)
title('兄盯原始语音信号')
subplot(2,1,2)
plot(zz)
title('译码后的语音信号')
sound(zz,Fs)
figure
plot(x,'b')
hold on
plot(zz,'r')
legend('编译码前的语音信号','编译码后的语音信号')
title('编,译码前后的语音信号')
function y = pcm_encode( x )
y=zeros(length(x),8) %存储矩阵(全零)
z=sign(x) %判断x的正负
x=abs(x)%取绝伏仿对值
%%段落码判断段区间的取值范围为前开后闭区间
for k=1:length(x)
%符号位的判断
if z(k)>0
y(k,1)=1
elseif z(k)<0
y(k,1)=0
end
if x(k)>128 &x(k)<=2048 %在第五段与第八段之间,段位码第一位都为“1”
y(k,2)=1
end
if (x(k)>32 &x(k)<=128) || (x(k)>512 &x(k)<=2048)
y(k,3)=1 %在第三四七八段内,段位码第二位为“1”
end
if (x(k)>16&x(k)<=32)||(x(k)>64&x(k)<=128)||(x(k)>256&x(k)<=512)||(x(k)>1024&x(k)<=2048)
y(k,4)=1 %在二四六八段内,段位码第三位为“1”
end
end
%段内码判断程序
N=zeros(1,length(x))
for k=1:length(x)
N(k)=y(k,2)*4+y(k,3)*2+y(k,4)+1 %找到x位于第几段
end
a=[0,16,32,64,128,256,512,1024] %量化间隔
b=[1,1,2,4,8,16,32,64]%除以16,得到每段的最小量化间隔
for m=1:length(x)
q=ceil((x(m)-a(N(m)))/b(N(m))) %求出在段内的位置
if q==0
y(m,(5:8))=[0,0,0,0] %如果输入为零则输出“0”
else k=num2str(dec2bin(q-1,4)) %编码段内码为二进制
y(m,5)=str2num(k(1))
y(m,6)=str2num(k(2))
y(m,7)=str2num(k(3))
y(m,8)=str2num(k(4))
end
end
%将N行8列矩阵转换为1行8*N列的矩阵
y=y'
y=reshape(y,1,length(x)*8)
end
function x=pcm_decode(y,max)
%将1行8*N列的矩阵转换为N行8列矩阵
y=reshape(y,8,length(y)/8)
y=y'
%PCM译码
n=size(y,1) %求出输入码组的个数
a=[0,16,32,64,128,256,512,1024] %段落起点值
b=[1,1,2,4,8,16,32,64] %每段的最小量化间隔
for k=1:n
t1=y(k,1)%取符号
t2=y(k,2)*4+y(k,3)*2+y(k,4)+1%判断段落位置
t3=y(k,5)*8+y(k,6)*4+y(k,7)*2+y(k,8) %判断段内位置
if t3==0 %段内码为零时
m(k)=(a(t2)+1+0.5*b(t2))/2048*max
else
m(k)=(a(t2)+b(t2)*t3+0.5*b(t2))/2048*max%还原出量化后的电平值
end
%判断符号位
if t1==0
x(k)=-m(k)
else
x(k)=m(k)
end
end
end
首先要世带袭纠正一下,相对多的输入转化成为相对少的输出,一般叫编码器;相对少的输入转化成为相对多的输出,一般叫译码器。所以,确切地说你要做的应该是8-3编码器搜兄(还是3-8译码器)。//8-3编码器verilog代码
module enc_83(datain, dataout, en)
input [7:0] datain
input en
output [2:0] dataout
reg [2:0] dataout
always @*
if (~en)
dataout = 3'b0
else case(datain)
8'b0000_0001: dataout = 3'b000
8'b0000_0010: dataout = 3'b001
8'b0000_0100: dataout = 3'b010
8'b0000_1000: dataout = 3'b011
8'b0001_0000: dataout = 3'b100
8'b0010_0000: dataout = 3'b101
8'b0100_0000: dataout = 3'b110
8'b1000_0000: dataout = 3'b111
default: dataout = 3'b000
endcase
endmodule
//关于上述代码的说明
1) 上述实现的是8-3普通编行昌码器,即同一时刻输入保证只有一位有效;
2) 上述代码没有输出编码有效的标识位,如果需要的话可自行添加。
欢迎分享,转载请注明来源:内存溢出
评论列表(0条)