matlab如何求解待定系数的问题

作者&投稿:应狗 (若有异议请与网页底部的电邮联系)
matlab如何求解待定系数的问题~

【1】把函数变形
>> f=sym('2*x^3+3*x^2+21*x+4-(3*a*x^3+b*x^2+c*x+d)=0')
f =
2*x^3+3*x^2+21*x+4-(3*a*x^3+b*x^2+c*x+d)=0
【2】用collect合并同类项
>>ff=collect(f):
(2-3*a)*x^3+(3-b)*x^2+(21-c)*x+4-d = 0

【3】用maple提取多项式系数,如果很多,可以用循环语句。
>> c3=maple('coeff',ff,x,3)
c3 =2-3*a

>> c1=maple('coeff',ff,x,1)
c1 =21-c

>> c2=maple('coeff',ff,x,2)
c2 =3-b

>>c0=maple('coeff',ff,x,0)
c0 =4-d

补充:
这次变成这样,程序倒是通了,我不是很满意,咱们共同把它整理好,如何?
syms a b c d x
%【1】把函数变形
f=sym('2*x^3+3*x^2+21*x+4-(3*a*x^3+b*x^2+c*x+d)')

N=3;
for i=0:N
temp=maple('coeff',f,x,N-i);
cp(1,i+1)={temp};
end
celldisp(cp);


再补充:这次总算是解决了,就是看起来很笨,不是很理想,凑合着用,当然我相信可以修改的很漂亮。
syms a b c d x
f=sym('2*x^3+3*x^2+21*x+4-(3*a*x^3+b*x^2+c*x+d)')
N=3;
for i=0:N
temp=maple('coeff',f,x,N-i);
temp1(i+1)=temp;
end
cp=temp1
a=solve(cp(1)), b=solve(cp(2)), c=solve(cp(3)), d=solve(cp(4))

运行结果:
a =2/3
b =3
c =21
d =4

假设系数矩阵为A,结果为b 则x=A\b

【1】把函数变形
>>
f=sym('2*x^3+3*x^2+21*x+4-(3*a*x^3+b*x^2+c*x+d)=0')
f
=
2*x^3+3*x^2+21*x+4-(3*a*x^3+b*x^2+c*x+d)=0
【2】用collect合并同类项
>>ff=collect(f):
(2-3*a)*x^3+(3-b)*x^2+(21-c)*x+4-d
=
0
【3】用maple提取多项式系数,如果很多,可以用循环语句。
>>
c3=maple('coeff',ff,x,3)
c3
=2-3*a
>>
c1=maple('coeff',ff,x,1)
c1
=21-c
>>
c2=maple('coeff',ff,x,2)
c2
=3-b
>>c0=maple('coeff',ff,x,0)
c0
=4-d
补充:
这次变成这样,程序倒是通了,我不是很满意,咱们共同把它整理好,如何?
syms
a
b
c
d
x
%【1】把函数变形
f=sym('2*x^3+3*x^2+21*x+4-(3*a*x^3+b*x^2+c*x+d)')
N=3;
for
i=0:N
temp=maple('coeff',f,x,N-i);
cp(1,i+1)={temp};
end
celldisp(cp);
再补充:这次总算是解决了,就是看起来很笨,不是很理想,凑合着用,当然我相信可以修改的很漂亮。
syms
a
b
c
d
x
f=sym('2*x^3+3*x^2+21*x+4-(3*a*x^3+b*x^2+c*x+d)')
N=3;
for
i=0:N
temp=maple('coeff',f,x,N-i);
temp1(i+1)=temp;
end
cp=temp1
a=solve(cp(1)),
b=solve(cp(2)),
c=solve(cp(3)),
d=solve(cp(4))
运行结果:
a
=2/3
b
=3
c
=21
d
=4

If all coefficients are numerical, use function: sym2poly()

>> syms x
>> f=2*x^3+3*x^2+21*x+4;
>> sym2poly(f)

ans =

2 3 21 4

If there is any symbolic coefficient, use function: coeffs()

>> syms x a;
>> f=3*a*x^2+2*x^3+21*x+4;
>> [c,l]=coeffs(f,'x')

c =

[ 4, 21, 2, 3*a]

l =

[ 1, x, x^3, x^2]

Function expand() is sometimes useful.

>> syms x a
>> [c l]=coeffs(expand((a+x)^2),'x')

c =

[ a^2, 2*a, 1]

l =

[ 1, x, x^2]

