e3to15b.m

最小二乘法仿真y=20*sin(100*pi*t+pi/3)+4*sin(200*pi*t)+10*sin(300*pi*t)+2*sin(400*pi*t)+6*sin(500*pi*t)波形

close all
clear all
clc
T=0.02;N=20;Ts=T/N;
t=0:0.02:6;
y=20*sin(100*pi*t+pi/3)+4*sin(200*pi*t)+10*sin(300*pi*t)+2*sin(400*pi*t)+6*sin(500*pi*t)+20*exp(-t/0.03);


f=inline('a(1)*sin(a(2)*t+a(13))+a(3)*sin(a(4)*t)+a(5)*sin(a(6)*t)+a(7)*sin(a(8)*t)+a(9)*sin(a(10)*t)+a(11)*exp(a(12)*t)','a','t');

%建立函数原型，则可以根据他来进行下面的求取系数的计算

[a,res]=lsqcurvefit(f,[20,100*pi,4,200*pi,10,300*pi,2,400*pi,6,500*pi,20,-1/0.03,pi/3],t,y);

a',res
subplot(111);
plot(t,y,'-r');
xlabel('t/ms');ylabel('幅值/pu');axis([0,6,10,40]);
title('最小二乘法计算结果');