一元线性回归模型的拟合优度检验的matlab代码

一元线性回归模型的拟合优度检验的matlab代码,第1张

主要是用regress函数来进行:给你举个例子来说明吧。

x=[0 1 2 3 4 ]';y=[10 13 15,20 23]';

x=[ones(5,1),x]; %给出两个数组元素

[b,bint,r,rint,stats]=regress(y,x,005); %对x和y进行一元线性回归,并得到相关系数,其中,stats中第一个数即为相关系数,大于09就认为拟合很好。

结果:stats =

09829 1719474 00010 00063

即为09829

就sym('3')与sym(3)而言,二者没有区别:

>> sym('3')==sym(3)

ans =

     1

 但如果不是整数,结果可能不一样:

>> sym('sqrt(2)-1')

ans =

sqrt(2)-1

>> sym(sqrt(2)-1) 

ans = 

74618081806211082^(-54)

前者以字符串的形式传给符号运算内核,可以保留完整的精度;而后者经过浮点数运算之后再转换为符号类型,存在精度损失。

可这样:

a=[1 2 4];

fai=[0 pi/3 2pi/3];

t=0:pi/180:2pi;

for k=3:-1:1

x=a(k)(2cos(t)-cos(2t))-a(k);

y=a(k)(2sin(t)-sin(2t));

[s,r]=cart2pol(x,y);

polar(s+fai(k),r)

hold on

axis equal

end

可以对ezplot返回线的句柄进行操作

或者用plot代替ezplot绘制,在plot中可以直接指定线形

h = ezplot()

这里h就是返回的句

柄操作用可以用set函数

h = ezplot(@(x)x^2);

set(h,'color','r','linewidth',2)

这就是一个高次的三元方程, 在空间上对应平面图形

要画出来需要点技巧, 先解出函数表达式来, 再画

命令如下:

Y=solve('(x^2 + (9/4)y^2 + z^2 - 1)^3 - x^2z^3 - (9/80)y^2z^3 = 0','y')

ezmesh(Y(1),200);

hold on

ezmesh(Y(2),200)

axis auto

axis equal

还有几个解应该是复空间里面的, 画不出来图形

这是笛卡尔的著名的“心形线”他花出来应该是颗心的形状

(下面就是那个故事)

1956年,斯德哥尔摩的街头,52岁的笛卡尔邂逅了18岁的瑞典公主克里斯汀。

那时,落魄、一文不名的笛卡尔过着乞讨的生活,全部的财产只有身上穿得破破烂烂的衣服和随身所带的几本数学书籍。

  一个宁静的午后,笛卡尔照例坐在街头。突然,有人来到他旁边,拍了拍他的肩膀:“你在干什么呢?”扭过头,笛卡尔看到一张年轻秀丽的脸庞,一双清澈的眼睛如湛蓝的湖水,楚楚动人。她就是瑞典的小公主,国王最宠爱的女儿克里斯汀。

  她蹲下身,拿过笛卡尔的数学书和草稿纸,和他交谈起来。言谈中,他发现,这个小女孩思维敏捷,对数学有着浓厚的兴趣。

  几天后,他意外地接到通知,国王聘请他做小公主的数学老师。满心疑惑的笛卡尔跟随前来通知的侍卫一起来到皇宫,在会客厅等候的时候,他看到前几天在街头偶遇的女孩子。从此,他当上了公主的数学老师。

  公主的数学在笛卡尔的悉心指导下突飞猛进,他们之间也开始变得亲密起来。每天形影不离也使他们彼此产生了爱慕之心。

  然而,没过多久,他们的恋情传到了国王的耳朵里。国王大怒,下令马上将笛卡尔处死。在克里斯汀的苦苦哀求下,国王将他放逐回国,公主被软禁。

  笛卡尔回到法国后不久,便染上重病。在生命进入倒计时的那段日子,他日夜思念的还是街头偶遇的那张温暖的笑脸。在笛卡尔给克里斯汀寄出第十三封信后,他永远地离开了这个世界。这最后一封信上没有写一句话,只有一个方程:r=a(1-sinX)。

  国王不忍看着心爱的女儿每天闷闷不乐,便把这封信给了她。拿到信的克里斯汀立即明白了恋人的意图,找来纸和笔,着手把方程图形画了出来,一颗心形图案出现在眼前,克里斯汀不禁流下感动的泪水,这条曲线就是著名的“心形线”。

  这封享誉世界的另类情书,至今还保存在欧洲笛卡尔的纪念馆里。

欢迎分享,转载请注明来源:表白网

原文地址:https://h5.hunlipic.com/biaobai/4056651.html

(0)
打赏 微信扫一扫微信扫一扫 支付宝扫一扫支付宝扫一扫
上一篇 2024-04-13
下一篇2024-04-13

发表评论

登录后才能评论

评论列表(0条)

    保存