最小二乗法直線フィットのC言語とmatlab実現

最小二乗法直線フィットのC言語とmatlab実現
ここは直接コードに上がって、理論は導いて普通の最小の二乗法の導出を参照してC言語の実現を証明します

#include 
//            y=a+bx   a b
void LFM(double x[],double y[],int n)
{
	int i;
	double a=0,b=0;
	double sum1=0,sum2=0;
	double mean_x=0,mean_y=0;
	//     
	for(i=0;i

MATLAB実装

clc;clear;close all
%%     
x=[-0.5,-1.0,-1.5,-2.0,-2.5,-3,-3.5,-5,-7,-10,-0.1,-0.2,-0.3,-0.4];
y=[0.528453,1.04233,1.51821,2.03788,2.50965,3.01868,3.48255,4.99680,6.98300,9.98355,0.133300,0.232100,0.325205,0.422272];
%%      
%       y=a+bx   a b
para=polyfit(x,y,1);
b=para(1);
a=para(2);
%    
y_hat=a+b*x;
%%   
fig= figure(1);
ax = axes('Parent',fig);
p = plot(x,[y',y_hat']);
set(p(1),'DisplayName','    ','Marker','*','LineStyle','none','Color',[0 0 0]);
set(p(2),'DisplayName','    ','Color',[1 0 0]);
ylabel('    (  :KHz)');
xlabel('    (  :V)');
title('           ');
leg = legend(ax,'show');
set(leg,'Position',[0.75 0.8 0.1 0.1]);
axis tight;
axis([-10 0 0 10]);

好きならいいね、知らないあなたを助けたい: