Info
이 질문은 마감되었습니다. 편집하거나 답변을 올리려면 질문을 다시 여십시오.
quad, Warning: Infinite or Not-a-Number function value encountered.
조회 수: 1 (최근 30일)
이전 댓글 표시
Good evening, Sorry, my previous message contains some french word
I would like to plot (t,f(t)), with
f(t)=integral_0^ (2*pi*1.0e+12)*75.2933 y(t,w) dw
For that, i used 'quad' to compute my integral, but there is a problem when matlab run: Warning: Infinite or Not-a-Number function value encountered.
Thanks i advantage for you help.
function y=toltallA(w,T) It not necessary to see between begin parameter and end parameter.
%%Begin parameter L=2.9e-6; v= 21.6e+3; a=0.142e-9; N_0=6.022e+23; A= 3*sqrt(3)*a^2*N_0/4;% w_max=(2*pi*1.0e+12)*75.2933; delta_M=3; B_N=3.85e-25; %% B_U=7.7e-25;%% alpha=3; %% theta=1000; c_d=1.0e-6; h=6.626e-34; hbar= h/(2*pi); kB=1.38*1.0e-23; coef=hbar^2/(4*pi*kB*v^2); %end of my parameter %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
tau=1./(v/L+ (A*delta_M^2*2*pi*c_d)/(2*pi*v^2*w_max^2)*w.^4+(B_N+B _U*exp(-theta/(alpha*T)))*T^3*w.^2);%%
%% this my integrand %%%%%%%y(t,w)%%%%%%%%%%%%%
y= coef*v^2/T^2*tau.*w.^3.*exp(hbar*w./(kB*T))./(exp(hbar*w./(kB*T))-1).^2;
%%% MAin Program : I use quad to approximate my integral.
T=1:10:1000;
QlA=zeros(size(T));
for k=1:length(T)
QlA(k)= quad(@(Z)toltallA(Z,T(k)),1.e-20,(2*pi*1.0e+12)*75);
end
thanks
댓글 수: 1
답변 (0개)
이 질문은 마감되었습니다.
참고 항목
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!