matlab解高阶微分方程怎么解混合一元高次方程?见方程

MATLAB怎么解一元二次方程?
首先一元二次方程需要用Matla解吗?呵呵,初中的知识吧
如果用Matla可以用fsolve函数或者Solve函数
具体使用方法的话键入help solve就可以啦
一般来说解一元二次方程有以下做法:
比如解x^2-4x+3=0
解法一:因式分解法
x^2-4x+3=(x-3)(x-1)=0
于是x=3或x=1
解法二:配方...
一元二次方程的求根公式是
解:设第二次购书x本,则第一次购书(x-10)本
由题意,得[100/(x-10)]+0.5=150/x
整理得x²-110x+3000=0,解得x1...
要表示指数: e^(n)写成 exp(n)对数:
log(n)底数e 对数
希望这次对了。
答: AMD四核640 (685元参考北京)技嘉880 (799元参考北京)2G内存(DDR)(80元参考金士顿2GB DDR3 G硬盘 (238元参...
答: 可以用路游器
答: 炒股,买卖期货,开网络店铺。呵呵。。
大家还关注
Copyright &
Corporation, All Rights Reserved
确定举报此问题
举报原因(必选):
广告或垃圾信息
激进时政或意识形态话题
不雅词句或人身攻击
侵犯他人隐私
其它违法和不良信息
报告,这不是个问题
报告原因(必选):
这不是个问题
这个问题分类似乎错了
这个不是我熟悉的地区当前位置:
&用matlab解一个含参数的一元三次方程,高手帮忙看看这个程序有什么问题,谢谢!
用matlab解一个含参数的一元三次方程,高手帮忙看看这个程序有什么问题,谢谢!
作者 xmuxiaoyu
用matlab解一个含参数的一元三次方程:
C3*P(z)^3+C1*P(z)+D(z)-P_av=0&&
z是一线性数组,C3,C1常数,P_av是P在不同z的根的平均值。
思路是首先假设一个初始值P_av0, 然后解得P的一组根,取平均值和P_av0比较,如果收敛,则为所求的根,否则把新的平均值代回方程,再求根,直到收敛。
程序如下:
function P_real=P_distribution(a,b,e,q,d,N,m,z,P_r,Dsc)
format long
% define parameters
a=-6.36e7;& && && && && && && && && && && && && && && && && &
b=1.38e9;& && && && && && && && && && && && && && && && && &
e=8.85e-11;& && && && && && && && && && && && && && && && && && &
q=1.61e-19;& && && && && && && && && && && && && && && && &&&
d=1e-7;& && && && && && && && && && && && && && && && && &&&
% coeffecients
% start value
P_av0=0.2;& && && & % initial value for the P average
P_real=zeros(1,11);
%start the caculation
& &while 1
& && & z=m*1e-8;
& && & Dsc=-q*N*(d/2-z);
& && & c0=Dsc-P_av0;
& && & P_eq=[c3,0,c1,c0];
& && & P_r=roots(P_eq);& &%get the roots
& && & for k=1:length(P_r)
& && && &&&P_rr=isreal(P_r(k));& &%get the real roots
& && & end
& && & pp=P_r(P_rr);
& && & P_real(1,m+1)=abs(pp);& && & %assign the real roots to the matrix
& && & m=m+1;& && && &&&%solve the equation again for another value
& && & if m&10
& && && &&&break
& && & end
& && & P_av1=mean(P_real);& & %get the average value
& && & if abs(P_av1-P_av0)&1e-11& &%less than the tolerance
& && && &&&break
& && & else
& && && &&&P_av0=P_av1;
& && & end
程序在matlab中无法运行,不知道问题出在哪里,高手帮忙看看,多谢了!
CODE:if abs(P_av1-P_av0)&1e-11& &%less than the tolerance中精度1e-11设置的太高了,关键与前面的rootsCODE:P_r=roots(P_eq);& &%get the roots相比精度太高了,
建议把root换成fzero,因为fzero可以直接控制误差的
引用回帖:: Originally posted by wangww2011 at
if abs(P_av1-P_av0)&1e-11& &%less than the tolerance
中精度1e-11设置的太高了,关键与前面的roots
P_r=roots(P_eq);& &%get the roots
相比精度太高了,
建议把root换成fzero,因为fzero可以直接 ... 尝试了降低精度,和用fzero,
但是还是无法在matlab中运行,matlab也没有报错,不知道问题出在哪里
引用回帖:: Originally posted by csgt0 at
end m&10已经给出条件了,
所以并不需要m=0呀
引用回帖:: Originally posted by xmuxiaoyu at
m&10已经给出条件了,
所以并不需要m=0呀... 1.是可以运行的
2.是因为你原意是将每次preal赋给下一次的P_av1,但是你的程序是把所有的preal均值赋给pav,因为你的m一直增大,preal就越来越大。降低精度就可以看到很快就算完了。
引用回帖:: Originally posted by xmuxiaoyu at
m&10已经给出条件了,
所以并不需要m=0呀... 楼上说的对&&你的m是一直增大的 应该加上m=0才行的
你可以在两个while循环之间加上m=0,不需要改精度
引用回帖:: Originally posted by csgt0 at
1.是可以运行的
2.是因为你原意是将每次preal赋给下一次的P_av1,但是你的程序是把所有的preal均值赋给pav,因为你的m一直增大,preal就越来越大。降低精度就可以看到很快就算完了。... 我的意思是把P_real数组的平均值赋给P_av1。
m增大,只是求不同z值的根,也就是一组方程的根。一共11个值。
P_real不会越来越大的,只是给数组里的每一个数赋予一个根。
或者我没有理解你的意思,这个程序能在你的机上运行吗?十分感谢,
24小时热帖
下载小木虫APP
与700万科研达人随时交流为什么用Matlab有的5次方程能得到根式解,
如x^5+5x/4-1=0
就能得到根式解
而x^5-x-1=0就只有近似值??
  解答:
  因为根据Galois 理论:一般一元五次方程式?]有根式解(即不是所有的一元五次方程式都有根式解)。且可得出一元五次方程式有根式解的充分必要条件(比较复杂,这里就不一一叙述)。
所有代数方程的解都是根式形式的,整数解、有理数解、无理数解、虚数解实际上都是由根式得到的。
你所说的根式形式的解,应该是指代数方程的无理数解是不是一定可以表示为...
对不起大家。
理论上应该有三个值
应该是覆盖了
其三个值分别为:2*e^(-pi*i/3)
2*e^(pi*i/3)
2*e^(pi*i)=-2
下面的程序可以得到...
一元五次方程,至少有一个实数根,(虚数根成对),有且只有五个根,1.求实数根:1)因式分解,2)算系数==》系数和等于0,有+1根,偶数项系数和等于奇数项系数和...
5x+2.4x=8.88
得出7.4X=8.88
X=8.88除以7.4=1.2
答: 1736X(1+10%)=1909.6(元)
答: x->0:lim(1+x)^(-1/x)
=1/[x->0:lim(1+x)^(1/x)
x->∞:limxsin(1/x)
=1/x->0:lim[...
答: 数学:甲数、乙数与丙数的和是1400,甲数是乙数的2倍,丙数是乙数的二分之一,求甲、乙、丙各多少?
答: 对于那些有志于穷尽数学奥秘的学生,他总是循循善诱地予以启发和教育,而对于那些急功近利、在学习上不肯刻苦钻研的人,则毫不客气地予以批评
大家还关注
Copyright &
Corporation, All Rights Reserved
确定举报此问题
举报原因(必选):
广告或垃圾信息
激进时政或意识形态话题
不雅词句或人身攻击
侵犯他人隐私
其它违法和不良信息
报告,这不是个问题
报告原因(必选):
这不是个问题
这个问题分类似乎错了
这个不是我熟悉的地区&>&MATLAB解一元多次方程程序
MATLAB解一元多次方程程序
上传大小:495B
MATLAB解一元多次方程程序MATLAB解一元多次方程程序MATLAB解一元多次方程程序
综合评分:4(14位用户评分)
下载个数:
{%username%}回复{%com_username%}{%time%}\
/*点击出现回复框*/
$(".respond_btn").on("click", function (e) {
$(this).parents(".rightLi").children(".respond_box").show();
e.stopPropagation();
$(".cancel_res").on("click", function (e) {
$(this).parents(".res_b").siblings(".res_area").val("");
$(this).parents(".respond_box").hide();
e.stopPropagation();
/*删除评论*/
$(".del_comment_c").on("click", function (e) {
var id = $(e.target).attr("id");
$.getJSON('/index.php/comment/do_invalid/' + id,
function (data) {
if (data.succ == 1) {
$(e.target).parents(".conLi").remove();
alert(data.msg);
$(".res_btn").click(function (e) {
var parentWrap = $(this).parents(".respond_box"),
q = parentWrap.find(".form1").serializeArray(),
resStr = $.trim(parentWrap.find(".res_area_r").val());
console.log(q);
//var res_area_r = $.trim($(".res_area_r").val());
if (resStr == '') {
$(".res_text").css({color: "red"});
$.post("/index.php/comment/do_comment_reply/", q,
function (data) {
if (data.succ == 1) {
var $target,
evt = e || window.
$target = $(evt.target || evt.srcElement);
var $dd = $target.parents('dd');
var $wrapReply = $dd.find('.respond_box');
console.log($wrapReply);
//var mess = $(".res_area_r").val();
var mess = resS
var str = str.replace(/{%header%}/g, data.header)
.replace(/{%href%}/g, 'http://' + window.location.host + '/user/' + data.username)
.replace(/{%username%}/g, data.username)
.replace(/{%com_username%}/g, data.com_username)
.replace(/{%time%}/g, data.time)
.replace(/{%id%}/g, data.id)
.replace(/{%mess%}/g, mess);
$dd.after(str);
$(".respond_box").hide();
$(".res_area_r").val("");
$(".res_area").val("");
$wrapReply.hide();
alert(data.msg);
}, "json");
/*删除回复*/
$(".rightLi").on("click", '.del_comment_r', function (e) {
var id = $(e.target).attr("id");
$.getJSON('/index.php/comment/do_comment_del/' + id,
function (data) {
if (data.succ == 1) {
$(e.target).parent().parent().parent().parent().parent().remove();
$(e.target).parents('.res_list').remove()
alert(data.msg);
//填充回复
function KeyP(v) {
var parentWrap = $(v).parents(".respond_box");
parentWrap.find(".res_area_r").val($.trim(parentWrap.find(".res_area").val()));
评论共有5条
这个代码不错,求解的方法值得学习
需要修改才能运行!
还得琢磨一下,楼主能给加个算法说明就更好了,谢谢
kingsudung
综合评分:
积分/C币:3
综合评分:
积分/C币:3
综合评分:
积分/C币:3
综合评分:
积分/C币:3
综合评分:
积分/C币:10
VIP会员动态
CSDN下载频道资源及相关规则调整公告V11.10
下载频道用户反馈专区
下载频道积分规则调整V1710.18
spring mvc+mybatis+mysql+maven+bootstrap 整合实现增删查改简单实例.zip
资源所需积分/C币
当前拥有积分
当前拥有C币
输入下载码
为了良好体验,不建议使用迅雷下载
MATLAB解一元多次方程程序
会员到期时间:
剩余下载个数:
剩余积分:0
为了良好体验,不建议使用迅雷下载
积分不足!
资源所需积分/C币
当前拥有积分
您可以选择
程序员的必选
绿色安全资源
资源所需积分/C币
当前拥有积分
当前拥有C币
(仅够下载10个资源)
为了良好体验,不建议使用迅雷下载
资源所需积分/C币
当前拥有积分
当前拥有C币
为了良好体验,不建议使用迅雷下载
资源所需积分/C币
当前拥有积分
当前拥有C币
您的积分不足,将扣除 10 C币
为了良好体验,不建议使用迅雷下载
你当前的下载分为234。
你还不是VIP会员
开通VIP会员权限,免积分下载
你下载资源过于频繁,请输入验证码
您因违反CSDN下载频道规则而被锁定帐户,如有疑问,请联络:!
若举报审核通过,可返还被扣除的积分
被举报人:
举报的资源分:
请选择类型
资源无法下载
资源无法使用
标题与实际内容不符
含有危害国家安全内容
含有反动色情等内容
含广告内容
版权问题,侵犯个人或公司的版权
*详细原因:
MATLAB解一元多次方程程序}

我要回帖

更多关于 matlab解高阶方程 的文章

更多推荐

版权声明:文章内容来源于网络,版权归原作者所有,如有侵权请点击这里与我们联系,我们将及时删除。

点击添加站长微信