24.实现线性拟合和相关系数（matlab程序）

1.简述

      

1. 基本语法
1.1 corr函数基本语法
语法    说明
rho = corr(X)    返回输入矩阵X中每对列之间的两两线性相关系数矩阵。
rho = corr(X, Y)    返回输入矩阵X和Y中每对列之间的两两相关系数矩阵。
[rho, pval] = corr(X, Y)    返回pval，一个p值矩阵，用于检验无相关性的假设与非零相关性的备择假设。
[rho, pval] = corr(___, Name, Value)    除了前面语法中的输入参数外，还指定使用一个或多个Name-Value对参数的选项，如：‘Type’, ‘Kendall’ …
1.2 corrcoef函数基本语法
语法    说明
R = corrcoef(A)    返回 A 的相关系数的矩阵，其中 A 的列表示随机变量，行表示观测值。
R = corrcoef(A, B)    返回两个随机变量 A 和 B 之间的系数。
[R, P] = corrcoef(___)    返回相关系数的矩阵和 p 值矩阵，用于测试观测到的现象之间没有关系的假设（原假设）。
此语法可与上述语法中的任何参数结合使用。
如果 P 的非对角线元素小于显著性水平（默认值为 0.05），则 R 中的相应相关性被视为显著。
如果 R 包含复数元素，则此语法无效。
[R, P, RL, RU] = corrcoef(___)    包括矩阵，这些矩阵包含每个系数的 95% 置信区间的下界和上界。如果 R 包含复数元素，则此语法无效。
___ = corrcoef(___, Name, Value)    在上述语法的基础上，通过一个或多个 Name, Value 对组参数指定其他选项以返回任意输出参数。
例如，corrcoef(A, ‘Alpha’, 0.1) 指定 90% 置信区间，corrcoef(A, ‘Rows’, ‘complete’) 省略 A 的包含一个或多个 NaN 值的所有行。
 

2.代码

function chengxu()
x1=[0.18
    0.30
    0.10
    0.22
    0.34
    0.14
    0.26]*1e-2;
x2=[28
    22
    32
    26
    20
    30
    24];
x3=[48
    45
    42
    39
    36
    33
    30];
y=[0.43
   0.41
   0.39
   0.49
   0.4
   0.42
   0.42];
A=[x1,x2,x1.^2,x2.^2,x3.^2];
b=y;
x=A\b;

yy=A*x;

plot(y)
hold on
plot(yy,'ro-')
legend({'实验值 ','拟合值'})
r=XGXS(y,yy)

end

function r=XGXS(X,Y)
%计算相关系数
X1=mean(X);
Y1=mean(Y);
fenzi=sum((X-X1).*(Y-Y1));
fenmu=sqrt(sum((X-X1).^2)*sum((Y-Y1).^2));
r=fenzi/fenmu;  %相关系数
end

3.运行结果