24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1684  |  回复: 6

bingmeng1987

木虫 (正式写手)

[求助] 拟合方程并作图 已有2人参与

我想对两个数组x,y进行线性拟合,输出拟合方程并作图,并且把拟合方程注明在图上。程序如下:
x=[1 2 3];
y=[2 4 6];
aa=polyfit(x,y,1)
plot(x,y,'g*',x,polyval(aa,x),'b-')
y=poly2sym(aa)
text(2,5,'y=poly2sym(aa)')
运行结果怎么是:y=2*x + 6879362064066739/5070602400912917605986812821504,不是应该是y=2*x吗?
图中注明的是y=poly2sym(aa),怎样才能让图中注明的是y=2*x?
Matlab不是很熟,请大家帮帮我,谢谢!
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Jakson_Tom

铜虫 (初入文坛)

【答案】应助回帖


感谢参与,应助指数 +1
bingmeng1987: 金币+1, 有帮助 2014-02-17 21:21:57
这是因为数据拟合过程中是有些误差的,这些误差是由于程序设计语言中数据类型中有效数字的取舍导致的,所以算出来的系数肯定不如大脑算出来的准确了,数据越多计算出来的系数越准确
2楼2014-02-17 15:34:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bucheron

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
bingmeng1987: 金币+9, ★★★★★最佳答案, 非常感谢 2014-02-17 21:22:09
1 设置一个精度,如果多项式某个系数小于该值,强制令其为零。
2 需要把symbol格式转换为string格式,使用char命令。

x=[1 2 3];
y=[2 4 6];
aa=polyfit(x,y,1);
for i=1:length(aa)
    if(aa(i))<eps('single')
        aa(i)=0;
    end
end
plot(x,y,'g*',x,polyval(aa,x),'b-')
y=poly2sym(aa)
text(2,5, strcat('y=',char(y)))
Born to suffer !Born to burn!
3楼2014-02-17 17:12:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bingmeng1987

木虫 (正式写手)

引用回帖:
3楼: Originally posted by bucheron at 2014-02-17 17:12:04
1 设置一个精度,如果多项式某个系数小于该值,强制令其为零。
2 需要把symbol格式转换为string格式,使用char命令。

x=;
y=;
aa=polyfit(x,y,1);
for i=1:length(aa)
    if(aa(i))<eps('single')
  ...

你好,如果程序为
aa=[127.3697 0]';
y=poly2sym(aa)
计算结果为什么是
y =(4481422917642905*x)/35184372088832
怎么样改程序能让结果显示
y=127.3697*x呢?
4楼2014-02-17 21:45:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bucheron

金虫 (小有名气)

aa=[127.3697 0]';
y=poly2sym(aa)
后边加一句
y=vpa(y)
就可以了
Born to suffer !Born to burn!
5楼2014-02-17 22:15:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bucheron

金虫 (小有名气)

vpa可以设定精度,譬如
y=vpa(y,7)

» 本帖已获得的红花(最新10朵)

Born to suffer !Born to burn!
6楼2014-02-17 22:16:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bingmeng1987

木虫 (正式写手)

送红花一朵
引用回帖:
6楼: Originally posted by bucheron at 2014-02-17 22:16:52
vpa可以设定精度,譬如
y=vpa(y,7)

谢谢!
7楼2014-02-18 09:00:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 bingmeng1987 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 070300化学学硕311分求调剂 +17 梁富贵险中求 2026-04-04 18/900 2026-04-08 22:53 by may_新宇
[考研] 求调剂 +7 吃口冰激凌 2026-04-07 7/350 2026-04-08 22:49 by bljnqdcc
[考研] 材料工程085601,270求调剂 +22 @ASDF1234 2026-04-08 23/1150 2026-04-08 22:03 by zhouyuwinner
[考研] 270求调剂 +3 031127 2026-04-06 4/200 2026-04-08 21:00 by 逆水乘风
[考研] 一志愿华东理工085601材料工程303分求调剂 +15 a1708 2026-04-06 15/750 2026-04-08 16:23 by luoyongfeng
[考研] 264求调剂 +11 麦小叮当 2026-04-07 11/550 2026-04-08 16:05 by 一只好果子?
[考研] 0702物理学学硕299求调剂 +4 祁柒连 2026-04-06 4/200 2026-04-08 13:56 by wutongshun
[考研] 305求调剂 +3 77Qi 2026-04-07 3/150 2026-04-07 23:18 by jp9609
[考研] 求调剂 +4 电气小神童 2026-04-04 6/300 2026-04-07 00:14 by guanxin1001
[考研] 求调剂 +5 小沢 2026-04-03 5/250 2026-04-06 22:45 by 875465
[考研] 材料调剂 +5 小刘同学吖吖 2026-04-06 5/250 2026-04-06 18:34 by sherry_1901
[考研] 一志愿C9的化学工程(085602) 340分,感觉校内调剂无望,求调剂 +12 万事宜臻 2026-04-04 12/600 2026-04-06 07:46 by 无际的草原
[考研] 一志愿211,化学学硕,310分,本科重点双非,求调剂 +13 努力奋斗112 2026-04-04 13/650 2026-04-06 07:13 by jj987
[考研] 求调剂 +10 Hll胡 2026-04-04 10/500 2026-04-05 20:09 by nepu_uu
[考研] 329求调剂 +17 miaodesi 2026-04-02 20/1000 2026-04-05 18:33 by 蓝云思雨
[考研] 081200-11408-276学硕求调剂 +4 崔wj 2026-04-04 5/250 2026-04-05 14:06 by imissbao
[考研] 359求调剂 +7 hhhhaaaa$ 2026-04-04 7/350 2026-04-04 18:49 by imissbao
[考研] 本9一志愿2 0854低分专硕286求调剂 +9 芒种111 2026-04-04 9/450 2026-04-04 11:01 by tangruihua
[考研] 262求调剂 +6 励志一定发文章 2026-04-02 7/350 2026-04-03 09:54 by linyelide
[考研] 326求调剂 +10 崽崽仔 2026-04-02 10/500 2026-04-03 09:08 by 帕尔马拉特
信息提示
请填处理意见