From: parul agrawal on 19 Mar 2010 02:53 syms z Ps = max_signal_power*(10^(-total_loss_s/(10*z))); Pp = max_pumping_power*(10^(-total_loss_p/(10*z))); Ps_desh = Ps/Ps_saturation; Pp_desh = Pp/Pp_saturation; A = Nt*(Pp_desh + Ps_desh); B = (1+(2*Ps_desh)+Pp_desh); N2 = A/B; N1 = Nt-N2; L = input('enter the length in meter'); L_desh = [0:1:L]; l = length(L_desh); F1 = int(N2,z,0,0); F2 = int(N1,z,0,0); F = ((sigma_signal_emission*F1)-(sigma_signal_absorption*F2)); G = gama_s*exp(F); for i =L_desh(1,2):L_desh(1,l); F1 = int(N2,z,0,i); F1 = double(F1); F2 = int(N1,z,0,i); F2 = double(F2); F = ((sigma_signal_emission*F1)-(sigma_signal_absorption*F2)); F = double(F); G = [G,gama_s*exp(F)]; G = double(G); end plot(L_desh,G); grid; its giving first taht explicit integral is not found ...when i tried to debug it....i found that problem is in reshape function(inside coding of int function)...when i checked the value of' ps' and tried to change it 'double' it is giving : >>double(Ps) ??? Error using ==> reshape To RESHAPE the number of elements must not change. help plz
|
Pages: 1 Prev: using case statement on string inputs when number of case choices Next: simulink question |