MATLAB的答案

替代使用polyfit和polyval功能? ?

29 Ansichten(30天)的
本性欲
本性欲 3月12日。2021年
Kommentiert: 本性欲3月12日。2021年
% %这段代码看起来找到是否二次或线性情节给较低的错误对于一组给定的数据。
% %首先清除任何先前的数据脚本
clc;清晰的所有;关闭所有;
% %开始
%的输入值
d10 = (0.0024;0.0025;0.0056;0.0083;0.0262;0.0302;0.0413;0.0428);
(k = 2.29 e-04;2.28 e-04;2.40 e-04;1.21 e-04;2.58 e 03;4.92 e 03;5.16 e 03;7.00 e 03];
%情节线性和多项式
p1 = polyfit (d10, k, 1);
p2 = polyfit (d10, k, 2);
情节(d10, k,“o”);
持有在
情节(d10 polyval (p1, d10),“b -”);
情节(d10 polyval (p2, d10),“r——”);
包含(“d10 (mm)”,“字形大小”,20);
ylabel (的k(毫米/秒),“字形大小”,20);
xlim ([0, 0.05]);
ylim ([0, 0.01]);
%残留误差的线性
kfit1 = polyval (p1, d10);
kresid1 = k - kfit1;
SSresid1 =总和(kresid1。^ 2);
disp ([' R的平方值:num2str (SSresid1)]);
%残留误差二次
kfit2 = polyval (p2, d10);
kresid2 = k - kfit2;
SSresid2 =总和(kresid2。^ 2);
disp ([' R的平方值:num2str (SSresid2)]);
如果SSresid2 > SSresid1
disp ([的线性有更好的SSresid值:num2str (SSresid1)]);
其他的
disp ([“二次有更好的SSredid价值:”num2str (SSresid2)]);
结束
这是我写我的代码,发现一个线性或二次情节是否更好的为给定的数据。我不希望使用polyfit和polyval内置函数,谁能帮助找到一个替代?
1 Kommentar
1月
1月 3月12日。2021年
请代码格式。这次我为你这样做。

Melden您西奇,嗯祖茂堂kommentieren。

Melden您西奇一个,嗯这Frage祖茂堂beantworten。

Antworten (1)

1月
1月 3月12日。2021年
Bearbeitet:1月 3月12日。2021年
POLYFIT和POLYVAL是简单的功能。如果你的老师让你写他们自己的,我认为,这是一个主题在你的课。
这是一些配件的代码:
%使用:x, y, n =多项式次数
V =(元素个数(x) n + 1);%构建范德蒙矩阵:
为j = n: 1:1
V (:, j) = (:, j + 1)。* x (:);
结束
(Q, R) = qr (V, 0);%解决最小二乘问题:
p = ' (R \(转置(Q) * y (:)));%一样:(V \ y) '
对于POLYVAL使用霍纳方案。

Kategorien

Mehr祖茂堂描述性统计法登您在帮助中心和文件交换

标签

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!

翻译的