%% % FFT实践及频谱分析 % %% %1.正弦波*% fs=100;%设定采样频率 N=128; n=0:N-1; t=n/fs; f0=10;%设定正弦信号频率 x=sin(2*pi*f0*t); %生成正弦信号 figure(1); subplot(231); plot(t,x);%作正弦信号的时域波形 xlabel(’t’); ylabel(‘y’); title(‘正弦信号y=2*pi*10t时域波形’); grid; %进行FFT变换并做频谱图 y=fft(x,N);%进行fft变换 mag=abs(y);%求幅值 f=(0:length(y)-1)'*fs/length(y);%进行对应的频率转换 figure(1); subplot(232); plot(f,mag);%做频谱图 axis([0,100,0,80]); xlabel(‘频率(Hz)'); ylabel(‘幅值’); title(‘正弦信号y=2*pi*10t幅频谱图N=128’); grid; %求均方根谱 sq=abs(y); figure(1); subplot(233); plot(f,sq); xlabel(‘频率(Hz)'); ylabel(‘均方根谱’); title(‘正弦信号y=2*pi*10t均方根谱’); grid; %求功率谱 power=sq.^2; figure(1); subplot(234); plot(f,power); xlabel(‘频率(Hz)'); ylabel(‘功率谱’); title(‘正弦信号y=2*pi*10t功率谱’); grid; %求对数谱 ln=log(sq); figure(1); subplot(235); plot(f,ln); xlabel(‘频率(Hz)'); ylabel(‘对数谱’); title(‘正弦信号y=2*pi*10t对数谱’); grid; %用IFFT恢复原始信号 xifft=ifft(y); magx=real(xifft); ti=[0:length(xifft)-1]/fs; figure(1); subplot(236); plot(ti,magx); xlabel(’t’); ylabel(‘y’); title(‘通过IFFT转换的正弦信号波形’); grid;

Continue reading

矩阵基础回顾

 加法: 矩阵的加法中,双方必须是同型矩阵(A行列与B行列全部相等)才可进行运算。 矩阵加法只要对应位直接相加即可,Aij + Bij = Cij  乘法: 若A是M x S (M行S列)矩阵,乘以B(S x N)矩阵,则会得到C(M x N)矩阵。 乘法中,不要求两个矩阵同型。但第一个矩阵的S列数一定要等于第二个矩阵的S行数。乘法是用A的行1中A11 A12 A13…Ams 乘以B矩阵的列1中B11 B21 B31…Bsn 作为C的C11,然后用A11 A12 A13乘以B的第二列B12 B22 B32作为C的C12。即A的i行乘以B的j列得出C中对应的Cij (ij<=mn)

Continue reading

Author's picture

LuLU

Love coding and new technologies

Cloud Solution Consultant

Canada