t120_2.m

一个谐波处理的matlab软件。只要将里面的数据换成自己的数据即可。

fs=6400;%设定采样频率
N=128;
n=0:N-1;
t=n/fs;
f0=50;%设定正弦信号频率
%生成正弦信号
%x=sin(2*pi*f0*t)+sin(4*pi*f0*t);
%x=sin(2*pi*f0*t)+sin(4*pi*f0*t)+sin(8*pi*f0*t);
%x=sin(2*pi*f0*t)+0.5 * sin(4*pi*f0*t)+sin(8*pi*f0*t);

%x1=[0.0792,1.2475,1.8988,1.8169,1.9050,2.0459,2.0853,2.1112,2.1102,2.0708,2.0023,1.8935,1.7392,1.5548,1.3475,1.1112,0.8360,0.5600,0.2662,0.0933,0.1062,0.0963,0.0947,0.0812,0.0846,0.0771,0.0698,0.0634,0.0563,0.0326,0.0087,-0.0029,0.0160,0.0237,0.0191,0.0123,0.0030,-0.0035,-0.2522,-0.2564,-0.2667,-0.2719,-0.2822,-1.6732,-2.0008,-2.0324,-2.0983,-2.2106,-2.2723,-2.3392,-2.3863,-2.3889,-2.3599,-2.2893,-2.1541,-1.9651,-1.7504,-1.5039,-1.2553,-0.9879,-0.7226,-0.4222,-0.3823,-0.3838,-0.1528,-1.3796,-1.8563,-1.8039,-1.9267,-2.0346,-2.0863,-2.1215,-2.1236,-2.0904,-2.0262,-1.9101,-1.7604,-1.5744,-1.3547,-1.1307,-0.8615,-0.5486,-0.2025,-0.0884,-0.0968,-0.0895,-0.0864,-0.0771,-0.0800,-0.0719,-0.0657,-0.0594,-0.0512,-0.0268,-0.0025,0.0066,-0.0160,-0.0190,-0.0170,-0.0124,-0.0051,0.0025,0.3115,0.2574,0.2636,0.2719,0.2791,1.6750,1.9630,2.0283,2.0904,2.1971,2.2667,2.3308,2.3791,2.3878,2.3599,2.2949,2.1661,1.9824,1.7611,1.5066,1.2548,0.9915,0.7009,0.3797,0.3843,0.3838];
%x2=[-0.0811,-0.0929,-0.0688,-0.0688,-0.0729,-0.0606,-0.0565,-0.0503,-0.0441,-0.0400,-0.0349,-0.0267,-0.0215,-0.0164,-0.0071,-0.0030,0.0020,0.0030,0.0020,0.0065,0.0091,0.0220,1.4057,1.8482,1.8436,1.9838,2.0825,2.1359,2.1790,2.2119,2.2242,2.2037,2.1010,1.9674,1.8051,1.6037,1.3612,1.1023,1.0571,0.8229,0.4597,0.3611,0.3796,1.7683,2.0968,2.1287,2.1903,2.3006,2.3557,2.4184,2.4585,2.4575,2.4200,2.3455,2.2027,2.0069,1.7845,1.5297,1.2708,1.0006,0.7263,0.4202,0.3749,0.3688,0.0421,0.0955,0.0667,0.0688,0.0678,0.0565,0.0513,0.0441,0.0380,0.0339,0.0267,0.0195,0.0133,0.0092,0.0010,-0.0028,-0.0066,-0.0071,-0.0048,-0.0064,-0.0128,-0.0236,-1.3294,-1.8562,-1.8138,-1.9335,-2.0527,-2.1030,-2.1405,-2.1780,-2.1862,-2.1533,-2.0342,-1.9006,-1.7311,-1.5236,-1.2950,-1.0489,-1.0756,-0.7900,-0.4140,-0.3585,-0.3770,-1.7660,-2.0557,-2.1194,-2.1739,-2.2775,-2.3402,-2.3999,-2.4408,-2.4400,-2.4088,-2.3372,-2.1996,-2.0105,-1.7837,-1.5236,-1.2647,-0.9955,-0.6996,-0.3760,-0.3791,-0.3678];
%x3=[1.4752,3.2736,1.1728,1.1560,1.6640,1.4496,1.4144,1.4592,1.4032,1.3504,1.3000,1.2176,1.1584,1.1200,1.0592,0.9096,0.6968,0.6080,0.6280,0.8128,0.8250,0.4488,0.4013,0.0152,0.0688,0.0192,-0.0608,-0.1040,-0.1872,-0.3152,-0.4160,-0.4864,-0.5040,-0.4988,-0.4684,-0.3520,-0.2936,-0.4080,-0.7936,-0.3168,-0.9952,-1.3664,-1.2864,-1.9133,-1.2360,-1.2904,-1.3964,-1.4527,-1.5120,-1.6072,-1.6800,-1.7280,-1.7600,-1.7456,-1.6304,-1.4624,-1.2695,-1.1160,-1.1484,-1.3792,-1.4160,-1.4316,-1.7920,-1.4684,-1.2224,-2.7040,-1.0344,-1.0224,-1.2960,-1.1808,-1.2128,-1.2064,-1.1488,-1.0784,-0.9976,-0.8896,-0.8120,-0.7488,-0.7488,-0.6972,-0.6048,-0.5120,-0.5976,-0.9371,-0.7392,-0.4792,-0.4520,-0.0944,-0.1787,-0.2366,-0.1016,-0.0552,-0.0059,0.0992,0.1888,0.2304,0.1928,0.1874,0.1088,-0.0288,-0.0096,0.3056,0.8144,0.1828,0.8712,1.0752,0.7312,1.4228,1.1312,1.1928,1.2624,1.3440,1.3745,1.4613,1.5408,1.5940,1.6150,1.5908,1.4826,1.2928,1.0944,1.0096,1.0088,1.0884,1.3392,1.7216,2.2528,1.6944];
%x4=[-0.1024,-0.5344,-0.0064,-0.0704,-0.2236,-0.0432,-0.0096,0.0048,0.0800,0.1312,0.1584,0.1673,0.1312,0.0128,-0.0320,0.1552,0.6000,0.7988,0.8208,0.9696,1.1048,1.1328,1.8704,1.0752,1.1304,1.2632,1.2776,1.3154,1.3800,1.4592,1.5072,1.5125,1.4120,1.2904,1.1496,0.9792,0.9184,1.0144,1.3312,0.8704,1.3408,2.1656,1.8948,2.3077,1.2460,1.4744,1.5432,1.5649,1.5873,1.6340,1.6164,1.5632,1.4976,1.4384,1.3280,1.2640,1.1456,1.1008,0.9632,0.7568,0.5824,0.5280,0.7424,0.4304,0.1024,0.7512,-0.0120,0.0208,0.0704,-0.0544,-0.0480,-0.0896,-0.1728,-0.2390,-0.2872,-0.3008,-0.2560,-0.1624,-0.1376,-0.2104,-0.4672,-0.8896,-1.2670,-1.5576,-1.6572,-1.5552,-2.3128,-1.0996,-1.1496,-1.3888,-1.3728,-1.3940,-1.4664,-1.5440,-1.5984,-1.6000,-1.4832,-1.3932,-1.2544,-1.0720,-1.0048,-1.1448,-1.4520,-0.8480,-1.3408,-2.0288,-1.4560,-1.9496,-1.2269,-1.3272,-1.3065,-1.3568,-1.3343,-1.3632,-1.3312,-1.2680,-1.1901,-1.0816,-0.9524,-0.8424,-0.7680,-0.7784,-0.7352,-0.6400,-0.5120,-0.4696,-0.7136,-0.3072];
%x1 = sin(2*pi*f0*t);
%x2 = sin(4*pi*f0*t);
%x3 = sin(2*pi*f0*t);
%x4 = sin(4*pi*f0*t);
x1=[2450,16221,18521,18205,19144,20231,20904,21164,21184,20796,20023,18801,17258,15573,13724,11610,9184,6398,3154,838,979,937,927,792,792,734,667,626,563,481,315,211,232,294,222,128,-2397,-2475,-2585,-2646,-2687,-2730,-2833,-16117,-19495,-20407,-21278,-22325,-23143,-23599,-23909,-23846,-23484,-22850,-21620,-20066,-18169,-16024,-13664,-10987,-8262,-5651,-3755,-3143,-803,-12304,-19144,-18334,-19432,-20530,-21194,-21454,-21444,-21089,-20376,-19154,-17672,-16014,-14127,-11776,-9019,-6283,-2885,-823,-947,-927,-879,-761,-761,-698,-647,-606,-532,-242,-35,87,46,-118,-160,3288,2407,2491,2532,2585,2698,2781,2865,15838,19291,20399,21273,22346,23298,23733,24036,23972,23640,23060,21951,20459,18573,16439,14025,11382,8595,5071,3776,3932];
x2=[-113,-729,-708,-678,-688,-565,-554,-503,-441,-400,-318,-297,-236,-184,-133,-71,-61,-10,10,35,91,156,12667,18133,18051,19314,20188,20619,21071,21297,21297,21092,20414,19160,17516,15770,16119,13982,11475,8640,4864,3508,3708,16879,20342,21228,22078,23105,23866,24277,24528,24431,24061,23334,22037,20393,18379,16201,13777,11085,8311,5599,3626,3010,707,863,714,667,719,585,544,472,426,421,349,308,226,174,92,51,51,30,10,-10,-71,-174,-11794,-17357,-17383,-18431,-19525,-19998,-20373,-20727,-20701,-20383,-19574,-18061,-16284,-17897,-15030,-12667,-10304,-7469,-4017,-3647,-3852,-16694,-20198,-21256,-22078,-23116,-24030,-24431,-24621,-24503,-24143,-23475,-22325,-20783,-18847,-16592,-14105,-11403,-8537,-4993,-3667,-3785];
x3=[16200,11140,14720,12008,16464,14848,14544,14024,13344,12616,11936,10914,10378,10372,10496,9984,9152,8032,9352,11280,10768,7536,8656,1776,145,200,-1016,-1856,-3456,-4992,-6064,-7096,-8288,-9072,-9088,-8480,-7936,-5068,-5440,-5472,-9040,-9568,-9472,-16704,-11152,-12906,-13440,-14316,-15392,-16248,-16928,-17168,-17352,-17536,-17008,-15840,-13984,-11952,-10488,-13008,-15856,-12120,-12656,-17144,-16593,-32768,-12412,-8896,-13232,-11560,-11520,-10928,-10008,-9100,-8192,-7104,-6392,-6016,-6224,-6456,-6720,-7640,-8856,-15360,-15072,-11136,-12208,-5376,-3280,-4088,-2224,-1344,-880,87,688,960,1162,704,-112,1200,96,3248,3088,2304,4528,8416,5248,15616,10512,11980,12400,13328,13968,14720,15400,15664,15744,15752,14937,13536,11552,9580,7938,7760,8532,12880,20448,16185];
x4=[-5280,92,-1824,-862,-2416,-768,-256,624,1616,2192,2560,2784,2520,1504,40,-480,768,3584,3152,5824,7424,5072,16736,11648,10144,11888,11696,12288,13152,13520,13696,14144,13920,12672,11144,9728,9216,7832,9216,9472,12904,20196,20048,25856,14372,14882,15136,15124,15680,15400,14944,14176,13352,12704,11992,11408,10548,9648,9792,8600,6512,6064,7458,9984,7935,13257,1000,-408,1224,-640,-832,-1824,-3040,-3856,-4352,-4592,-4248,-3456,-1904,-2910,-5806,-5968,-5480,-8936,-10688,-7760,-18680,-13820,-12144,-14224,-14040,-13824,-14263,-14865,-15140,-15168,-14592,-13280,-12016,-13112,-11192,-13648,-13776,-13136,-14624,-19024,-12856,-23204,-13296,-13268,-13312,-13180,-13276,-12504,-12000,-11340,-10480,-9423,-8047,-6512,-5055,-4256,-4864,-6048,-6144,-5424,-8928,-10039];

