%%%程序编写者 西北工业大学自动化学院 Email: yincwxa2013@mail.nwpu.edu.cn
%% All rights reserved
clear
clc
x=input('输入节点坐标x=')
y=input('输入节点坐标函数值y=')
w=input('输入权重向量w=')
n=input('输入拟合最高次数n=')
for i=1:n+1
for j=1:n+1
p(i,j)=(w.*(x.^(i-1)))*(x.^(j-1))';
end
end
for i=1:n+1
q(i)=(w.*y)*(x.^(i-1))';
end
C=inv(p)*q'
f=poly2sym(C)
%%%
输入节点坐标x=[1 1 2 3 4 ]
x =
1 1 2 3 4
输入节点坐标函数值y=[2 3 5 6 9]
y =
2 3 5 6 9
输入权重向量w=2
w =
2
输入拟合最高次数n=2
n =
2
C =
0.9231
1.5641
0.1026
f =
32477881928153/35184372088832*x^2+1719748956265/1099511627776*x+7217307095145/70368744177664
>>
转载本文请联系原作者获取授权,同时请注明本文来自殷春武科学网博客。
链接地址:http://blog.sciencenet.cn/blog-748007-751138.html
上一篇:数值分析--雅克比迭代法、高斯赛德尔迭代法、超松弛迭代法matlab
下一篇:数值分析--非线性方程组不动点迭代法matlab程序