声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

12
返回列表 发新帖
楼主: zhlong

[编程技巧] 多项式求解问题

[复制链接]
 楼主| 发表于 2007-8-3 10:20 | 显示全部楼层
还要麻烦一下各位。

附件p是系数(^_^,这次里面有四个参数了)(ps:如下p(1)/4、 P(2)/3、 p(3)/2、 p(4)才是多项式系数)

x=[1:20000];
Fx=p(1)/4*x.^4+p(2)/3*x.^3+p(3)/2*x.^2+p(4)*x;
plot(Fx)

画出多项式如图所示:

1.jpg

我原来的问题是:如何由给定的Fx值,如Fx=[0:2*pi/20:2*pi*3000]求出相应的x值,并且保证随着Fx增长,x也单调增长。

下面按照花如月版主给出的解决办法进行:
syms x a;
Fx=p(1)/4*x.^4+p(2)/3*x.^3+p(3)/2*x.^2+p(4)*x-a
k=solve(Fx)

for i=1:length(k)
    subplot(2,2,i);
    fplot(inline(char(k(i))),[0 4000*pi]);
end

得出结果如下图:
2.jpg

我的问题是: 第二幅图的结果基本为直线,这样的结果对吗?

p.mat

209 Bytes, 下载次数: 3

回复 支持 反对
分享到:

使用道具 举报

发表于 2007-8-3 11:31 | 显示全部楼层
mat文件个人觉得实在不方便,龙大侠,再弄就直接给出系数的值。程序应该没有问题,图画成那个样子应该是数据范围取得太小吧。所以没有看到全貌
fplot(inline(char(k(i))),[- 4000*pi 4000*pi]);试试看
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-5-18 07:12 , Processed in 0.047934 second(s), 19 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表