MATLAB中利用cftool导出代码实现一张图中拟合多条平滑曲线

版权声明:本文为博主原创文章,未经博主允许不得转载。 https://blog.csdn.net/zyj1286076714/article/details/50375935

MATLAB中自带的cftool拟合工具箱不能将多条曲线同时画在同一副图中,而常规的plot()函数又不能拟合平滑直线,接下来总结一种可以利用cftool导出的代码,在一张图中拟合多条平滑曲线。

首先输入所要拟合的数据,如x, y, x1, y1, x2, y2等等。

之后打开cftool工具箱,使用数据拟合曲线,在拟合方式一栏选择Smoothing Spline。可以得到图像。

之后在文件菜单栏中点击Generate Code,之后将在工作空间里看到导出的m文件。

function [fitresult, gof] = createFit(x, y)
%CREATEFIT(X,Y)
%  Create a fit.
%
%  Data for 'untitled fit 1' fit:
%      X Input : x
%      Y Output: y
%  Output:
%      fitresult : a fit object representing the fit.
%      gof : structure with goodness-of fit info.
%
%  另请参阅 FIT, CFIT, SFIT.

%  由 MATLAB 于 22-Dec-2015 00:27:49 自动生成


%% Fit: 'untitled fit 1'.
[xData, yData] = prepareCurveData( x, y );

% Set up fittype and options.
ft = fittype( 'smoothingspline' );

% Fit model to data.
[fitresult, gof] = fit( xData, yData, ft );

% Plot fit with data.
figure( 'Name', 'untitled fit 1' );
h = plot( fitresult, xData, yData );
legend( h, 'y vs. x', 'untitled fit 1', 'Location', 'NorthEast' );
% Label axes
xlabel( 'x' );
ylabel( 'y' );
grid on

也就是说,其实也可以通过编写函数来实现拟合平滑曲线,但是由cftool工具箱导出的代码可以减少自己的工作量,减少由于MATLAB语言不熟练导致效率低下的工作。

之后在这个m文件的适当位置加入拟合多条曲线的代码,就可以实现与cftool拟合平滑曲线效果相同的多条曲线了。

以下是我修改的代码,注意注释。


function [fitresult, gof] = e1(x, y, x1, y1, x2, y2)     %将需要绘制图像的数据加入参数中
%CREATEFIT(X,Y)
%  Create a fit.
%
%  Data for 'untitled fit 1' fit:
%      X Input : x
%      Y Output: y
%  Output:
%      fitresult : a fit object representing the fit.
%      gof : structure with goodness-of fit info.
%
%  另请参阅 FIT, CFIT, SFIT.

%  由 MATLAB 于 21-Dec-2015 23:23:33 自动生成


%% Fit: 'untitled fit 1'.
[xData, yData] = prepareCurveData( x, y );
[xData1, yData1] = prepareCurveData( x1, y1 );<span style="white-space:pre">	</span>%第二条曲线需要的代码,注意改变参数
[xData2, yData2] = prepareCurveData( x2, y2 );<span style="white-space:pre">	</span>%第三条曲线需要的代码,注意改变参数

% Set up fittype and options.
ft = fittype( 'smoothingspline' );<span style="white-space:pre">		</span>%顾名思义,将拟合方式设置为光滑曲线

% Fit model to data.
[fitresult, gof] = fit( xData, yData, ft );
[fitresult1, gof1] = fit( xData1, yData1, ft );<span style="white-space:pre">	</span>%第二条曲线需要的代码,注意改变参数
[fitresult2, gof2] = fit( xData2, yData2, ft );<span style="white-space:pre">	</span>%第三条曲线需要的代码,注意改变参数

% Plot fit with data.
figure( 'Name', 'fit 1' );
h = plot( fitresult, xData, yData );
hold on;<span style="white-space:pre">					</span>%加入hold on 可以使之后的曲线画在同一张图中
h1 = plot( fitresult1, xData1, yData1, 'o'  );<span style="white-space:pre">	</span>%绘制第二条曲线,注意参数
hold on;<span style="white-space:pre">					</span>%作用同上
h2 = plot( fitresult2, xData2, yData2, 'o'  );<span style="white-space:pre">	</span>%绘制第三条曲线,注意参数
hold on;

% Label axes
xlabel( 'x' );
ylabel( 'y' );
grid on


保存之后,在命令行窗口中运行该m函数文件,即可得到在一张图中同时拟合三条光滑曲线。






同理,其他拟合方式也可以使用相同的方式达到多条曲线同时拟合的图像。



猜你喜欢

转载自blog.csdn.net/zyj1286076714/article/details/50375935