机械测试MATLAB的应用

1、首先，采样点数不变，改变采样频率信号频率取50Hz采样频率取50、100、200、1000四种情况。采样点数N=128.程序：fs=50（100、200、1000）; %设定采样频率N=128;n=0:N-1;t=n/fs;f0=50; %设定正弦信号频率%生成正弦信号x=sin(2*pi*f0*t);figure(1);subplot(111);plot(t,x);%作正弦信号的时域波形xlabel('t');ylabel('y');title('正弦信号y=2*pi*50t时域波形');grid;f0=50Hz时，产生混叠 f0=100Hz时

2、，理论上恰好不产生混叠，但实际上已经错误。f0=200Hz时，比较尖，不平滑f0=1000Hz时，比较平滑下面做FFT变换两种：一种是MATLAB自带的FFT变换程序，另一种是自编的FFT变换程序。先介绍自编程序：y=myfft（x，n）function y= myfft(xr,n)p=0:n-1;% 开始倒位序nu=log2(n);p1=p;b=zeros(1,n);for t=1:nu; p2=floor(p1/2); b=b*2+(p1-2*p2); p1=p2;end;yr(p+1)=xr(b+1);xr=yr; % 倒位序结束t=0:n/2-1; %计算因子w开始 (只计算w0到w

3、n/2-1)for v=0:n/2-1;w=exp(-2*i*pi*t/n);end; %计算因子w结束for m=1:nu;% 计算x(k)开始 h=2(m-1); k=1; while(k<n+1) for t=1:h; y=bitshift(k-1,nu-m,nu)+1; %求w的幂次数 xch(k)=xr(k)+w(y)*xr(k+h); k=k+1; end; for t=1:h; y=bitshift(k-1-h,nu-m,nu)+1; %求w的幂次数 xch(k)=xr(k-h)-xr(k)*w(y); k=k+1; end; end; xr=xch;end;% 计算x(k

4、)结束y=xr%输出变换后的结果给出画图程序：fs=(自取);%设定采样频率N=(自取);n=0:N-1;t=n/fs;f0=50;%设定正弦信号频率%生成正弦信号x=sin(2*pi*f0*t);%进行FFT变换并做频谱图y=fft(x,N)（或者是myfft）;%进行fft变换mag=abs(y);%求幅值f=(0:length(y)-1)'*fs/length(y);%进行对应的频率转换figure(1);subplot(111);plot(f,mag);%做频谱图axis(0,100,0,80（可改）);xlabel('频率(Hz)');ylabel('幅值');title('正弦信号y=sin(2*pi*f0*t)幅频谱图');grid;给定信号频率50Hz，采样频率1000Hz，改变采样点数NN=128N=512N=1024给定信号频率50Hz，采样频率200Hz，采样点数N=128给定信号频