24小时热门版块排行榜

返回列表

杨小胖

金虫 (正式写手)

应助: 2 (幼儿园)
金币: 1007
散金: 100
红花: 6
帖子: 936
在线: 379小时
虫号: 454210
注册: 2007-11-08
性别: MM
专业: 测量与地图学

[求助] matlab曲线拟合 basic fitting 与polyfit

想请高手帮忙改个代码。
其实基本实现了，但怎么越改越乱。
问题
<1>我不知道怎么在图中把拟合的公式显示出来。
<2>当我只用一组数据时直接在命令窗口输入，那如果我有很多组这样的数据，那要怎么改呢？
x=[ 10.5 13 13.8 16 18 19 20.8];
y=[1190 1040 1000 800 600 430 200];
<3>为什么这样拟合出来的结果是曲线，但如果通过普通的polyfit拟合，出来是折线？

想要实现的是设置不一样的拟合次数，将拟合的系数都显示出来，并在每一个拟合系数求出来以后把表达式显示出来。

这是我从matlab宝典一书中找到的代码。

CODE:

function [th,err,yi]=polyfits(x,y,N,xi,r)

%x,y  :数据点系列

%N    :多项式拟合的系统

%r    :加权系数的逆矩阵

M=length(x);

x=x(:);

y=y(:);

%判断调用函数的格式，考虑权阵

if nargin==4

%当用户调用函数格式为(x,y,N,r)

    if length(xi)==M

        r=xi;

        xi=x;

%当用户调用函数格式为(x,y,N, xi)，权阵为1

    else r=1;

    end

%当用户调用函数格式为(x,y,N)，不考虑权阵

elseif nargin==3

    xi=x;

    r=1;

end

%求解系数矩阵

A(:,N+1)=ones(M,1);

for n=N:-1:1

    A(:,n)=A(:,n+1).*x;

end

if length(r)==M

    for m=1:M

        A(m,:)=A(m,:)/r(m);

        y(m)=y(m)/r(m);

    end

end

%计算拟合系数

th=(A\y)';

ye=polyval(th,x);

err=norm(y-ye)/norm(y);

yi=polyval(th,xi);

CODE:

clear;

clf;

x=input('请输入x值:');

y=input('请输入y值:');

% x=[ 10.5  13  13.8  16  18  19 20.8];

% y=[1190  1040 1000 800 600 430 200];

% eb=[0.2*ones(5,1);ones(5,1)];

[x,i]=sort(x);

y=y(i);

xi=min(x)+[0:100]/100*(max(x)-min(x));

 for i=1:4

 N=2*i-1;

[th,err1,yi]=polyfits(x,y,N,xi);

  subplot(220+i);

plot(x,y,'k*')

hold on 

plot(xi,yi,'g:','LineWidth',1.5)

title(['The ',num2str(N),'th Polynomial Curve Fitting'])

grid on

 end