> I can't read your f1.
==========
sorry,f1 is:
f1(z,theta) = I(z)*sin(theta)*exp(jk*z*cos(theta))
==========
> Then you say that z is a vector
> z= linspace(0,0.7,100);
> then later say that
> I(z) = is a given vector same length with
> z => length(I(z)) = length(z)
> is I(z) a function of z?
> If so, what is the definition of this function I(z) ?
==========
yes,I(z) was a function from z (but no explicity).you suppose I have a vector from z and
I(z) that are same size.
in fact I note that I(zi) & k(zi) are correspond with zi , where i=1,2,3,...,length(z);
I have to plot f2 versus theta in polar diagram.
in fact my f2 is as following (I decompose it):
f2(theta) = integral[f1(z1,theta),dz,0,z1]+integral[f1(z2,theta),dz,z1,z2]+...
integral[f1(z3,theta),dz,z2,z3]+...;
where:
f1(z1,theta)=I(z1)*sin(theta)*exp(jk(z1)*z1*cos(theta))
f1(z2,theta)= I(z2)*sin(theta)*exp(jk(z2)*z2*cos(theta))
...
where:
zi = i*delta_z = i*[0.7/length(z)] i=1,2,3,...,lenght(z)
I(zi) = ith element in vector of value I(z) ; i=1,2,3,...,lenght(z)
k(zi) = ith element in vector of value k ; i=1,2,3,...,lenght(z)
I have all of them.
now,I want to calculate and next plot f2(theta) versus theta.
how integrate from f2(theta) relation and next with theta=linspace(0,2*pi,360) plot it?
for example you suppose we have:
m=100;
z=linspace(0,0.7,m);
real_k = linspace(10,100,m);
imag_k = linspace(1,5,m);
k = real(k)+1j*imag(k);
I = linspace(12,120,m); % I= I(z)
theta = linspace(0,2*pi,360);
NOW,how integrate from f2 relation?
