[yfsbits]=wavread（‘C:\Users\ibm\Desktop\a.wav‘）;??????%語音信號的采集
n1=length?（y）?;????%求出語音信號的長度
Y=fft（yn1）;???????%傅里葉變換
f=8000000;?%載波頻率
M=size（y）;
t=0:1:M-1;
s=sin（f*2*3.14*t）;%構造載波序列
s=s‘;??%矩陣倒置
c=[ss];
A0=1;
A=ones（M）;%構造直流分量
y1=（y+A）.*c;%同步調制，給信號加密
n2=length（y1）;
Y1=fft（y1n2）;???????%傅里葉變換
wavwrite（y1?‘C:\Users\ibm\Desktop\加密.wav‘）;%?加密后的音頻信號名為“加密.wav”
y2=y1.*c;?%與載波相乘
n3=length（y2）;
Y2=fft（y2n3）;
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）;
[n11wn11]=buttord（wpws150‘s‘）;????%求低通濾波器的階數和截止頻率
[b11a11]=butter（n11wn11‘s‘）;????????%求S域的頻率響應的參數?
[num11den11]=bilinear（b11a110.5）;??%利用雙線性變換實現頻率響應S域到Z域的變換?
z11=filter（num11den11y2）;????????%低通濾波后的信號
n4=length（z11）;
m11=fft（z11n4）;??????????%求濾波后的信號頻譜
figure;
subplot（211）;
plot（z11）;
subplot（212）;
plot（abs（m11））;
[ab]=size（z11）;
B=ones（ab）;
y3=z11-0.5*B;???%減去直流分量
figure;
subplot（321）;
plot（y）;
title（‘原始信號波形‘）;
subplot（322）;
plot（abs（Y））;?title（‘原始信號頻譜‘）;
subplot（323）;
plot（y1）;
title（‘加密信號波形‘）;
subplot（324）;
plot（abs（Y1））;?
title（‘加密信號頻譜‘）;
subplot（325）;
plot（y3）;
title（‘解密信號波形‘）;
n4=length（y3）;
Y3=fft（y3n4）;%傅里葉變換
subplot（326）;
plot（abs（Y3））;
title（‘解密信號頻譜‘）;
wavwrite（y3?‘C:\Users\ibm\Desktop\解密.wav‘）;%?解密后的音頻信號名為“解密.wav”

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2015-08-03?20:37??基于matlab的語音加密器的設計\
?????文件?????2297900??2000-06-12?08:51??基于matlab的語音加密器的設計\a.wav
?????文件?????2297900??2000-06-12?11:17??基于matlab的語音加密器的設計\加密.wav
?????文件??????145920??2015-08-03?18:00??基于matlab的語音加密器的設計\基于matlab的語音加密器的設計?設計報告.doc
?????文件??????854016??2015-08-03?18:00??基于matlab的語音加密器的設計\基于matlab的語音加密器的設計.ppt
?????文件?????2297900??2000-06-12?11:18??基于matlab的語音加密器的設計\解密.wav
?????文件????????1534??2014-06-19?12:08??基于matlab的語音加密器的設計\語音加密.asv
?????文件????????1534??2014-06-19?12:13??基于matlab的語音加密器的設計\語音加密.m