[yfsnbits]=wavread（‘G:\cfy\2017備課\數(shù)字信號處理-原理、實(shí)現(xiàn)及應(yīng)用\實(shí)驗(yàn)\語音實(shí)驗(yàn)\信息科學(xué)?（online-audio-converter.com）.wav‘）;
sound（yfsnbits）;

N=length（y）;
Y=fft（yN）;
subplot（211）;
plot（y）;title（‘原始信號波形‘）;
subplot（212）;
%?plot（1:Nfftshift（abs（Y））./max（max（abs（Y）））‘k‘‘linewidth‘2）;
%?title（‘原始信號頻譜‘）
plot（abs（Y））;title（‘原始信號頻譜‘）;
figurespecgram（y（:1）512fs100）;
title（‘語譜圖‘）


S2?=?y?*?4;?????????????????????????%對語音信號進(jìn)行放大
wavwrite（S2fsnbits‘louder.wav‘）;%生成放大后的語音信號
sound（S2）


fc=10000;????????????????????????????%載波頻率
y1=modulate（yfcfs‘fm‘）;??????????%對原語音信號調(diào)制
Y1=fft（y1N）;????????????????????%采樣點(diǎn)數(shù)40000
figure;
subplot（211）;plot（y1）;title（‘調(diào)制后信號波形‘）;
%?subplot（212）;plot（1:Nfftshift（abs（Y1））./max（max（abs（Y1）））‘r‘‘linewidth‘2）;
%?title（‘調(diào)制后信號頻譜‘）;
subplot（212）;plot