Advanced Mathematics and Mechanics Applications Using MATLAB phần 8 pptx

... 2.0202e-002 7. 084 8e+000 | 4 3.0303e-002 6.6 688 e+000 | 5 4.0404e-002 6.2776e+000 Material deleted for publication | 296 2.9596e+000 6.2776e+000 | 297 2.9697e+000 6.6 688 e+000 | 2 98 2.9798e+000 7. 084 8e+000 | ... indx=[indx;2*ones(length(frqso),1)]; 83 : [frqs,k]=sort(frqs); modes=modes(:,:,k); 84 : indx=indx(k); cptim=cptim+cpto; 85 : end 86 : 87 : if nargin==6, return, end 88 :...

Ngày tải lên: 08/08/2014, 11:21

... close; 78: 79: % Plot the motion 80 : for j=1:nm 81 : % Compute and scale the applied force 82 : f=f1*cos(w*t(j))+f2*sin(w*t(j)); 83 : f=.5*f; fa=abs(f); sf=sign(f); 84 : xj=x(j,:); xmaxj=max(xj); 85 : ... Edition Advanced Mathematics and Mechanics Applications Using Howard B. Wilson University of Alabama Louis H. Turcotte Rose-Hulman Institute of Technology David Halper...

Ngày tải lên: 08/08/2014, 11:21

... tim=tim/ngquad+tbpwf-tclock; pe=100*(v./e-1); 84 : G=[v,100*ones(5,1),pe,tim/5*ones(5,1)]; 85 : 86 : format short e 87 : disp(’ ’) 88 : disp(’ Results Using Function quadl’) 89 : disp( 90: ’ Integral Function ... integration using a 76: % vector integrand. Save results in array G 77: 78: ngquad=0; tim=0; t0=clock; 79: while tim<secs 80 : v=ftest(bp,6)*wf; 81 : ngquad=ngquad+...

Ngày tải lên: 08/08/2014, 11:21

... axis’) 78: ylabel(’y axis’), zlabel(’z axis’) 79: colormap([1 1 1]); grid on, figure(gcf) 80 : end 81 : 82 : %============================================== 83 : 84 : function mat=rotatran(v) 85 : % 86 : ... ratcof. 144: % © 2003 by CRC Press LLC 86 : % 87 : % m,c,k - mass, damping, and stiffness matrices 88 : % f1,f2 - amplitude vectors for the sine and cosine 89 : % forcing...

Ngày tải lên: 08/08/2014, 11:21

... c(n/2) and 181 : % c(-n/2). For a particular value of n, the 182 : % highest usable harmonic is n/2-1. 183 : % 184 : % User m functions called: none 185 : % 186 : 187 : x=x(:); n=length(c); 188 : if ... GEOMETRIC PROPERTY APPLICATIONS 195 182 : v=quadl(@fshift,0,pi,tol,trace,varargin{:}); 183 : end 184 : 185 : %========================================= 186 : 187 : function u...

Ngày tải lên: 08/08/2014, 11:21

... and stiffness matrices 282 : % 283 : % User m functions called: none. 284 : % 285 : 286 : m=diag(masses); 287 : b=flipud(cumsum(flipud(masses(:))))* 288 : gravty./lngths; 289 : n=length(masses); k=zeros(n,n); ... Integrator of Order 6 Figure 8. 2: Stability Zone for Explicit Integrator of Order 6 © 2003 by CRC Press LLC 284 ADVANCED MATH AND MECHANICS APPLICATIONS USING...

Ngày tải lên: 08/08/2014, 11:21

... 84 .2714 286 : 14.7960 34.1543 53.3 483 46.2 980 66.1476 85 . 682 5 287 : 16.2234 35.4999 54.6719 47 .81 78 67.5905 87 . 087 0 288 : 17.6159 36 .83 36 55. 988 5 49.3204 69.0240 88 . 484 6 289 : 18. 980 1 38. 1563 57.2 984 50 .80 72 ... 28. 6266 47.9743 67.1 982 60.2703 79.9960 283 : 33. 988 7 30.0337 49.33 08 68. 5339 61.7567 81 .4291 284 : 11.7916 31.42 28 50.6 782 43.19 98 63...

Ngày tải lên: 08/08/2014, 11:21

... off; figure(gcf); 78: print -deps apprxdat 79: 80 : %============================================= 81 : 82 : function z=vibfun(x) 83 : % 84 : % z=vibfun(x) 85 : % ~~~~~~~~~~~ 86 : % 87 : % This function ... u=reshape(u,n,m); 84 : if nargin==1 | noplot==0 85 : z=reshape(z,n,m); 86 : surf(real(z),imag(z),u); xlabel(’x axis’); 87 : ylabel(’y axis’); 88 : title([’Mixed Boundary Valu...

Ngày tải lên: 08/08/2014, 11:21

... CRC Press LLC 84 : % end slope and makes the discontinuity 85 : % of y’’’ at xd(end-1) small. 86 : % endc=[4,right_slope] imposes the right 87 : % end slope and makes the discontinuity 88 : % of y’’’ ... LLC 6 08 ADVANCED MATH AND MECHANICS APPLICATIONS USING MATLAB Routine Chapter Description splinerr 2 Program showing differential geome- try properties of a space curve...

Ngày tải lên: 08/08/2014, 11:21

... 143: sprintf(’R = %8. 3f’,R), 144: sprintf(’I = %8. 3f’,I), 145: sprintf(’A = %8. 3f’,A), 146: sprintf(’q0 = %8. 3f’,q0), 147: sprintf(’q1 = %8. 3f’,q1), 1 48: sprintf(’q2 = %8. 3f’,q2)); 149: w=axis; ... R2=[real(u2);imag(u2);z2]; 177: 1 78: % Get curve properties from crvprp3d 179: [T,N,B,kappa]=crvprp3d(R1,R2); 180 : 181 : %============================================== 182 : 18...

Ngày tải lên: 08/04/2014, 09:57