24小时热门版块排行榜    

查看: 2046  |  回复: 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的回帖

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的回帖
查看全部 7 个回答

Jakson_Tom

铜虫 (初入文坛)

【答案】应助回帖

★
感谢参与,应助指数 +1
bingmeng1987: 金币+1, ★有帮助 2014-02-17 21:21:57
这是因为数据拟合过程中是有些误差的,这些误差是由于程序设计语言中数据类型中有效数字的取舍导致的,所以算出来的系数肯定不如大脑算出来的准确了,数据越多计算出来的系数越准确
2楼2014-02-17 15:34:31
已阅   回复此楼   关注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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 YXM1FDCSE8Bn 2026-11-01 6/300 2026-11-02 09:57 by dcc2yOIRGGTT
[找工作] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 YXM1FDCSE8Bn 2026-10-31 5/250 2026-11-02 09:42 by dcc2yOIRGGTT
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 YXM1FDCSE8Bn 2026-10-31 6/300 2026-11-02 09:36 by dcc2yOIRGGTT
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-02 09:31 by dcc2yOIRGGTT
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 xvJS1BYDxvBS 2026-10-31 11/550 2026-11-02 06:00 by xzcv0Tehd8lQ
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 1uaZt2kAZj2S 2026-10-30 16/800 2026-11-02 05:47 by xzcv0Tehd8lQ
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 14/700 2026-11-02 05:41 by xzcv0Tehd8lQ
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 21/1050 2026-11-02 05:41 by xzcv0Tehd8lQ
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +5 1uaZt2kAZj2S 2026-10-30 18/900 2026-11-02 05:34 by xzcv0Tehd8lQ
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +6 X5bHxOULxFqb 2026-10-30 27/1350 2026-11-02 05:28 by xzcv0Tehd8lQ
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 1uaZt2kAZj2S 2026-10-30 16/800 2026-11-02 02:47 by xzcv0Tehd8lQ
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 8/400 2026-11-01 21:35 by xzcv0Tehd8lQ
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 1uaZt2kAZj2S 2026-10-30 11/550 2026-11-01 21:35 by xzcv0Tehd8lQ
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +6 X5bHxOULxFqb 2026-10-30 16/800 2026-11-01 21:19 by xzcv0Tehd8lQ
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 X5bHxOULxFqb 2026-10-30 22/1100 2026-11-01 16:35 by BCdvBAEpt3zm
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 14:48 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 1uaZt2kAZj2S 2026-10-30 17/850 2026-11-01 12:24 by BCdvBAEpt3zm
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-30 16/800 2026-11-01 12:14 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 X5bHxOULxFqb 2026-10-30 18/900 2026-11-01 12:11 by BCdvBAEpt3zm
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 03:59 by YXM1FDCSE8Bn
信息提示
请填处理意见