查看单个帖子
旧 2010-10-13, 20:35   #1
bbtw
初级会员
 
注册日期: 2008-08-27
年龄: 42
帖子: 1
声望力: 0
bbtw 正向着好的方向发展
默认 求助,怎样对输出结果进行化简,谢谢了

求助:怎么样进行函数化简
求助:算出来的结果长的令人发指,怎么样进行简化。
代码如下:
function qiuhe
clear;
clc;
syms x y z
j0 = 1.11e-3;
u0 = 18.53e3;
eplong0 = 8.85e-12;
k = 2.1e-4;
b = 0.09;
h = 0.082;
c = 0.075;
alfa = 2*j0/(pi*k*eplong0);
beta = u0^2/b^2;
for i=1:2
for j=1:22
E0(i,j) = sqrt([y-0.0425-(i-1)*0.15].^2+[z-0.0625-(j-1)*0.075].^2).*[alfa.*x+beta.*(1-x/h)]/c
+ 2*j0*h/(k*eplong0).*(h/x+1/pi-1).*(1-sqrt([y-0.0425-(i-1)*0.15]^2+[z-0.0625-(j-1)*0.075]^2)/c);
E0(i,j) = vpa(E0(i,j),2);
end
for k = 1:21
E1(i,k) = sqrt([y-0.1175-(i-1)*0.15]^2+[z-0.1-(k-1)*0.075]^2)*[alfa*x+beta*(1-x/h)]/c
+ 2*j0*h/(k*eplong0)*(h/x+1/pi-1)*(1-sqrt([y-0.1175-(i-1)*0.15]^2+[z-0.1-(k-1)*0.075]^2)/c);
E1(i,k) = vpa(E1(i,k),2);
end
end
E0 = sum(E0();
E1 = sum(E1();
E2 = E0+E1;
E = sqrt(E2)
bbtw 当前离线   回复时引用此帖