matlab中,关于多元函数的拟合matlab程序如何书写?(4自变量1因变量)

已知x1=[0.550.650.650.650.650.550.550.550.65]x2=[1.61.61.41.61.41.41.61.41.4]x3=[202020... 已知x1=[0.55 0.65 0.65 0.65 0.65 0.55 0.55 0.55 0.65]x2=[1.6 1.6 1.4 1.6 1.4 1.4 1.6 1.4 1.4] x3=[20 20 20 10 10 20 20 20 20]x4=[10 10 10 10 10 10 12 12 12]Y=[1.848 3.1453.3373.0223.1881.971.631.6212.534]求问怎么用matlab拟合出f(x1,x2,x3,x4)来?1.大概知道是用nlinfit或者lsqcurvefit来做,但是并不太了解区别?2.用多项式拟合的话需要先选用合适的模型,但是看了一下x与y的关系没有发现适合的函数,所以可以随便选择吗?3.我看有些说需要给初值,初值是怎么得来的呢?4.如何看拟合后的曲线相关性?是用这个[b,bint,r,rint,stats]=regress(Y,X)吗?
我如图这样书写之后,错误提示Not enough input arguments.
请问是因为要求的系数太多,但是数据给的不够吗?
展开
 我来答
lhmhz
高粉答主

2019-01-24 · 专注matlab等在各领域中的应用。
lhmhz
采纳数:7263 获赞数:17001

向TA提问 私信TA
展开全部

题主的(4自变量1因变量)多元函数的拟合matlab程序。可以用nlinfit非线性回归函数来做(也可以用lsqcurvefit函数)。实现代码:

x1=[0.55 0.65 0.65 0.65 0.65 0.55 0.55 0.55 0.65]';

x2=[1.6 1.6 1.4 1.6 1.4 1.4 1.6 1.4 1.4]';

x3=[20 20 20 10 10 20 20 20 20]';

x4=[10 10 10 10 10 10 12 12 12]';

Y=[1.848 3.145 3.337 3.022 3.188 1.97 1.63 1.621 2.534]';

X=[x1 x2 x3 x4];

n=length(x1);

a0=rand(1,5);

func=@(a,X)(a(1)+a(2)*X(:,1)+a(3)*X(:,2)+a(4)*X(:,3)+a(5)*X(:,4));

[a,r,J] = nlinfit(X,Y,func,a0);a

Y1=func(a,X);

[Y  Y1]

运行上述代码,可以得到

a1= -1.52778571412534;a2= 11.4823809522886;a3=  -0.359047619087202 ;a4= 0.00441190476147387;a5=-0.23363095238449

多元函数表达式,y=a1+a2*x1+a3*x2+a4*x3+a5*x4

nlinfit与lsqcurvefit两者区别并不太大,前者用回归的方法来求解,而后者用最小二乘法来求解,两者都可以用于非线性函数和线性函数。

追问
还有两个问题:1.我报错的原因是因为试验数据不够吗?
2.决定系数是哪一句可以查看显示的呢?谢谢
追答
1、报错的原因,不是因为试验数据不够,而是数据应采用列向量的形式来表示,拟合表达式中的x1应用X(:,1)来表示,其他变量类同。

2、如何看拟合后的曲线相关性?可以根据线性函数的相关性公式(可以查看概率统计的书籍)来求。即用原始数据与拟合数据来计算。
Sievers分析仪
2024-10-13 广告
是的。传统上,对于符合要求的内毒素检测,最终用户必须从标准内毒素库存瓶中构建至少一式两份三点标准曲线;必须有重复的阴性控制;每个样品和PPC必须一式两份。有了Sievers Eclipse内毒素检测仪,这些步骤可以通过使用预嵌入的内毒素标准... 点击进入详情页
本回答由Sievers分析仪提供
暗送秋浡365
2019-01-24 · TA获得超过4660个赞
知道大有可为答主
回答量:6401
采纳率:78%
帮助的人:299万
展开全部
如果不确定应变量和自变量的关系,可以考虑使用神经网络来拟合
MATLAB有自带的神经网络工具箱,可以自己研究下,不需要编码,按照界面的要求自己一步步来就可以了。
已赞过 已踩过<
你对这个回答的评价是?
评论 收起
推荐律师服务: 若未解决您的问题,请您详细描述您的问题,通过百度律临进行免费专业咨询

为你推荐:

下载百度知道APP,抢鲜体验
使用百度知道APP,立即抢鲜体验。你的手机镜头里或许有别人想知道的答案。
扫描二维码下载
×

类别

我们会通过消息、邮箱等方式尽快将举报结果通知您。

说明

0/200

提交
取消

辅 助

模 式