%150


figure(1);
subplot(2,6,1);
plot(t,x1);%作正弦信号的时域波形
xlabel('t');
ylabel('y');
%title('正弦信号y=2*pi*10t时域波形');
grid;

%进行FFT变换并做频谱图
y1=fft(x1,N);%进行fft变换
mag=abs(y1);%求幅值
m=length(y1);
f=(0:m/2-1)'*fs/m;%进行对应的频率转换
figure(1);
subplot(2,6,2);
bar(f,mag(1:m/2));%做频谱图
%axis([0,6500,0,2]);
axis([0,1000,0,1280000]);
xlabel('频率(Hz)');
ylabel('幅值');
%title('正弦信号y=2*pi*10t幅频谱图N=128');
grid;

z1 = sqrt((mag(3)*mag(3) + mag(4)*mag(4) + mag(6)*mag(6) + mag(8)*mag(8) + mag(10)*mag(10) + mag(12)*mag(12) + mag(14)*mag(14)) / (mag(2) * mag(2)));
mag1 = mag;

y1(2) = 0;
y1(128) = 0;
xifft=ifft(y1);
magx=real(xifft);
ti=[0:length(xifft)-1]/fs;
figure(1);
subplot(2,6,3);
plot(ti,magx);
xlabel('t');
ylabel('y');
%title('通过IFFT转换的正弦信号波形');
grid;

