matlab实现回归分析

一、实验目的

      回归分析是一种寻求多个变量之间关系的方法。通过实验要求掌握回归基本原理和处理办法。

二、实验原理

 

        (3)F检验

 三、实验内容

      用x光机检查镁合金铸件内部缺陷时,为了获得最佳的灵敏度,透视电压y应随透视件的厚度x而改变,经实验获得下表所示一组数据,假设透视件的厚度x无误差,试求透视电压y随厚度x变化的经验公式。

     

 四、实验总结

1、实验程序

clc
close all
clear

x = [12 13 14 15 16 18 20 22 24 26];
y = [52 55 58 61 65 70 75 80 85 91];

N = 10;
xy = 0;
for i=1:1:N
    xy = xy+x(i)*y(i);
end

b = (N*xy - sum(x)*sum(y))/(N*sum(x.^2)-sum(x)^2);
b0 = (sum(x.^2)*sum(y)-sum(x)*xy)/(N*sum(x.^2)-sum(x).^2);

ave_x = 1/N*sum(x);
ave_y = 1/N*sum(y);

lxx = sum((x-ave_x).^2);
lxy = sum((x-ave_x).*(y-ave_y));
lyy = sum((y-ave_y).^2);

%方差分析
S = lyy;
U = b*lxy;
Q = lyy - b*lxy;

%F检验
F = U/(Q/(N-2));

%绘制散点图
x1 = 10:1:30;
y1 = b0 + b*x1;

figure
plot(x,y,'ro');
hold on
plot(x1,y1,'k');
title(['y=',num2str(b0),'+','x*',num2str(b)]);

2、给出运行结果,画出散点图和回归曲线

  • 1
    点赞
  • 11
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

薄情书生

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值