Labfans是一个针对大学生、工程师和科研工作者的技术社区。 | 论坛首页 | 联系我们(Contact Us) |
![]() |
![]() |
#1 |
初级会员
注册日期: 2010-05-21
年龄: 36
帖子: 2
声望力: 0 ![]() |
![]()
function jrl = Jr_l(eps)
if length(eps)>1 error('epsilon must be s scale,not a vector!'); end if eps>1 Fai_L=pi; else Fai_L = acos(1-2*eps); %积分限。 end syms f int_f = ((1-(1-cos(f))/(2*eps))^(1.1))*cos(f)/(2*pi); %被积函数 jrl =double(int(int_f,f,-Fai_L,Fai_L)); %积分 jrl=real(jrl); 上面这个程序是可运行的,输入一个变量就会出一个结果。 我想画它的曲线图 eps = 0.1:0.1:5; J = zeros(size(eps)); for i=1:length(eps) J(i) = Jr_l(eps(i)); end clear i plot(eps,J); xlabel 'eps' ylabel 'Jr' 却总是显示错误,请各位帮帮忙:'( |
![]() |
![]() |