Labfans是一个针对大学生、工程师和科研工作者的技术社区。 论坛首页 | 联系我们(Contact Us)
MATLAB爱好者论坛-LabFans.com
返回   MATLAB爱好者论坛-LabFans.com > 工程数学软件 > MATLAB论坛
MATLAB论坛 一切MATLAB相关问题在此讨论。
回复
 
主题工具 显示模式
旧 2008-11-30, 10:52   #1
ldd00
初级会员
 
注册日期: 2008-11-30
年龄: 41
帖子: 3
声望力: 0
ldd00 正向着好的方向发展
微笑 solve错在哪里?? 求助!!!

solve错在哪里?? 求助!!!
clear all
syms x
A=[1.2 0;0 1.3];
B=A.';
C=[0.5 0;0 0.6];
D=C.';
K=[5 0;0 -3]*B*[exp(4*x-4) 0;0 exp(-x+1)]*D*C*[exp(4*x-4) 0;0 exp(-x+1)]*A;
L=B*[exp(4*x-4) 0;0 exp(-x+1)]*D*C*[exp(4*x-4) 0;0 exp(-x+1)]*A*[5 0;0 -3];
E=K+L;
y=[10 10]*[exp(x) 0;0 exp(-2*x)]*E*[exp(x) 0;0 exp(-2*x)]*[10;10];
s=solve(y,x)
s =
5/8-1/16*log(500/507)
这个程序运行后的结果很正确 可是如果我把矩阵A C改成变量的形式 为什么运行就不对了呢 ???
clear all
syms a b c d m n p q x
A=[a b;c d];
B=A.';
C=[m n;p q];
D=C.';
K=[5 0;0 -3]*B*[exp(4*x-4) 0;0 exp(-x+1)]*D*C*[exp(4*x-4) 0;0 exp(-x+1)]*A;
L=B*[exp(4*x-4) 0;0 exp(-x+1)]*D*C*[exp(4*x-4) 0;0 exp(-x+1)]*A*[5 0;0 -3];
E=K+L;
y=[10 10]*[exp(x) 0;0 exp(-2*x)]*E*[exp(x) 0;0 exp(-2*x)]*[10;10];
s=solve(y,x)
运行后出现错误为:
??? Error using ==> solve
Unable to find closed form solution.
Error in ==> sym.solve at 49
[varargout{1:max(1,nargout)}] = solve(S{:});

请高手帮忙啊!!! 想要实现y=0时 x的值 用a b c d m n p q 表示形式
ldd00 当前离线   回复时引用此帖
旧 2008-11-30, 11:39   #2
debateshang
初级会员
 
注册日期: 2008-11-30
年龄: 40
帖子: 20
声望力: 17
debateshang 正向着好的方向发展
默认 回复: solve错在哪里?? 求助!!!

可能没有显示解我把y值直接带入solve中求解出错信息是
solve('10*(10*exp(x)*((5*m*exp(4*x-4)*a+5*n*exp(-x+1)*c)*(m*exp(4*x-4)*a+n*exp(-x+1)*c)+(5*p*exp(4*x-4)*a+5*q*exp(-x+1)*c)*(p*exp(4*x-4)*a+q*exp(-x+1)*c)+5*(m*exp(4*x-4)*a+n*exp(-x+1)*c)^2+5*(p*exp(4*x-4)*a+q*exp(-x+1)*c)^2)+10*exp(-2*x)*((-3*m*exp(4*x-4)*b-3*n*exp(-x+1)*d)*(m*exp(4*x-4)*a+n*exp(-x+1)*c)+(-3*p*exp(4*x-4)*b-3*q*exp(-x+1)*d)*(p*exp(4*x-4)*a+q*exp(-x+1)*c)+5*(m*exp(4*x-4)*a+n*exp(-x+1)*c)*(m*exp(4*x-4)*b+n*exp(-x+1)*d)+5*(p*exp(4*x-4)*a+q*exp(-x+1)*c)*(p*exp(4*x-4)*b+q*exp(-x+1)*d)))*exp(x)+10*(10*exp(x)*((5*m*exp(4*x-4)*a+5*n*exp(-x+1)*c)*(m*exp(4*x-4)*b+n*exp(-x+1)*d)+(5*p*exp(4*x-4)*a+5*q*exp(-x+1)*c)*(p*exp(4*x-4)*b+q*exp(-x+1)*d)-3*(m*exp(4*x-4)*a+n*exp(-x+1)*c)*(m*exp(4*x-4)*b+n*exp(-x+1)*d)-3*(p*exp(4*x-4)*a+q*exp(-x+1)*c)*(p*exp(4*x-4)*b+q*exp(-x+1)*d))+10*exp(-2*x)*((-3*m*exp(4*x-4)*b-3*n*exp(-x+1)*d)*(m*exp(4*x-4)*b+n*exp(-x+1)*d)+(-3*p*exp(4*x-4)*b-3*q*exp(-x+1)*d)*(p*exp(4*x-4)*b+q*exp(-x+1)*d)-3*(m*exp(4*x-4)*b+n*exp(-x+1)*d)^2-3*(p*exp(4*x-4)*b+q*exp(-x+1)*d)^2))*exp(-2*x)')
??? Error using ==> solve at 162
Unable to find closed form solution
debateshang 当前离线   回复时引用此帖
旧 2008-11-30, 11:40   #3
debateshang
初级会员
 
