clc;clear all;
x=0:0.01:4*pi;
y=(sin(x)).^2;
plot(x,y)
hold on;
grid on;
x=[1978 ... 1989];
y=[0.1825 ... 0.1504];
x1=[ones(size(x)) x x.^2];
b=regress(y,x1)
x,y自己赋值,得到的b就是对应的a,b,c
ezplot('y=(sin(x))^2')

二次回归方程应该用,regress是线性回归
p=polyfit(x,y,n)
将方程式编写为一个函数,比如fun函数;
将x值放在某一个数组中,比如X数组:x=[x1,x2,x3,...,xn];
将y值存放在另一个数组中,比如Y数组:y=[y1,y2,y3,...,yn];
则有以下编程:
x=[x1,x2,x3,...,xn];
n=length(x);
y=zeros(1,n);
for k=1:n
y(k)=fun(x(k));
end
把x,y中数据多的去掉一些
例如y数据多
for i=1:length(x)
y1(i)=y(i);
end
plot(x,y1)
貌似有点糊