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

... 20 03 by CRC Press LLC 131 : 132 : case 2 133 : % Slope given at both ends 134 : [b,c]=unmkpp(spline(xd,[endc(2);yd;endc (3) ])); 135 : 136 : case 3 137 : % Slope at left end given. Compute right end 138 : ... ellipse with 30 : % semi-diameters of 2 and 1. 31 : if nargin==0 32 : m=21; th=linspace(0,2*pi,m); 33 : x=2*cos(th); y=sin(th); x(m)=2; y(m)=0; 34 : end 35 : 36 : % U...

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

... F2*SIN(W*T)’) 29: disp( 30 : ’ WITH ANIMATION OF THE RESPONSE’) 31 : disp(’ ’) 32 : 33 : % Example data used when nargin > 0 34 : if nargin > 0 35 : m=1; c= .3; k=1; f1=1; f2=0; w=2; x0=0; v0=2; 36 : tmax=25; ... relat- © 20 03 by CRC Press LLC 3. 5.3TheStructuralDynamicsEquation 3. 6ComputingNaturalFrequenciesforaRectangularMembrane 3. 7ColumnSpace,NullSpace,OrthonormalBases...

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

... and location 30 : IntSupX= [1; 2]; IntSupDelta=[0; 0]; 31 : % End (left and right) conditions 32 : EndCondVal= [0; 0; 0; 0]; % magnitude 33 : % 1=shear,2=moment ,3= slope,4=delta 34 : EndCondFunc= [3; ... interior supports 31 2: ns=size(IntSup,1); 31 3: if ns > 0 31 4: ysprt=IntSup(:,1); r=IntSup(:,2); 31 5: snf=sngf(x,r,1); 31 6: else 31 7: ysprt=[]; r=[]; snf=zeros(NoSegs,...

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

... df and nf 30 : rb=df(1); rt=df(2); h=df (3) ; 31 : [xf,yf,zf]=frus(rb,rt,h,nf); zf=zf 35 *h; 32 : 33 : % Plot four figures combined together 34 : hold off; clf; close; 35 : surfmany(x1,y1,z1,x2,y2,z2,x3,y3,z3,xf,yf,zf) 36 : ... equation  x a  2 +  y b  2 +  z c  2 =1 © 20 03 by CRC Press LLC 131 : x=varargin{k}; y=varargin{k+1}; 132 : z=varargin{k+2}; 133 : x=x(:); y=y(:);...

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

... 131 : 000]’; 132 : yft=[ 133 : 0 0.92 -0.25 1.00 -0.29 134 : 0.46 -0.16; -0.97 -0.49 -0. 83 135 : 0.95 0.86 -0.76 0.85; -0.55 136 : 0 .36 -0.52 -0 .38 0.02 -0.19 137 : 0.08; -0.26 0.24 0.00 0 138 : ... -b2 38 : % c = [-2,2] 39 : % 40: % Then the exact value is 128*pi/5. © 20 03 by CRC Press LLC 184 ADVANCED MATH AND MECHANICS APPLICATIONS USING MATLAB −2 −1 0 1...

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

... Implicit ’, 33 4: ’4th Order Integrator’]); 33 5: xlabel(’time’); 33 6: ylabel(’solution error measure’); 33 7: lg1=[’h= ’, num2str(h4), 33 8: ’, relative cputime= ’, num2str(tr4)]; 33 9: lg2=[’h= ... cputime= ’, num2str(Tr2)]; 32 8: legend(lg1,lg2,2); figure(gcf); 32 9: disp(’Press [Enter] to continue’); pause 33 0: % print -deps deislne2 33 1: 33 2: plot(t4,er4,’-’,T4,Er4,’ ’);...

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

... 79.9960 2 83: 33 .9887 30 . 033 7 49 .33 08 68. 533 9 61.7567 81.4291 284: 11.7916 31 .4228 50.6782 43. 1998 63. 231 3 82.8 535 285: 13. 3 237 32 .7958 52.0172 44.75 93 64.6947 84.2714 286: 14.7960 34 .15 43 53. 34 83 46.2980 ... 38 .68 43 32.1897 52.0077 71.4 639 63. 6114 272: 20.8070 39 .9526 33 .7166 53. 435 2 72.8506 65.15 93 2 73: 22.0470 41.2 135 35 .2187 54.8517 74. 230...

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

... numerical 130 : % integration 131 : % 132 : % F - The value of the Cauchy integral 133 : % corresponding to matrix argument z 134 : % 135 : % User m functions called: gcquad splined 136 : % 137 : 138 : n=length(fb); ... abs(C(2:ntrms+1)).^2); 31 : ster=100*abs(stfaprx-stfexct)/stfexct; 32 : 33 : % Plot circumferential and normal stresses at 34 : % the boundary 35 : th=180/pi*...

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

... mat=[null(v’),v]; 33 2: if det(mat)<0, mat(:,1)=-mat(:,1); end 33 3: 33 4: %=========================================== 33 5: 33 6: function [xmin,fmin,m,ntype]=nelmed( 33 7: F,x0,dx,epsx,epsf,M,ifpr,varargin) 33 8: ... R=[X(j);Y(j);Z(j)]; 535 : 536 : %================================================= 537 : 538 : function R=rads 539 : % R=rads 540: % Radii for the problem solut...

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

... and UsingSymbolicMath 5.6NumericalandSymbolicResultsfortheExample 5.7GeometricalPropertiesofaPolyhedron 5.8EvaluatingIntegralsHavingSquareRootTypeSingularities 5.8.1ProgramListing 5.9GaussIntegrationofaMultipleIntegral 5.9.1Example:EvaluatingaMultipleIntegral 6FourierSeriesandtheFastFourierTransform 6.1DeÞnitionsandComputationofFourierCoefÞcients 6.1.1TrigonometricInterpolationandtheFastFou...

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