Matlab,怎么和上次那位同学问的同一个问题啊?
syms a b c d x
f=sym('2*x^3+3*x^2+21*x+4-(3*a*x^3+b*x^2+c*x+d)')
N=3;
for i=0:N
temp=maple('coeff',f,x,N-i);
temp1(i+1)=temp;
end
cp=temp1
a=solve(cp(1)), b=solve(cp(2)), c=solve(cp(3)), d=solve(cp(4))

运行结果:
a =2/3
b =3
c =21
d =4


matlab方程求解
matlab是美国mathworks公司出的商业数学软件,用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境,主要包括MATLAB和Simulink两大部分。但由于其功能的强大,对于初学者显得有些困难,下面,为广大初学者介绍下matlab如何求解方程。matlab方程求解有两种方法,符号求解和数值求解。1....

怎么用Matlab解方程?
我们还可以引入常量k1与k2表达通解。7、对于线性非齐次的解,可以使用linsolve(A,b)。其中A是系数矩阵,b是非齐次项(如果b是多列矩阵,意味着解多个砼系数不同齐次项的方程组)。对于图中的方程组,可已这样求解。特别提示 matlab求解方程的方法有很多,要找到适合自己需要的方法,需要多加练习。

大佬救命这matlab题怎么做?
大佬救命这matlab题怎么做?【问题】常微分方程数值解问题。用预估校正Euler法,求解初值问题。求出步长h=0.1的所有点的值,并绘制图形。【求解方法】用预估校正Euler法来求解,其方法是:第一步,根据y(0)=1边界值

MATLAB求解微分方程问题
解微分方程有两种方式,分别是解析解和数值解。解析解是通过使用dsolve函数进行求解。(1)求解析解 (2)初值问题 (3)边界问题 (4)高阶方程 (5)方程组问题 数值解则是通过使用ode函数进行求解。:% 求解语句 ode45 ode113 ode15s ode23s ode23t ode23tb :% 输出语句 odeplot odeprint ...

matlab求解方程组?
题主给出的复杂的分式方程组,可以用vpasolve()函数得到其数值解。求解方法如下:syms x y eq1=x-(107.1+0.2*(4*y+3*x)*(3.83-107.1)\/(4*y+3*3.83));eq2=y-(83.7+y*(4*y+3*x)*(1.28-83.7)\/(3*x*(3*y+2*1.28)+4*y*(2*y+3*1.28)));[x,y]=vpasolve...

如何用matlab求解线性方程组Ax=0.
现在我们要找出这个正交矩阵Q以及对应的对角阵Λ。A是一个2x2的矩阵,直接计算特征值λ1, λ2,把它作为对角阵Λ的元素。然后把这些特征向量v1, v2分别乘以单位长度组合成列向量q1, q2构成正交矩阵Q。假设我们找到了这样的Q和Λ,有以下关系成立:A = Q * Λ * QT 接下来看第二个问题:求解Ax...

如何用matlab求解线性方程组的符号解?用matlab解符号方程组方法介绍_百 ...
2、上面线性方程组可以表示为矩阵的形式,如下图所示:3、采用矩阵除法的方法求解相信方程组,具体计算代码及结果如下如所示:4、采用solve指令的一般代数方程解法,具体计算代码及结果如下如所示:教程结束,以上就是关于如何用matlab求解线性方程组的符号解?用matlab解符号方程组方法介绍,是不是很简单呢...

用matlab怎么求解微分方程?
Simulink是matlab下基于框图的用于复杂动态系统建模的工具。今以 y"+4y'+4y=cos2t(0≤t≤10),y(0)=1,y'(0)=-2 二阶微分方程为例,用Simulink求解的步骤如下:第一步、分析微分方程的特点确定模块 我们将y"+4y'+4y=cos2t方程改写成y"=cos2t-4y'-4y,从右边可知,式子是由cos2t外加载...

如何用matlab求解微分方程并画图
如何用matlab求解微分方程并画图,可以先用dsolve()或ode()求出其微分方程(组)的解析解或数值解,然后用plot()绘制其图形。例如:解微分方程 y'=y-2t\/y,y(0)=1,0<t<4 1、用dsolve()求解,代码及结果如下 2、用ode45()求解,代码及结果如下 3、当然喽,使用dsolve()或ode(...

matlab如何求解线性方程组?
拓展:1、MATLAB 是美国MathWorks公司出品的商业数学软件,用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境,主要包括MATLAB和Simulink两大部分。2、MATLAB是matrix&laboratory两个词的组合,意为矩阵工厂(矩阵实验室)。是由美国mathworks公司发布的主要面对科学计算、可视化以及交互...

