想用if语句编写裂纹综合模型的程序，总出错，

竹语随风

我编写的程序如下，不知道为什么总出错，显示theta没有定义，初学MATLAB，请不吝赐教。谢谢代码：

1. function f=fcrack(theta)
   
2. alfa=pi/6;
   
3. if -pi/2+alfa<=theta && theta<=pi/2-alfa
   
4.     f=1;
   
5. elseif pi/2-alfa<theta && theta<=pi/2+alfa
   
6.     f=1/2+1/2*cos(pi*(theta-pi/2+alfa)./(2*alfa));
   
7. elseif pi/2+alfa<theta && theta<=3*pi/2-alfa
   
8.     f=0;
   
9. elseif 3*pi/2-alfa<=theta && theta<3*pi/2+alfa
   
10.     f=1/2+1/2*cos(pi*(theta-3*pi/2-alfa)./(2*alfa));   
    
11. end

复制代码

ChaChing

求助完整格式：出错代码和出错提示!!!

竹语随风

ChaChing 发表于 2013-9-6 00:06

                               
登录/注册后可看大图

求助完整格式：出错代码和出错提示!!!

??? Input argument "alfa" is undefined.

Error in ==> fcrack at 3
if -pi/2+alfa<=theta & theta<=pi/2-alfa

这是错误提示
我想将它调用到这个函数里面

1. function dx=crackandliner(t,x,c,w)
   
2. e=0.1;
   
3. k=0.5;
   
4. beta=0;
   
5. phi0=0;
   
6. alfa=pi/6;
   
7. phi=atan(x(1)/x(3));
   
8. theta=w*t+phi0+beta-phi;
   
9. n=sin(w*t);
   
10. m=cos(w*t);
    
11. dx=zeros(4,1);
    
12. dx(1)=x(2);
    
13. dx(2)=e*w^2.*cos(w*t+phi0)-c*x(2)-x(1)+fcrack(theta)*k*m^2*x(1)...
    
14.     +fcrack(theta)*k*n*m*x(3);
    
15. dx(3)=x(4);
    
16. dx(4)=e*w^2.*sin(w*t+phi0)-c*x(4)-x(3)+fcrack(theta)*k*n^2*x(3)...
    
17.     +fcrack(theta)*k*n*m*x(1);

复制代码

主函数

1. j=1;
   
2.     c=0.15;
   
3.     for w=0.1:0.01:3
   
4.         x0=[0.01;0.01;0.01;0.01];
   
5.         options=odeset;options.RelTol=1e-4;
   
6.         [t,x]=ode45(@crackandliner,[0 400],x0,options,c,w);
   
7.         p= sqrt((x(end-400:end,1).^2)+(x(end-400:end,3).^2));
   
8.         Pmax(j,1)=max(p);
   
9.         j=j+1;
   
10.     end
    
11. w=0.1:0.01:3;
    
12. h=plot(w,Pmax);
    
13. title('p-w图');
    
14. xlabel('回转速度w');
    
15. ylabel('振幅p');

复制代码

ChaChing

1F=> 显示theta没有定义

3F=> Input argument "alfa" is undefined

是否有些冲突!?


到底LZ的版本是多少? 个人使用R2009a试跑下, 但出错提示不同呀!?

Error in ==> fcrack at 2
alfa=pi/6;

??? Output argument "f" (and maybe others) not assigned during call to
"C:\MATLABR11\work\fcrack.m>fcrack".

而这个错误乃因无返回值造成!
无返回值的原因, 检查下LZ的范围逻辑性是否合理! 或设下预设值吧

竹语随风

ChaChing 发表于 2013-9-6 15:02

                               
登录/注册后可看大图

是否有些冲突!?

不懂 啊，老师，我用07版本，和2012b都试过了。主要是theta=w*t+phi0+beta-phi，phi=atan(x(1)/x(3));我不知道fcrack能不能去出theta的范围，急啊

chybeyond