figure(1);
subplot(2,6,4);
plot(t,x2);%作正弦信号的时域波形
xlabel('t');
ylabel('y');
%title('正弦信号y=2*pi*10t时域波形');
grid;

%进行FFT变换并做频谱图
y2=fft(x2,N);%进行fft变换
mag=abs(y2);%求幅值
m=length(y2);
f=(0:m/2-1)'*fs/m;%进行对应的频率转换
figure(1);
subplot(2,6,5);
bar(f,mag(1:m/2));%做频谱图
%axis([0,6500,0,2]);
axis([0,1000,0,1280000]);
xlabel('频率(Hz)');
ylabel('幅值');
%title('正弦信号y=2*pi*10t幅频谱图N=128');
grid;

z2 = sqrt((mag(3)*mag(3) + mag(4)*mag(4) + mag(6)*mag(6) + mag(8)*mag(8) + mag(10)*mag(10) + mag(12)*mag(12) + mag(14)*mag(14)) / (mag(2) * mag(2)));
mag2 = mag;

y2(2) = 0;
y2(128) = 0;
xifft=ifft(y2);
magx=real(xifft);
ti=[0:length(xifft)-1]/fs;
figure(1);
subplot(2,6,6);
plot(ti,magx);
xlabel('t');
ylabel('y');
%title('通过IFFT转换的正弦信号波形');
grid;