海晏县13935238830: matlab如何求解待定系数的问题 -
门馨坦立: 【1】把函数变形>> f=sym('2*x^3+3*x^2+21*x+4-(3*a*x^3+b*x^2+c*x+d)=0') f =2*x^3+3*x^2+21*x+4-(3*a*x^3+b*x^2+c*x+d)=0 【2】用collect合并同类项>>ff=collect(f):(2-3*a)*x^3+(3-b)*x^2+(21-c)*x+4-d = 0 【3】用maple提取多项式系数,如果很多...

海晏县13935238830: matlab如何求解待定系数的问题
门馨坦立: 【1】把函数变形 &gt;&gt; f=sym('2*x^3+3*x^2+21*x+4-(3*a*x^3+b*x^2+c*x+d)=0') f = 2*x^3+3*x^2+21*x+4-(3*a*x^3+b*x^2+c*x+d)=0 【2】用collect合并同类项 &gt;&gt;ff=collect(f): (2-3*a)*x^3+(3-b)*x^2+(21-c)*x+4-d = 0 【3】用maple提取多项式系...

海晏县13935238830: 如何用matlab求解待定系数法 -
门馨坦立: 代入值,最小2乘法求解

海晏县13935238830: matlab待定系数法求函数系数 -
门馨坦立: 方法是先算出1/x的值,与y对应, x2=1./x; polyfit(x2,y,1)就可以得到最小二乘法拟合的a、b的系数.

海晏县13935238830: Matlab中怎么用最小二乘法求解待定系数?有两个题目求助!如图 -
门馨坦立: 假设系数矩阵为A,结果为b 则x=A\b

海晏县13935238830: matlab怎么待定系数
门馨坦立: 用非线性拟合就可以,我现在也正在弄这个问题,涉及m文件保存什么的,也总是出问题

海晏县13935238830: MATLAB曲线拟合待定系数 -
门馨坦立: 出现这种问题,跟你选择拟合函数有关,跟你选的初值有关,也跟你给出的x,y的对应数值偏差是否太大有关,甚至根本得不到较理想的拟合函数.你能把x,y的对应数值贴出来吗?

海晏县13935238830: 如何用MATLAB解Logistic模型里面的系数 -
门馨坦立: 给你一个例子,如何用MATLAB解Logistic模型里面的系数,希望对你有帮助.实现代码如下:t=[...]; y=[...]; fun=inline('a(1)./(1+exp(a(2)-a(3).*t))','a','t');%Logistic模型 b=[0 0 0];%初值 [a,r,J] = nlinfit(t,y,fun,b);vpa(a,10); x1=t; y1=fun(a,t); R2=corrcoef(y, y1) %R2≈1,可以认为拟合是有效的 [y' y1'] %显示已知值与拟合值

海晏县13935238830: matlab用最小二乘法求一形如y=t/(at+b)(a和b为待定系数)的多项式,使之与下列数据相拟合 -
门馨坦立: 1.使用非线性最小二乘拟合函数lsqcurvefit拟合t=[1 2 3 4 5 6 7 8]; y=[4.00 6.40 8.00 8.80 9.22 9.50 9.70 9.68]; fun=@(b,x)x./(b(1)*x+b(2)); x0=[0.1 0.1]; b=lsqcurvefit(fun,x0,t,y) 结果为: b =0.0811 0.1468 即a=0.0811 b=0.1468 2.绘图 plot(t,y,'ko'); hold on plot(t,fun(b,t),'r-'); legend('Original data','fitted curve ')

海晏县13935238830: 已知曲线x,y两组数据和R,k值,用matlab怎么求系数,需要详细步骤,谢谢! -
门馨坦立: 已知曲线x,y两组数据和R,k值,用matlab求解其系数的步骤如下: 1、数据.x=[...];y=[...]; 2、自定义函数.func=@(a,x)模型表达式.例如:func=@(a,x)a(1)*x^2+a(2)*x+a(3),a(1)、a(2)、a(3)分别对应函数中各个系数. 3、初定初始值.x0=[a10,a20,a30] 4、用nlinfit函数或lsqcurvefit函数求解,模型系数a1,a2,a3. 5、用决定系数R²,判断拟合出来的模型系数a1,a2,a3是否合理.一般要求,R²≈1

本站内容来自于网友发表,不代表本站立场,仅表示其个人看法,不对其真实性、正确性、有效性作任何的担保
相关事宜请发邮件给我们
© 星空见康网