clc
clear?all
%***************1.正弦波****************%
fs?=100;????????????????????????????????????????????%設(shè)定采樣頻率
N?=128;
n?=0:N?-1;
t?=?n/?fs;
f0?=10;?????????????????????????????????????????????%設(shè)定正弦信號頻率
%生成正弦信號
x?=?sin（2*pi*f0*t）;
figure（1）;
subplot（231）;
plot（tx）;?%作正弦信號的時域波形
xlabel（‘時間/?s‘）;
ylabel（‘幅值‘）;
title（‘時域波形‘）;
grid;
%進行FFT?變換并做頻譜圖
y?=?fft（xN）;???????????????????????????????????????%進行FFT?變換
mag?=?abs（y）;???????????????????????????????????%求幅值
f?=?（0:length（y）?-1）‘*fs/?length（y）;????????%進行對應的頻率轉(zhuǎn)換
subplot（232）;
plot（fmag）;????????????????????????????????????????%作頻譜圖
axis（[0100080]）;
xlabel（‘頻率/?Hz‘）;
ylabel（‘幅值‘）;
title（‘幅頻譜圖‘）;
grid;
%求均方根譜
sq?=?abs（y）;
subplot（233）;
plot（fsq）;
xlabel（‘頻率/?Hz‘）;
ylabel（‘均方根譜‘）;
title（‘均方根譜‘）;
grid;
%求功率譜
power?=?sq.^2;
subplot（234）;
plot（fpower）;
xlabel（‘頻率/?Hz‘）;
ylabel（‘功率譜‘）;
title（‘功率譜‘）;
grid;
%求對數(shù)譜
ln?=?log（sq）;
subplot（235）;
plot（fln）;
xlabel（‘頻率/?Hz‘）;
ylabel（‘對數(shù)譜‘）;
title（‘對數(shù)譜‘）;
grid;
%用IFFT?恢復原始信號
xifft?=?ifft（y）;
magx?=?real（xifft）;
ti?=?[0:length（xifft）-1]?/?fs;
subplot（236）;
plot（timagx）;
xlabel（‘時間/?s‘）;
ylabel（‘幅值‘）;
title（‘IFFT?后的信號波形‘）;
grid;
%****************2.白噪聲****************%
fs?=50;??????????????????????????????????????????????%設(shè)定采樣頻率
t?=?-5:0.1:5;
x?=?rand（1100）;
figure（2）;
subplot（231）;
plot（t（1:100）x）;???????????????????????????????%作白噪聲的時域波形
xlabel（‘時間（s）‘）;
ylabel（‘幅值‘）;
title（‘時域波形‘）;
grid;
%進行FFT?變換并做頻譜圖
y?=?fft（x）;?????????????????????????????????????%進行FFT?變換
mag?=?abs（y）;???????????????????????????????%求幅值
f?=?（0:length（y）?-1）‘*fs/?length（y）;?%進行對應的頻率轉(zhuǎn)換
subplot（232）;
plot（fmag）;????????????????????????????????????%作頻譜圖
xlabel（‘頻率/?Hz‘）;
ylabel（‘幅值‘）;
title（‘幅頻譜圖‘）;
grid;
%求均方根譜
sq?=?abs（y）;
subplot（233）;
plot（fsq）;
xlabel（‘頻率/?Hz‘）;
ylabel（‘均方根譜‘）;
title（‘均方根譜‘）;
grid;
%求功率譜
power?=?sq.^2;
subplot（234）;
plot（fpower）;
xlabel（‘頻率/?Hz‘）;
ylabel（‘功率譜‘）;
title（‘功率譜‘）;
grid;
%求對數(shù)譜
ln?=?log（sq）;
subplot（235）;
plot（fln）;
xlabel（‘頻率/?Hz‘）;
ylabel（‘對數(shù)譜‘）;
title（‘對數(shù)譜‘）;
grid;
%用IFFT?恢復原始信號
xifft?=?ifft（y）;
magx?=?real（xifft）;
ti?=?[0:length（xifft）-1]?/?fs;
subplot（236）;
plot（timagx）;
xlabel（‘時間/?s‘）;
ylabel（‘幅值‘）;
title（‘IFFT?后的信號波形‘）;
grid;