Engineering Matlab Problem Solving phần 3 docx

... abs(z3); phi3 = angle(z3); 78 >>z1 =3+ 4j; >> z2 = z1 * exp(j*pi/2) z2 = -4.0000+ 3. 0000i >>z3=j*z1 z3 = -4.0000+ 3. 0000i >> theta1 = (180/pi) * angle(z1) theta1 = 53. 130 1 >> ... z 2 )+z 3 = z 1 +(z 2 + z 3 ) (z 1 z 2 )z 3 = z 1 (z 2 z 3 ) Distributivity: Multiplication can be distributed across a sum without changing the result. z 1 (z 2 + z 3 )...

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

... variable. Example 8.1 Avoiding division by zero Consider the following: >> x=( -3: 3) /3 x= -1.0000 -0.6667 -0 .33 33 0 0 .33 33 0.6667 1.0000 >> sin(x)./x Warning: Divide by zero. ans = 0.8415 0.9276 ... number eps, which is approximately 2.2×10 −16 , resolving the divide-by-zero problem. >> x = x + (x==0)*eps x= -1.0000 -0.6667 -0 .33 33 0.0000 0 .33 33 0.6667 1.000...

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

... Precision . . . 33 3 Files and File Management 37 3. 1 FileManagementDefinitionsandCommands 37 3. 2 SavingandRestoringMatlabInformation 39 3. 3 ScriptM-Files 43 3.4 ErrorsandDebugging 47 3. 5 MatlabSearchPath,PathManagement,andStartup ... basic Matlab command is the mathematical expression, which has the following prop- erties: 15 Contents 1 Engineering Problem Solving 1 1.1...

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

... 5.0 833 333 333 333 34e+001 14 decimal digits plus exponent format short g 50. 833 3 better of format short or format short e (default), switching for ans > 1000 format long g 5.0 833 333 333 333 34e+001 ... ratio = 610./12 ratio = 50. 833 3 Matlab Command ratio Comments format short 50. 833 3 4 decimal digits format long 50. 833 333 333 333 34 14 decimal digits format short e 5.0 833 e+...

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

... is used to separate rows. >>f=[1 23; 456] f= 1 23 456 >>g=f’ g= 14 25 36 >>h=[1 23 456 789] h= 1 23 456 789 88 >> F = 3. 0.^A F= 9 2 43 729 Array-array mathematics on other than ... 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 −2 0 2 4 6 8 10 12 Multiple Axis Plot x 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 −10 0 10 20 30 40 50 60 Figure 5 .3: Plot with two polynomial curves, .....

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

... conv(a,b) c= 4 -14 11 -35 58 -9 26 -32 114 −1 −0.5 0 0.5 1 1.5 2 2.5 3 3.5 4 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 x y Figure 6.14: Plot of the function y =0.025x 5 − 0.0625x 4 − 0 .33 3x 3 + x 2 The result ... the polynomial. To specify C(s)=3A(s) for A(s) above, the following commands are used: >>a=[1 -30 -12]; >>c =3* a c= 3- 9 0 -3 6 Thus C(s)=3s 4 − 9s 3 − 3s +6 • Multipli...

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

... =  13 −15  the determinant is |A| =1·5 3 · (−1)=8 For a 3 3 matrix A, the determinant is |A| = a 11 (a 22 a 33 − a 23 a 32 ) − a 12 (a 21 a 33 − a 23 a 31 )+a 13 (a 21 a 32 − a 22 a 31 ) = ... a 22 a 31 ) = a 11 a 22 a 33 − a 11 a 23 a 32 − a 12 a 21 a 33 + a 12 a 23 a 31 + a 13 a 21 a 32 − a 13 a 22 a 31 Thus, for B =    130 −152 121   ...

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

... form:     t 5 f t 4 f t 3 f 5t 4 f 4t 3 f 3t 2 f 20t 3 f 12t 2 f 6t f        a 1 a 2 a 3    =    θ 1 (t f ) − θ 1 (0) 0 0    2 03 b = T\d % % Equations of motion L1=4; L2 =3; t = linspace(0,2,401); tq ... minimum splmin = -8.26 83 kmin = 32 >> fmin = fi(kmin) % frequency at minimum fmin = 3. 2828e+0 03 Thus, by this analysis, the human ear is most sens...

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

... = 400/R^2*(pi*R ^3- 750)+400*pi*R >> Rmins = solve(dCdR,R) % Solve dC/dR for R: Rmin Rmins = [ 5/pi *3^ (1 /3) *(pi^2)^(1 /3) ] [ -5/2/pi *3^ (1 /3) *(pi^2)^(1 /3) +5/2*i *3^ (5/6)/pi*(pi^2)^(1 /3) ] [ -5/2/pi *3^ (1 /3) *(pi^2)^(1 /3) -5/2*i *3^ (5/6)/pi*(pi^2)^(1 /3) ] >> ... s >> A = s^4 -3* s ^3 -s +2; >> C = 3* A C= 3* s^4-9*s ^3- 3*s+6 • Multiplication: >&g...

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

... = 1 2 3 »b = [1;2 ;3] b = 1 2 3 »[ 1 3 5 2 4 6 3 6 9] ans = 1 3 5 2 4 6 3 6 9 Exercise: Input the following matrices:: 1 23 538 2 .3 5.6 10 , 1 23 100 105 There are several useful shortcuts ... 57 17 42 -61 -65 Answer: 57X 4 + 17X 3 + 42X 2 - 61X - 65 Exercise: What is (3X 3 + 2X + 5) * (X 3 + 2X 2 - 2) ? answer: 3X 6 + 6X 5 + 2X 4 + 3X 3 + 10X 2 - 4X...

Ngày tải lên: 06/08/2014, 13:22