x = [1 2 3 4 5 6 7 8 9 ]; y = [9 7 6 3 -1 2 5 7 20]; p=polyfit(x,y,3); %数字代表拟合函数的阶数 xi=0:0.01:10; yi=polyval(p,xi); plot(x,y,\'r>\',xi,yi) for i=1:9 text(x(i)+0.1,y(i)+0.02,num2str(i)) ; %加上0.1使标号和点不重合，标记每个点 end
也可以使用GUI工具箱来产生多项式拟合
选定matlab 工作空间中x,y的值
然后点击matlab菜单栏的绘图，选择如下 或者直接在命令行窗口输入cftool
也可以指定函数进行拟合
拟合函数 f(t)=acos(kt)e^(wt)
syms t x=[0;0.4;1.2;2;2.8;3.6;4.4;5.2;6;7.2;8;9.2;10.4;11.6;12.4;13.6;14.4;15]; %数据必须为列向量 y=[1;0.85;0.29;-0.27;-0.53;-0.4;-0.12;0.17;0.28;0.15;-0.03;-0.15;-0.071;0.059;0.08;0.032;-0.015;-0.02];%数据必须为列向量 subplot(1,2,1) plot(x,y,\'r<\'); % 拟合函数 f(t)=acos(kt)e^(wt) f=fittype(\'a*cos(k*t)*exp(w*t)\',\'independent\',\'t\',\'coefficients\',{\'a\',\'k\',\'w\'}); polyfun=fit(x,y,f); %显示拟合函数 xi=0:0.01:20; yi=polyfun(xi); subplot(1,2,2) plot(x,y,\'rs\',xi,yi,\'.\');
拟合函数为 f(t)=0.9987cos(1.001t)e^-0.2066t
请发表评论