24小时热门版块排行榜    

查看: 1655  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 271材料工程求调剂 +6 .6lL 2026-03-18 6/300 2026-03-19 15:41 by haoshis
[考研] 本人考085602 化学工程 专硕 +17 不知道叫什么! 2026-03-15 19/950 2026-03-19 15:06 by 尽舜尧1
[考研] 281求调剂(0805) +9 烟汐忆海 2026-03-16 19/950 2026-03-19 11:42 by laoshidan
[考研] 0703化学调剂 +8 妮妮ninicgb 2026-03-15 12/600 2026-03-19 05:55 by anny19840123
[考研] 0703化学 305求调剂 +4 FY_yy 2026-03-14 4/200 2026-03-19 05:54 by anny19840123
[考研] 085601专硕,总分342求调剂,地区不限 +5 share_joy 2026-03-16 5/250 2026-03-18 14:48 by haxia
[考研] 一志愿西南交大,求调剂 +4 材化逐梦人 2026-03-18 4/200 2026-03-18 14:22 by 007_lilei
[考研] 收复试调剂生 +4 雨后秋荷 2026-03-18 4/200 2026-03-18 14:16 by elevennnne
[考研] 0703化学调剂 +3 妮妮ninicgb 2026-03-17 3/150 2026-03-18 10:29 by macy2011
[基金申请] 被我言中:新模板不强调格式了,假专家开始管格式了 +4 beefly 2026-03-14 4/200 2026-03-17 22:04 by 黄鸟于飞Chao
[考研] 268求调剂 +7 好运连绵不绝 2026-03-12 8/400 2026-03-17 20:28 by xilongliang
[考研] 278求调剂 +3 Yy7400 2026-03-13 3/150 2026-03-17 08:24 by laoshidan
[基金申请] 今年的国基金是打分制吗? 50+3 zhanghaozhu 2026-03-14 3/150 2026-03-16 17:07 by 北京莱茵润色
[考研] 085600材料与化工 求调剂 +13 enenenhui 2026-03-13 14/700 2026-03-16 15:19 by 了了了了。。
[考研] 327求调剂 +6 拾光任染 2026-03-15 11/550 2026-03-15 22:47 by 拾光任染
[考研] 求老师收留调剂 +4 jiang姜66 2026-03-14 5/250 2026-03-15 20:11 by Winj1e
[考研] 材料与化工 323 英一+数二+物化,一志愿:哈工大 本人本科双一流 +4 自由的_飞翔 2026-03-13 5/250 2026-03-14 19:39 by hmn_wj
[考研] 一志愿哈工大材料324分求调剂 +5 闫旭东 2026-03-14 5/250 2026-03-14 14:53 by 木瓜膏
[考研] 330求调剂 +3 ?酱给调剂跪了 2026-03-13 3/150 2026-03-14 10:13 by JourneyLucky
[考研] 329求调剂 +3 miaodesi 2026-03-12 4/200 2026-03-13 20:53 by 18595523086
信息提示
请填处理意见