matlab 개요와 응용 연습문제 풀이 및 m-file 제3장-1 (매트랩 독학, 솔루션)

퇴근하고 몇일만에 매트랩 문제풀기 숑숑~

어릴 적 눈높이 문제집 풀던 느낌으로 차근차근 배우고 있다. 별거 아닌데 재밌넹 ㅎㅎ 

일주일동안 11장까지 문제 다 푸는게 목푠데 아직 3장

% problem 1.

x=[-2:1:5]

y=(2*x.^2-5*x+4).^3./x.^2  %inf=무한대

x =

    -2    -1     0     1     2     3     4     5

y =

   1.0e+03 *

    2.6620    1.3310       Inf    0.0010    0.0020    0.0381    0.2560    0.9756

% problem 2.

t=0:8

y=5*sqrt(t)-(t+2).^2/0.5./(t+1)+8

t =

     0     1     2     3     4     5     6     7     8

y =

         0    4.0000    4.4044    4.1603    3.6000    2.8470    1.9617    0.9788   -0.0801

% problem 3.

g=9.81; h0=2;  n=1:8; v1=(0.85*sqrt(2*g*h0)); 

hn=0.85.^(2*n)*2  %h1=0.85^2*2 , h2=0.85^2*h1=0.85^4*2 ....

hn =

    1.4450    1.0440    0.7543    0.5450    0.3937    0.2845    0.2055    0.1485

% problem 4.

g=9.81; cd=0.5; rho=1.2; m=0.624; r=0.117; A=pi*r^2;

t=0:10

vt=sqrt(2*m*g/rho/A/cd)*(1-exp(-sqrt(rho*g*cd*A/2/m).*t))

t =

     0     1     2     3     4     5     6     7     8     9    10

vt =

         0    7.8984   12.9328   16.1417   18.1870   19.4907   20.3217   20.8513   21.1889   21.4041   21.5413

% problem 5. 

u=[14 25 -10]; 

  %a

ud=sqrt(u(1)^2+u(2)^2+u(3)^2)

  %b

sqrt(sum(u.*u))

ud =

   30.3480

ans =

   30.3480

% problem 6.

g=9.81; v0=100; t=0:2:20; theta=79;

x=v0*cosd(theta)*t;

y=v0*sind(theta).*t-1/2*g.*t.^2;

rt=sqrt(x.^2+y.^2)

rt =

         0  180.7792  323.3089  427.9925  495.4815  526.8909  524.2756  491.7770  438.6131  386.7074  381.6201

% problem 7.

u=[4 9 -5]; v=[-3; 6; -7];

a=u*v

b=dot(u,v)    % 내적 dot(u,v)  외적 cross(u,v)

b2=dot(u,v')  % dot 함수 사용시 행벡터, 열벡터 상관없음

a =

    77

b =

    77

b2 =

    77

% problem 8. 

x=2:2:10; y=3:3:15;

z=(y./x).^2+(x+y).^((y-x)./x)

z =

    4.4861    5.4123    6.1230    6.7221    7.2500

% problem 9.

h=0.7; k=8.85;

x=1:5; y=2.1:-0.1:1.7; z=2.0:0.5:4.0;

G=(h*x+k*y)./(x+y).^h+exp(h*y./z)./z.^(y./x)

G =

    9.2215    7.9378    6.9952    6.2876    5.7379

% problem 10.

format long

x=[1 0.5 0.1 0.01 0.001 0.00001 0.0000001];

y=(exp(x)-1)./x

yn1=y-1

y =

  열 1 ~ 5

   1.718281828459046   1.297442541400256   1.051709180756477   1.005016708416795   1.000500166708385

  열 6 ~ 7

   1.000005000006965   1.000000049433680

yn1 =

  열 1 ~ 5

   0.718281828459046   0.297442541400256   0.051709180756477   0.005016708416795   0.000500166708385

  열 6 ~ 7

   0.000005000006965   0.000000049433680

% problem 11.

n1=0:100;

n2=0:10000;

n3=0:1000000;

s1=(-1).^n1./(2*n1+1);

s2=(-1).^n2./(2*n2+1);

s3=(-1).^n3./(2*n3+1);

sum100=4*sum(s1)

sum10000=4*sum(s2)

sum1000000=4*sum(s3)

>> ex_part3

sum100 =

    3.1515

sum10000 =

    3.1417

sum1000000 =

    3.1416

% problem 12.

n1=0:50;

n2=0:500;

n3=0:5000;

s1=1./(2*n1+1)./(2*n1+2);

s2=1./(2*n2+1)./(2*n2+2);

s3=1./(2*n3+1)./(2*n3+2);

sum50=sum(s1)

sum500=sum(s2)

sum5000=sum(s3)

ln2=log(2)

sum50 =

    0.6883

sum500 =

    0.6926

sum5000 =

    0.6931

ln2 =

    0.6931

% problem 13.

k=[0.25 0.5 0.75];

lmax=50; tau=0.5; t=2;

l=lmax.*(1-exp(-k.*(t+tau)))

l =

   23.2369   35.6748   42.3323

% problem 14. 15 생략

 % 행렬 A+B에 대하여

  % 1) A+B=B+A  

  % 2) A+(B+C)=(A+B)+C  

  % 3) 5(A+C)=5A+5C 

  % 4) A*(B+C)=A*B+A*C

  % 5) A*B~=B*A

  % 6) A*(B*C)=(A*B)*C

  % 7) (A*B)'=B'*A'

  % 8) (A+B)'=A'+B'

16 ~ 19번은 앞에 예제 좀 풀어보고 다음 포스팅에~후후

오늘은 여기까지만 하고 치인트 보러가야쥐~

아래는 m-file인데 시간 없으면 이거 받아서 공부하셔도 좋을 것 같아요.

댓글다시면 암호 알려드릴게요 :)

ex_part3.egg