TBE_Legend 发表于 2009-11-27 21:41:52

Maple版联合Matlab版,大家来做题拿分啦,加分,大大地有 ^_^

本帖最后由 TBE_Legend 于 2009-11-27 21:59 编辑

http://forum.simwe.com/thread-908305-1-1.html

maplelab 发表于 2009-11-28 10:07:19

1# TBE_Legend

restart; assign(p = 0.20e-1, q = .219, m = 10369.6, C1 = 800, C2 = 2, C3 = 6, C4 = 8, k = 3, g = 4, Txing = 8, T2 = 12); A := m*p*(p+q)^2*exp(-(p+q)*u)/(p+q*exp(-(p+q)*u))^2; B := m*p*(p+q)^2*exp(-(p+q)*Txing)/(p+q*exp(-(p+q)*Txing))^2; W1 := -C2*(int((int(A, u = t .. T2))/(g*(p+q)), t = t2 .. T2)); W2 := C4*B*(1-k*(p+q)); W3 := C3*(k*(p+q)*(int(A, u = T1 .. Txing))+B*k*(p+q)*t2-k*(p+q)*Txing); eq1 := W1+W2+W3 = 0; eq2 := C2*T1^2/(2*g*(p+q))-k*(p+q)*C3*(t2-T1)-C4*(1-k*(p+q)) = 0; fsolve({eq1, eq2});



这样的题目不是难,而是输入麻烦.建议多搞些输入不麻烦的.
当然也可以先将题目的代码贴出.再让回复的人自己想办法解决题目.

maplelab 发表于 2009-11-28 10:21:04

本帖最后由 maplelab 于 2009-11-28 10:22 编辑

将fsolve改为solve更好:restart; assign(p = 0.20e-1, q = .219, m = 10369.6, C1 = 800, C2 = 2, C3 = 6, C4 = 8, k = 3, g = 4, Txing = 8, T2 = 12); A := m*p*(p+q)^2*exp(-(p+q)*u)/(p+q*exp(-(p+q)*u))^2; B := m*p*(p+q)^2*exp(-(p+q)*Txing)/(p+q*exp(-(p+q)*Txing))^2; W1 := -C2*(int((int(A, u = t .. T2))/(g*(p+q)), t = t2 .. T2)); W2 := C4*B*(1-k*(p+q)); W3 := C3*(k*(p+q)*(int(A, u = T1 .. Txing))+B*k*(p+q)*t2-k*(p+q)*Txing); eq1 := W1+W2+W3 = 0; eq2 := C2*T1^2/(2*g*(p+q))-k*(p+q)*C3*(t2-T1)-C4*(1-k*(p+q)) = 0; solve({eq1, eq2});

wanglu 发表于 2009-11-28 11:55:34

maplelab 能否看看这个帖子,虽然帖子长,但问题输入不麻烦,谢谢!
http://forum.simwe.com/thread-908387-1-1.html
页: [1]
查看完整版本: Maple版联合Matlab版,大家来做题拿分啦,加分,大大地有 ^_^