bbtw
2010-10-13, 20:35
求助:怎么样进行函数化简
求助:算出来的结果长的令人发指,怎么样进行简化。
代码如下:
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)
求助:算出来的结果长的令人发指,怎么样进行简化。
代码如下:
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)