matlab非线性数据拟合程序编写

2025-03-29 00:58:12

推荐回答（1个）

回答1：

function fxnihe2()
format long
clc
clear
close all
x=[1970    1971    1972   1973  1974    1975    1976 1977     1978  ...
   1979 1980  1981    1982  1983   1984   1985    1986   1987    1988 ...
   1989       1990    1991      1992     1993     1994    1995    1996 ...
   1997   1998    1999     2000     2001   2002     2003     2004    2005 ...
   2006    2007    2008      2009      2010     2011     2012];
y=[82992, 85229 ,87177,89211, 90859 ,92420 , 93717, 94974  ,96259 , ...
    97542,98705 ,100072 , 101654 , 103008  ,104357 , 105851  ,107507 ,...
    109300,111026, 112704 , 114333,  115823, 117171 , 118517 ,119850 ,...
    121121, 122389, 123626 ,124761 , 125786  , 126743 ,127627 , 128453, ...
    129227 , 129988 , 130756 , 131448 , 132129  ,132802  ,133450  ,134091 ,...
    134735, 135404];
%x=x/mean(x);
%y=y/mean(y);
x0=[3.42*mean(y),100,3.73/mean(x)];
LB=[];
UB=[];
z=lsqcurvefit(@(x,t) curve(x,t),x0,x,y,LB,UB)
disp('参数如下')
zz=z;
yy=curve(zz,x);
plot(x,y,'k*')
hold on
plot(x,yy,'r') 
legend('原始数据散点图','拟合后函数曲线图')
%残差
disp('    原始数值             拟合数值             两者之间差值');
disp([y',yy',yy'-y'])
disp('残差平方和')
v=sum((yy-y).^2)


function y=curve(x,t)
a=x(1);b=x(2);r=x(3);
y=a./(1+b*exp(-r*t));

最新问答

水榭花园怎么样？好不好？值不值得买？

二氯化二硫（S 2 Cl 2 ）是广泛用于橡胶工业的硫化剂，其分子结构如图所示。常温下，S 2 Cl 2 是一种橙黄

她现在事多，心烦，我一说话他就烦，我该怎么办才能让他开心起来

轻度耳聋要双耳戴助听器吗？

怎样能计算出m20X1.5-6h 螺纹塞规的大径和中径

白度新闻

QQ浏览器,总是显示喔唷崩溃啦,重新加载,卸载后重装也不行,怎么解决

锐普德数控车床的补刀怎么用

电动推杆速度参数计算？

有个很像adidas的Logo的牌子叫什么名字？AD的是3条杠，它是5条杠的，是个美国鞋的牌子。