程序中一句代码（第6行）的含义是什么？

苏小贝

程序代码：

1. clear;clc;close all;
   
2. gear_par = gear_parameters
   
3. speed = 120;     % rpm
   
4. % R_par = stoch_par(speed);
   
5. wn    = gear_par.wn;
   
6. TT    = .01*2*pi/(speed*2*pi*1/60/wn);  % 20 shaft period
   
7. NN    = 1e5;
   
8. tt    = linspace(0,TT,NN);
   
9. dt    = tt(2)-tt(1);
   
10. randn('state',100);
    
11. dW    = sqrt(dt)*randn(1,NN);
    
12. y0    = [0,0,0.2,0.1]';
    
13. Xem   = zeros(4,NN);
    
14. Xtemp = y0;
    
15. for h = 1:NN
    
16.     Winc   = dW(h);
    
17.     th     = tt(h);
    
18.     [Ft,Gt] = drift_coeff(th,Xtemp,speed,R_par);
    
19.     Xtemp  = Xtemp + dt*Xtemp.*Ft + Gt*Winc;
    
20.     Xem(:,h) = Xtemp;
    
21.     h
    
22. end
    
23. plot(tt,Xem(3,:))   
    

复制代码

第6行代码如何理解？

犟牛

wn是什么参数?

苏小贝

犟牛 发表于 2015-3-17 22:14
wn是什么参数?

系统的固有频率

犟牛

speed*2*pi*1/60 是轴的转动角频率
speed*2*pi*1/60/wn 是角频率的无量纲化
2*pi/(speed*2*pi*1/60/wn) 是无量纲化的周期

jkjxl

神马东东？？？？