Labfans是一个针对大学生、工程师和科研工作者的技术社区。 | 论坛首页 | 联系我们(Contact Us) |
![]() |
![]() |
#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'
__________________
坚持就是胜利,努力就有奇迹。 |
![]() |
![]() |