基于matlab的語(yǔ)音加密器的設(shè)計(jì)

基于matlab的語(yǔ)音加密器的設(shè)計(jì)設(shè)計(jì)要求

選擇一段語(yǔ)音信號(hào)作為分析對(duì)象，或錄制一段自己的語(yǔ)音信號(hào)，對(duì)其進(jìn)行頻譜分析，然后對(duì)信號(hào)進(jìn)行比較分析。給出加密解密系統(tǒng)的框圖，并用matlab實(shí)現(xiàn)。設(shè)計(jì)思路

采用AM調(diào)制對(duì)語(yǔ)音信號(hào)調(diào)制解調(diào)實(shí)現(xiàn)語(yǔ)音信號(hào)的加密解密，在MATLAB中用命令語(yǔ)句實(shí)現(xiàn)功能。設(shè)計(jì)框圖如下：

LPFsinwtsinwtsinwtA0s（t）（t）MATLAB程序[y,fs,bits]=wavread('C:\Users\ibm\Desktop\a.wav');%語(yǔ)音信號(hào)的采集n1=length(y);%求出語(yǔ)音信號(hào)的長(zhǎng)度Y=fft(y,n1);%傅里葉變換f=8000000;%載波頻率M=size(y);t=0:1:M-1;s=sin(f*2*3.14*t);%構(gòu)造載波序列s=s';%矩陣倒置c=[s,s];MATLAB程序A0=1;A=ones(M);%構(gòu)造直流分量y1=(y+A).*c;%AM調(diào)制，給信號(hào)加密n2=length(y1);Y1=fft(y1,n2);%傅里葉變換wavwrite(y1,'C:\Users\ibm\Desktop\加密.wav');%加密后的音頻信號(hào)名為“加密.wav”MATLAB程序y2=y1.*c;%與載波相乘n3=length(y2);Y2=fft(y2,n3);Ft=8000;Fp=1000;Fs=1200;wp=2*pi*Fp/Ft;ws=2*pi*Fs/Ft;fp=2*Ft*tan(wp/2);fs=2*Fs*tan(wp/2);MATLAB程序[n11,wn11]=buttord(wp,ws,1,50,'s');%求低通濾波器的階數(shù)和截止頻率[b11,a11]=butter(n11,wn11,'s');%求S域的頻率響應(yīng)的參數(shù)[num11,den11]=bilinear(b11,a11,0.5);%利用雙線性變換實(shí)現(xiàn)頻率響應(yīng)S域到Z域的變換z11=filter(num11,den11,y2);%低通濾波后的信號(hào)n4=length(z11);m11=fft(z11,n4);%求濾波后的信號(hào)頻譜[a,b]=size(z11);B=ones(a,b);y3=z11-0.5*B;%減去直流分量MATLAB程序figure;subplot(3,2,1);plot(y);title('原始信號(hào)波形');subplot(3,2,2);plot(abs(Y));title('原始信號(hào)頻譜');subplot(3,2,3);plot(y1);title('加密信號(hào)波形');subplot(3,2,4);plot(abs(Y1));title('加密信號(hào)頻譜');MATLAB程序subplot(3,2,5);plot(y3);title('解密信號(hào)波形');n4=length(y3);Y3=fft(y3,n4);%傅里葉變換subplot(3,2,6);plot(abs(Y3));title('解密信號(hào)頻譜');wavwrite