你好！非常感谢你上次帮我解答的问题。我现在有一个类似的问题，麻烦你再帮我编一下好吗？
这个是我毕设的问题。我需要将E的这个式子画出其关系图形，并根据e1取不同值画出不同情况下的图，如e1=9时，e1=3时。下图是老师画出来的。麻烦你帮我编一下，万分感谢！把我所有的财富值都送给你啦！
谢谢！
请帮忙写一段matlab的程序，谢谢！
[email protected]
下面是ε1=3的情况，但结果很不同。
clear;
clc;
close all;
ep0=1;
ep1=3;
a=1; %是不是?
k1=(ep0-ep1)/(ep0+ep1);
k2=2*ep0/(ep0+ep1);
k3=2*ep1/(ep0+ep1);
q=1e-7;
z=1;
Eout=zeros(1,20);%取结果为20个数据，更精确的话多取
b=linspace(0.01,0.4,20);%b为0.01到0.04，根据需要改
t0=q/(4*pi*ep0);
for i=1:20
N=fix(a/2/b(i)+1);
n=0:N;
t1=sum(k1.^(2*n)./(z-2*n*b(i)).^2);
n=1:N;
t2=sum(k1.^(2*n-1)./(z-2*a+2*(n-1)*b(i)).^2);
n=N:N+9;
t3=sum(k3*k1.^(2*n)./(z-2*b(i)*n).^2);
n=1:10;
t4=sum(k2*k1.^(2*n-1)./(z-2*a-2*n*b(i)).^2);
Eout(i)=t0.*(t1+t2+t3-t4);
end
plot(b,Eout)追问
我给你发到邮箱了，谢谢！页数范围在53，我要仿真的是公式3.4.7
追答我不知道如何得到和上面一样的图，但是从论文看，我知道要得到什么结果，所有资料我发到你邮箱。
来自：求助得到的回答a=10,q(t)=0.0000001
追答这样的话，E=常数，要做哪两个变量间的关系图
追问a值其实是变化的，为了简化，我暂时去了a=10。a值是由远及近或由近及远的变化。
要求的关系图是阻隔介质厚度b和电场强度E的关系；以及探测距离a与电场强度E关系图。是两个图。谢谢！
或者说，是一个三维图，取了不同角度，得到两张关系图
关系式真的好复杂啊，累加和的总个数N还与a和b有关，累加和用什么方法来算，我用符号级数求和感觉有问题。
有一个比较麻烦的方法，现在只有a，b和E三个变量，假定a=常数10，现在来求E和b的关系
用for循环b的值，令b=0:0.001:0.05，每次循环前N=a/2*b+1可求出N值，N不是整数，这里应该要圆整为整数（？），然后每次循环可以得出E的一个值。
这样最后可以得到E和b的关系图。
你自己再看看吧，或者再等等高手吧
那能不能把写好的命令帮我贴一下，或者传给我啊QQ328568723。因为我觉得自己现在一点思路没有，写不出了。非常感谢。那E和a的图怎么得出呢？另外N确实要取成整数
追答function qiue(e,N1)
K1=(1-e)/(1+e);
K2=2/(1+e);
K3=(2*e)/(1+e);
E1=zeros(201,40);
a=10;
for z1=1:201
for b1=1:40
fg1=0;fg2=0;fg3=0;fg4=0;N=fix(a/(0.02*b1)+1);
for n=0:N
fg1=fg1+(K1^(2*n)/(0.1*(z1-1)-2*n*0.01*b1)^2);
end
for n=1:N
fg2=fg2+(K1^(2*n-1)/(0.1*(z1-1)-2*a+2*(n-1)*0.01*b1)^2);
end
for n=N:N1
fg3=fg3+(K3*K1^(2*n)/(0.1*(z1-1)-2*n*0.01*b1)^2);
end
for n=1:N1
fg4=fg4+(K2*K1^(2*n-1)/(0.1*(z1-1)-2*a-2*n*0.01*b1)^2);
end
E1(z1,b1)=0.0000001/(4*pi*1)*(fg1+fg2+fg3-fg4);
end
end
[Z1,B1]=meshgrid(0.01:0.01:0.4,0:0.1:20);
surf(Z1,B1,E1)
end
你同学编的程序基本没有问题，现在我觉得是物理背景可能没说清楚或公式有问题
现在我有几点不明白
首先z是什么，为什么要做Z B E这三者之间的关系图，不是只做B E之间的吗
其次z ，a，b这些等等参数的值的范围和单位是不是统一了
还有你老师做的图的取值是什么，及做的是哪两个量间的关系图
能否再发点更详细的资料