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' 却总是显示错误,请各位帮帮忙:'( |
![]() |
![]() |
![]() |
#2 |
高级会员
注册日期: 2008-11-07
住址: 湖南长沙
帖子: 233
声望力: 21 ![]() |
![]()
函数文件的变量有作用域,函数运行以后就没有作用了。所以解决的办法有两个,
1. 设置全局变量。但是这种办法在一般的程序里不可取,因为只要你不清除变量,它一直存在; 2. 将函数文件改写为脚本文件,这样就可以了 下面的代码是我稍微修改的,有错误,但是运行的顺序没有问题。你对着你问题再把循环改写好就可以了。 for eps = 0.1:0.1:5; J = zeros(size(eps)); if length(eps)>1 error('epsilon must be s scale,not a vector!'); if eps>1 Fai_L=pi; else Fai_L = acos(1-2*eps); %积分限。 end 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); plot(eps,J); hold on; end hold off; xlabel 'eps' ylabel 'Jr'
__________________
坚持就是胜利,努力就有奇迹。 |
![]() |
![]() |
![]() |
#3 |
初级会员
注册日期: 2008-09-28
年龄: 36
帖子: 4
声望力: 0 ![]() |
![]()
那个function文件里面的变量是不能在其他的函数段里面用的,它就仅仅是一个代号,没有实际意义的;你可以运行你那个function文件,结果出来肯定不是jrl=……,而是ans=……
你可以将画图的函数段也集成到function函数中去,直接对数组作用,画出图像。 fuction Jr_l(eps) J = zeros(size(eps)); if length(eps)>1 error('epsilon must be s scale,not a vector!'); 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); plot(eps,J); hold on; end hold off; xlabel 'eps' ylabel 'Jr' end |
![]() |
![]() |