fcrack这个函数是你自己定义的吗

竹语随风

chybeyond 发表于 2013-9-9 14:25

                               
登录/注册后可看大图

fcrack这个函数是你自己定义的吗

是的，裂纹的综合模型就是这样的，我想在程序中使用这个函数

竹语随风

竹语随风 发表于 2013-9-9 15:06

                               
登录/注册后可看大图

是的，裂纹的综合模型就是这样的，我想在程序中使用这个函数

C:\Documents and Settings\Administrator\桌面

chybeyond

那你应该把fcrack贴出来，不然我们运行就出错

竹语随风

chybeyond 发表于 2013-9-9 15:07

                               
登录/注册后可看大图

那你应该把fcrack贴出来，不然我们运行就出错

这就是我定义的fcrack

1. function f=fcrack(theta)
   
2. alfa=pi/6;
   
3. if -pi/2+alfa<=theta && theta<=pi/2-alfa
   
4.     f=1;
   
5. elseif pi/2-alfa<theta && theta<=pi/2+alfa
   
6.     f=1/2+1/2*cos(pi*(theta-pi/2+alfa)./(2*alfa));
   
7. elseif pi/2+alfa<theta && theta<=3*pi/2-alfa
   
8.     f=0;
   
9. elseif 3*pi/2-alfa<=theta && theta<3*pi/2+alfa
   
10.     f=1/2+1/2*cos(pi*(theta-3*pi/2-alfa)./(2*alfa));   
    
11. end

复制代码

竹语随风

竹语随风 发表于 2013-9-9 15:11

                               
登录/注册后可看大图

这就是我定义的fcrack
function f=fcrack(theta)
alfa=pi/6;

老师，我都程序都贴在上面了，我在这里

function dx=crackandliner(t,x,c,w)
e=0.1;
k=0.5;
beta=0;
phi0=0;
alfa=pi/6;
phi=atan(x(1)/x(3));
theta=w*t+phi0+beta-phi;
n=sin(w*t);
m=cos(w*t);
dx=zeros(4,1);
dx(1)=x(2);
dx(2)=e*w^2.*cos(w*t+phi0)-c*x(2)-x(1)+fcrack(theta)*k*m^2*x(1)...
    +fcrack(theta)*k*n*m*x(3);
dx(3)=x(4);
dx(4)=e*w^2.*sin(w*t+phi0)-c*x(4)-x(3)+fcrack(theta)*k*n^2*x(3)...
    +fcrack(theta)*k*n*m*x(1);
引用了我定义的fcrack的函数

竹语随风

竹语随风 发表于 2013-9-9 15:13

                               
登录/注册后可看大图

老师，我都程序都贴在上面了，我在这里

function dx=crackandliner(t,x,c,w)

就是这个模型这个程序我弄不对，我如果用简单的模型还是可以出来图的，就是这个分段的裂纹模型不知道怎么引用

chybeyond

不好意思，忘记了，我分析了一下出错的原因，在fcrack这个函数中你列出了四种情况，但是你写的函数刚好四种情况都没满足，也就是说f没有任何值，所以提示错误。在if最后再加一个else f=。。就可以了

ChaChing

可能个人在4F没说清楚吧!?

而这个错误乃因无返回值造成!
无返回值的原因, 检查下LZ的范围逻辑性是否合理! 或设下预设值吧

基本就是13F所说的, 不过建议LZ需谨慎使用预设值(类似13F提供的)

逻辑范围在 -pi/3~pi*5/3, LZ确定theta都在这范围吗?
建议LZ了解下atan,atan2,angle等函数间的差异

竹语随风

chybeyond 发表于 2013-9-9 15:34

                               
登录/注册后可看大图

不好意思，忘记了，我分析了一下出错的原因，在fcrack这个函数中你列出了四种情况，但是你写的函数刚好四种 ...

老师，我就想引用 这个函数不好意思，最后一个取值范围应该是3*pi/2-alfa到后面的值