e3to14.m

来自「最小二乘法仿真y=20*sin(100*pi*t+pi/3)+4*sin(200」· M 代码 · 共 24 行

close all
clear all
clc
T=0.02;N=20;Ts=T/N;
t=0:0.02:12;
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);


f=inline('a(1)*sin(a(2)*t+pi/3)+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','t');

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

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

a',res;
subplot(111);
plot(t,y,'r');
xlabel('t/s');ylabel('幅值');
title('最小二乘法计算结果');