| 查看: 727 | 回复: 2 | ||
[求助]
用MATLAB做符号运算，但是没有explicited solution
|
|
syms r1 r2 x y z a b c1 c2 Vzx Vzy Vzz Vxx Vxy Vyy Bx By Bz Mz Mx My u0 f1 f2 f3 f4 f5 f6 f7 f8 Bxx Bxy Bxz Byy Byz; r1=sqrt((x-a)^2+(y-b)^2+(z-c1)^2); r2=sqrt((x-a)^2+(y-b)^2+(z-c2)^2); Vzx=(x-a)*(1/r2^3-1/r1^3); Vzy=(y-b)*(1/r2^3-1/r1^3); Vzz=(c2-z)/r2^3-(c1-z)/r1^3; Vxx=((y-b)^2-(x-a)^2)/((y-b)^2+(x-a)^2)^2*((c1-z)/r1-(c2-z)/r2)-(x-a)^2/((y-b)^2+(x-a)^2)*((c1-z)/r2^3-(c2-z)/r1^3); Vxy=2*(x-a)*(y-b)/((y-b)^2+(x-a)^2)^2*((z-c1)/r1-(z-c2)/r2)-(x-a)*(y-b)/((y-b)^2+(x-a)^2)*((z-c1)/r2^3-(z-c2)/r1^3); Vyy=((x-a)^2-(y-b)^2)/((y-b)^2+(x-a)^2)^2*((c1-z)/r1-(c2-z)/r2)-(y-b)^2/((y-b)^2+(x-a)^2)*((c1-z)/r2^3-(c2-z)/r1^3); Bx=u0/(4*pi)*(Mx*Vxx+My*Vxy+Mz*Vzx); By=u0/(4*pi)*(Mx*Vxy+My*Vyy+Mz*Vzy); Bz=u0/(4*pi)*(Mx*Vzx+My*Vzy+Mz*Vzz); S=solve(Bx==u0/(4*pi)*(Mx*Vxx+My*Vxy+Mz*Vzx),By==u0/(4*pi)*(Mx*Vxy+My*Vyy+Mz*Vzy),Bz==u0/(4*pi)*(Mx*Vzx+My*Vzy+Mz*Vzz),Bxx==diff(Bx,x),... Bxy==diff(Bx,y),Bxz==diff(Bx,z),Byy==diff(By,y),a,b,c1,c2,Mx,My,Mz); 拜托一些大牛们，帮我看看 |
回复此楼» 猜你喜欢
浙大物院表面物理课题组长期招收普博及直博学生
已经有259人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
物理学I论文润色/翻译怎么收费?
已经有55人回复
-
25届硕士申请，材料计算方向
已经有2人回复
-
课题组招1-2名硕士
已经有1人回复
-
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
-
斯坦福大学与Elsevier世界前2%学者新排名查询
已经有1人回复
-
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有19人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
apl投稿状态
已经有2人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有0人回复
erichan911
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 3741.5
- 散金: 15
- 红花: 1
- 帖子: 202
- 在线: 46.4小时
- 虫号: 1913102
- 注册: 2012-07-27
- 性别: GG
- 专业: 工程热物理与能源利用
2楼2018-04-21 19:17:35
独孤神宇
版主 (知名作家)
- 应助: 490 (硕士)
- 贵宾: 0.008
- 金币: 31030.3
- 散金: 802
- 红花: 122
- 沙发: 1
- 帖子: 5603
- 在线: 856.8小时
- 虫号: 3522474
- 注册: 2014-11-06
- 性别: GG
- 专业: 机械动力学
- 管辖: 计算模拟
3楼2018-04-21 19:21:55
50