Advanced Mathematics and Mechanics Applications Using MATLAB phần 9 pps

... data 89: % vectors timdat and ydat are passed as global 90 : % variables. The function to be fitted is: 91 : % 92 : % y=a*exp(b*t)*cos(c*t+d) 93 : % 94 : % x - a vector defining a,b,c and d 95 : % 96 : ... Value Problem ’, 89: ’for a Circular Disk’]); 90 : grid; figure(gcf); %gra(.4), pause 91 : %print -deps mbvtest 92 : end 93 : 94 : %======================================...

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

... Done’) 87: 88: %============================================= 89: 90 : function zdot=topde(t,z,uz,c1,c2) 91 : % 92 : % zdot=topde(t,z,uz,c1,c2) 93 : % ~~~~~~~~~~~~~~~ 94 : % 95 : % This function defines the equation of motion 96 : % for a symmetrical ... externally defined time dependent 193 : % forcing function. This parameter 194 : % should be omitted if no forcing 195 : % functi...

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

... 24. 494 9 43.7155 38.1 598 57.6538 76 .96 99 69. 7230 276: 25.7051 44 .95 77 39. 6032 59. 04 09 78.3305 71.2205 277: 26 .90 74 21.2117 41.0308 60.4 194 52.6241 72.7065 278: 28.1024 22.7601 42.44 39 61.7 893 54.1856 ... 53.3483 46. 298 0 66.1476 85.6825 287: 16.2234 35. 499 9 54.67 19 47.8178 67. 590 5 87.0870 288: 17.61 59 36.8336 55 .98 85 49. 3204 69. 0240 88.4846 2 89: 18 .98 01 38.15...

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

... z=(r*ones(1,nt)).*(ones(nr,1)*exp(i*t)); 88: 89: % Use high point resolution for the 90 : % outer contour 91 : touter=pi/180*linspace(t1,t2,10*nt); 92 : zouter=r2*exp(i*touter); 93 : 94 : % Compute the series coefficients and 95 : % evaluate ... Edition Advanced Mathematics and Mechanics Applications Using Howard B. Wilson University of Alabama Louis H. Turcotte Rose-Hul...

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

... matrix based 89: % on the columns of v. 90 : % 91 : % v - a matrix having three rows and either 92 : % one or two columns which are used to 93 : % create an orthonormal triad [i,j,k] 94 : % returned ... of mat. The 95 : % third base vector k is defined as 96 : % v(:,1)/norm(v(:,1)). If v has two 97 : % columns then, v(:,1) and v(:,2) define 98 : % the xz plane with the direction...

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

... 1.3101 Results Using Function SRFV Volume = 59. 1056 Rg = [0 .91 016 0 .91 016 0.137 49] Inertia Tensor = 1.0e+003 * 0.5768 0.1167 -0.0 090 0.1167 0.5768 -0.0 090 -0.0 090 -0.0 090 1. 099 9 >> −5 0 5 −5 0 5 −0.5 0 0.5 1 1.5 x ... Results Using Function quadl’) 89: disp( 90 : ’ Integral Function Percent Computation’) 91 : disp( 92 : ’ values evaluations error seconds’) 93 : d...

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

... values 190 : % ys0,vs0 - Position and velocity at t=0 191 : % as - Acceleration ys’’(t), if this 192 : % quantity is required 193 : % 194 : % User m functions called: none. 195 : % 196 : 197 : if nft==0 © ... This 89: % should be an integer power of 2. 90 : % The default is 1024 91 : % 92 : % User m functions called: none. 93 : % 94 : 95 : if nargin<5, nft=1024; end; 96...

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

... publication | 296 2 .95 96e+000 6.2776e+000 | 297 2 .96 97e+000 6.6688e+000 | 298 2 .97 98e+000 7.0848e+000 | 299 2 .98 99e+000 7.5273e+000 | 300 3.0000e+000 7 .99 76e+000 Solution time was 0.55 secs. Reactions ... ’,num2str(neig)]); 94 : disp([’Highest frequency = ’, 95 : num2str(frqs(end))]), disp(’ ’) 96 : disp(’Press return to see modal plots.’) 97 : pause, plotmode(a,b,x,y,frqs,modes...

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

... end slope and makes the discontinuity 88: % of y’’’ at xd(2) small. 89: % 90 : if nargin<3, endc=1; end; 91 : type=endc(1); xd=xd(:); yd=yd(:); 92 : 93 : switch type 94 : 95 : case 1 96 : % y’’’(x) ... R=[X(j);Y(j);Z(j)]; 535: 536: %================================================= 537: 538: function R=rads 5 39: % R=rads 540: % Radii for the problem solutions 541: 542: R=[ 543:...

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

... simple 193 : % polynomial interpolation is employed 194 : % 195 : % xd,yd - data vectors determining the spline 196 : % curve produced by function spline 197 : % x - vector of values where the first or 198 : ... 2 Elementary Aspects of MATLAB Graphics 2.1 Introduction MATLAB s capabilities for plotting curves and surfaces are versatile and easy to understand. In fact, the effort...

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