请教吴老师：嵌套函数输出结果的问题

chunshui2003

吴老师，你好，又一次麻烦你，依然是之前的嵌套函数求解问题，代码如下：

1. 
   
2. 
   
3. function [t,u] = chunshui2003
   
4. 
   
5. %%%% 基本参数数值 %%%%%
   
6. mr = 600*10^3;            
   
7. De = 0.5*10^6;            
   
8. Ke = 0.5*10^10;           
   
9. delta_0 = 18*10^-3;      
   
10. e_0 = 2*10^-3;            
    
11. omega = 13.1;              %水轮发电机额定转速，单位rad/s
    
12. 
    
13. 
    
14. %%%% 电磁刚度参数数值 %%%%%
    
15. Rg = 6.24;      
    
16. L1 = 2.1;      
    
17. k_u = 1.102;   
    
18. mu_0 = 4*pi*10^(-7);   
    
19. f = 50;        
    
20. theta = 30.64/180*pi;
    
21. phi = acos(0.875);   
    
22. p = 40;               
    
23. Fsm = 19210;         
    
24. Fjm = 24214;         
    
25. Lambda_0 = mu_0./(k_u*delta_0);      
    
26. sigma = k_u*delta_0;      
    
27. omega_f = 2*pi*f./p;      
    
28. 
    
29. 
    
30. function uu = diff_equ(t,u)
    
31. %%%%% Fx,Fy 表达式 %%%%%%
    
32. 
    
33. Fx = mr.*omega.^2.*e_0.*cos(omega.*t);
    
34. Fy = mr.*omega.^2.*e_0.*sin(omega.*t);
    
35. 
    
36. 
    
37. %%%% 电磁刚度表达式 %%%%%
    
38. 
    
39. coeff_1 = Rg.*L1.*Lambda_0./(2*sigma.^2);
    
40. coeff_2 = coeff_1.*sigma;
    
41. 
    
42. B = @(alpha) Fsm.*cos(omega_f.*t-p.*alpha);                     %积分大括号内第一项
    
43. C = @(alpha) Fjm.*cos(omega_f.*t-p.*alpha+theta+phi+pi/2);      %积分大括号内第二项
    
44. 
    
45. C11 = @(alpha) 1+cos(2.*alpha);
    
46. C12 = @(alpha) sin(2.*alpha);
    
47. C22 = @(alpha) 1-cos(2.*alpha);
    
48. C_1 = @(alpha) cos(alpha);
    
49. C_2 = @(alpha) sin(alpha);                  %不同电磁刚度项的积分内容
    
50. %
    
51. D11 = @(alpha) coeff_1.*C11(alpha).*(B(alpha)+C(alpha)).^2;
    
52. D12 = @(alpha) coeff_1.*C12(alpha).*(B(alpha)+C(alpha)).^2;
    
53. D22 = @(alpha) coeff_1.*C22(alpha).*(B(alpha)+C(alpha)).^2;
    
54. D_1 = @(alpha) coeff_2.*C_1(alpha).*(B(alpha)+C(alpha)).^2;
    
55. D_2 = @(alpha) coeff_2.*C_2(alpha).*(B(alpha)+C(alpha)).^2;                   %电磁刚度积分表达式
    
56. %
    
57. K11 = quadgk(D11,0,2*pi);
    
58. K12 = quadgk(D12,0,2*pi);
    
59. K21 = quadgk(D12,0,2*pi);
    
60. K22 = quadgk(D22,0,2*pi);
    
61. K_1 = quadgk(D_1,0,2*pi);
    
62. K_2 = quadgk(D_2,0,2*pi);              %电磁刚度求解
    
63. 
    
64. uu = zeros(4,1);
    
65. uu(1) = u(2);
    
66. uu(2) = -De/(mr*omega)*u(2)-(K11+Ke)/(mr*omega^2)*u(1)-K12/(mr*omega^2)*u(3)+Fx/(mr*delta_0*omega^2)-K_1/(mr*delta_0*omega^2);
    
67. uu(3) = u(4);
    
68. uu(4) = -De/(mr*omega)*u(4)-K21/(mr*omega^2)*u(1)-(K22+Ke)/(mr*omega^2)*u(3)+Fy/(mr*delta_0*omega^2)-K_2/(mr*delta_0*omega^2);
    
69. end
    
70. 
    
71. 
    
72. y0 = [0.001 0.001 0.001 0.001];
    
73. period = 2*pi/13.1;
    
74. [t,u] = ode15s(@diff_equ,[0:period/100:5*period],y0);
    
75. end
    

复制代码

之前我曾经问过，如何显示计算这个带有嵌套函数的运动微分方程结果。您说在chunshui2003这个函数名称前添加[t,u]即可，即：

1. 
   
2. 
   
3. function    [t,u] = chunshui2003
   
4. 
   

复制代码

当运行后在变量参数处并没有显示时间t和相关运动项u的名称，只有ans以及主窗口中的一系列数值，这些应该是时间t的数值，也就是时间子步。
请问一下老师，这里面是哪个环节出现了问题，应该如何操作才能够获得我所希望的[t,u]呢。
麻烦吴老师了！

rocwoods

你是怎么运行的？不是在指令窗口运行

1. 
   
2. [t,u] = chunshui2003

复制代码

吗

chunshui2003

谢谢您，吴老师！之前可能是我理解出现了错误，直接运行"chunshui2003"这个函数了。对嵌套函数的概念理解不是很清楚，以为在"chunshui2003"里面有

1. 
   
2. 
   
3. [t,u] = ode15s(@diff_equ,(0:period/100:5*period),y0);
   
4. 
   

复制代码

就可以直接获得结果，我的想法是错误的。还是需要通过运行函数并赋予输出值才可以！
虽然花了一些时间才弄懂，但能够得到吴老师的指导和帮助很开心。
祝吴老师工作顺利，新书畅销！