figure(1);
subplot(2,6,7);
plot(t,x3);%作正弦信号的时域波形
xlabel('t');
ylabel('y');
%title('正弦信号y=2*pi*10t时域波形');
grid;

%进行FFT变换并做频谱图
y3=fft(x3,N);%进行fft变换
mag=abs(y3);%求幅值
m=length(y3);
f=(0:m/2-1)'*fs/m;%进行对应的频率转换
figure(1);
subplot(2,6,8);
bar(f,mag(1:m/2));%做频谱图
%axis([0,6500,0,2]);
axis([0,1000,0,1280000]);
xlabel('频率(Hz)');
ylabel('幅值');
%title('正弦信号y=2*pi*10t幅频谱图N=128');
grid;

z3 = sqrt((mag(3)*mag(3) + mag(4)*mag(4) + mag(6)*mag(6) + mag(8)*mag(8) + mag(10)*mag(10) + mag(12)*mag(12) + mag(14)*mag(14)) / (mag(2) * mag(2)));
mag3 = mag;

y3(2) = 0;
y3(128) = 0;

xifft=ifft(y3);
magx=real(xifft);
ti=[0:length(xifft)-1]/fs;
figure(1);
subplot(2,6,9);
plot(ti,magx);
xlabel('t');
ylabel('y');
%title('通过IFFT转换的正弦信号波形');
grid;

figure(1);
subplot(2,6,10);
plot(t,x4);%作正弦信号的时域波形
xlabel('t');
ylabel('y');
%title('正弦信号y=2*pi*10t时域波形');
grid;

%进行FFT变换并做频谱图
y4=fft(x4,N);%进行fft变换
mag=abs(y4);%求幅值
m=length(y4);
f=(0:m/2-1)'*fs/m;%进行对应的频率转换
figure(1);
subplot(2,6,11);
bar(f,mag(1:m/2));%做频谱图
%axis([0,6500,0,2]);
axis([0,1000,0,1280000]);
xlabel('频率(Hz)');
ylabel('幅值');
%title('正弦信号y=2*pi*10t幅频谱图N=128');
grid;

z4 = sqrt((mag(3)*mag(3) + mag(4)*mag(4) + mag(6)*mag(6) + mag(8)*mag(8) + mag(10)*mag(10) + mag(12)*mag(12) + mag(14)*mag(14)) / (mag(2) * mag(2)));
mag4 = mag;

y4(2) = 0;
y4(128) = 0;
xifft=ifft(y4);
magx=real(xifft);
ti=[0:length(xifft)-1]/fs;
figure(1);
subplot(2,6,12);
plot(ti,magx);
xlabel('t');
ylabel('y');
%title('通过IFFT转换的正弦信号波形');
grid;