应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » matlab 用solve解决不了的参数方程 怎么办？求大牛们指点指点哈
南方科技大学公共卫生及应急管理学院2025级博士研究生招生报考通知
51/1返回列表
查看: 943  |  回复: 4
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

pottery

金虫 (正式写手)

[求助] matlab 用solve解决不了的参数方程 怎么办？求大牛们指点指点哈已有2人参与

指令如下 5个方程组y1-5，求3个参数B2 B3 A3；求帮忙解决哈
syms B2 A3 B3 k1 k2 k3 a b z y1 y2 y3 y4 me2 me3
y1=int('(B2*sin(k2*z))^2 ',z,0,a)+int('(A3*exp(k3*z)+B3*exp(-k3*z))^2',z,a,b)-1;
y2=B2*sin(k2*a)-A3*exp(k3*a)-B3*exp(k3*a);
y3=B2*k2*cos(k2*a)/me2-(A3*k3*exp(k3*a)-B3*k3*exp(-k3*a))/me3;
y4=A3*exp(k3*b)+B3*exp(-k3*b);
y5=k3*sin(k2*a)/cos(k2*a)*me2+k2*me3*(exp(k3*b-k3*a)-exp(-k3*b+k3*a))/(exp(k3*b-k3*a)+exp(-k3*b+k3*a));
S=solve('y1','y2','y3','y4','y5','B2','B3','A3');
S.B2,S.B3,S.A3
先谢谢哈
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
学者一名
1楼2014-03-12 11:25:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pottery

金虫 (正式写手)
不太懂 好像不行吧。。
赞一下 回复此楼
学者一名
2楼2014-03-12 11:37:35
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hppdyx

木虫 (知名作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
pottery: 金币+10, ★有帮助 2014-03-15 18:16:36
CODE:
function pottery
%%
syms B2 A3 B3 k1 k2 k3 a b z y1 y2 y3 y4 me2 me3;
y1=int((B2*sin(k2*z))^2,z,0,a)+int((A3*exp(k3*z)+B3*exp(-k3*z))^2,z,a,b)-1;
y2=B2*sin(k2*a)-A3*exp(k3*a)-B3*exp(k3*a);
y3=B2*k2*cos(k2*a)/me2-(A3*k3*exp(k3*a)-B3*k3*exp(-k3*a))/me3;
% y4=A3*exp(k3*b)+B3*exp(-k3*b);
% y5=k3*sin(k2*a)/cos(k2*a)*me2+k2*me3*(exp(k3*b-k3*a)-exp(-k3*b+k3*a))/(exp(k3*b-k3*a)+exp(-k3*b+k3*a));
S = solve(y1==0, y2==0, y3==0, B2, B3, A3);
S = [S.B2; S.B3; S.A3];
display(S);
% S.B2; S.B3; S.A3;

里面有两个问题：1. y5里面不含未知数

2. y1, y2, y3, y4求解三个未知数，方程多了。。。我上面试的是前三个方程求解的，可以得到结果。但是使用不同的三个方程求得的结果不同（已试）。
赞一下 回复此楼
不以风骚惊天下，但求淫荡动世人
3楼2014-03-15 10:29:56
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pottery

金虫 (正式写手)
引用回帖:
3楼: Originally posted by hppdyx at 2014-03-15 10:29:56
function pottery
%%
syms B2 A3 B3 k1 k2 k3 a b z y1 y2 y3 y4 me2 me3;
y1=int((B2*sin(k2*z))^2,z,0,a)+int((A3*exp(k3*z)+B3*exp(-k3*z))^2,z,a,b)-1;
y2=B2*sin(k2*a)-A3*exp(k3*a)-B3*exp(k3*a);
y3=B2 ...

syms B2 A3 B3 k1 k2 k3 a b z y1 y2 y3 y4 me2 me3
y1=int('(B2*sin(k2*z))^2 ',z,0,a)+int('(A3*exp(k3*z)+B3*exp(-k3*z))^2',z,a,b)-1;
y2=B2*sin(k2*a)-A3*exp(k3*a)-B3*exp(-k3*a);
y3=B2*k2*cos(k2*a)/me2-(A3*k3*exp(k3*a)-B3*k3*exp(-k3*a))/me3;
y4=A3*exp(k3*b)+B3*exp(-k3*b);
y5=k3*sin(k2*a)/cos(k2*a)*me2+k2*me3*(exp(k3*b-k3*a)-exp(-k3*b+k3*a))/(exp(k3*b-k3*a)+exp(-k3*b+k3*a));
S=solve('y1','y2','y3','y4','B2','B3','A3');
S.B2,S.B3,S.A3
才发现y2漏了个负号，谢谢你。y5可以不用。确实 solve函数好像不能解决问题。
赞一下 回复此楼
学者一名
4楼2014-03-15 18:15:54
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xmjiao

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
pottery: 金币+20, ★有帮助 2014-03-17 22:08:30
领金币来了，哈哈哈哈。。。
赞一下 回复此楼
心有多大，舞台就有多大！
5楼2014-03-17 22:08:09
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 pottery 的主题更新
51/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