注册日期: 2008-11-30
年龄: 40
帖子: 20
声望力: 17
debateshang 正向着好的方向发展
默认 回复: solve错在哪里?? 求助!!!

clear all
clc
syms a b c d m n p q
A=[a b;c d];
C=[m n;p q];
syms x
% A=[1.2 0;0 1.3];
% C=[0.5 0;0 0.6];
k1=C*[exp(4*x-4) 0;0 exp(-x+1)]*A;
K=[5 0;0 -3]*k1.'*k1;
L=k1.'*k1*[5 0;0 -3];
E=K+L;
F=[exp(x) 0;0 exp(-2*x)]*[10;10];
y=F.'*E*F;
% solve(y)
debateshang 当前离线   回复时引用此帖
旧 2008-11-30, 15:15   #4
ldd00
初级会员
 
注册日期: 2008-11-30
年龄: 41
帖子: 3
声望力: 0
ldd00 正向着好的方向发展
默认 回复: solve错在哪里?? 求助!!!

你的这个还是不行啊???
ldd00 当前离线   回复时引用此帖
旧 2008-12-01, 11:22   #5
debateshang
初级会员
 
注册日期: 2008-11-30
年龄: 40
帖子: 20
声望力: 17
debateshang 正向着好的方向发展
默认 回复: solve错在哪里?? 求助!!!

我的意思是你这种设法有问题,可能matlab中无法处理
solve('a*b-d*exp(x-3)+exp(-4*x+2)=0','x');
??? Error using ==> solve at 162Unable to find closed form solution.
解这个方程matlab都解不出解析解,字母太多在matlab里可能就没法讨论了
这就像五次方程没有一般系数解
solve('a*x^5+b*x^4+c*x^3+d*x^2+m*x+n','x');
??? Error using ==> solve at 162
Unable to find closed form solution.
但是
solve('3*x^5+1*x^4+2*x^3+4*x^2+2*x+1')

ans =

.50436452606242968189338746040867+1.0358328721587116228614551628226*i
-.23178710110867900315683730233415+.48181007825073643046183194749188*i
-.87848818324083469080643364948239
-.23178710110867900315683730233415-.48181007825073643046183194749188*i
.50436452606242968189338746040867-1.0358328721587116228614551628226*i
就能解
debateshang 当前离线   回复时引用此帖
回复


发帖规则
不可以发表新主题
不可以发表回复
不可以上传附件
不可以编辑自己的帖子

启用 BB 代码
论坛启用 表情符号
论坛启用 [IMG] 代码
论坛禁用 HTML 代码


相似的主题
主题 主题作者 版面 回复 最后发表
[求助]采样simulink数据到工作空间的问题 vsunl MATLAB论坛 0 2009-02-23 21:31
matlab可以建气溶胶荧光光谱数据库吗 rongerlovesj MATLAB论坛 0 2009-02-06 15:20
[求助]matlab rgb图像滤波和robert、Laplacian锐化 doudouzhjm MATLAB论坛 0 2008-12-26 19:31
[求助]matlab怎么跟网络浏览器绑定? lhongzhong MATLAB论坛 1 2008-11-26 17:43
Simplify Instrument Communication: Using Instrument Drivers with the Instrument Contr TechnicalArticles MATLAB技术文章 0 2008-01-06 16:32


所有时间均为北京时间。现在的时间是 01:32


Powered by vBulletin
版权所有 ©2000 - 2025,Jelsoft Enterprises Ltd.