<div style="text-align:left;"> </div> <div>소인, 대감께 감히 여쭐것이 있어 이렇게 오게 되었습니다....</div> <div>부디 하루종일 붙잡고 머리 싸매는 까막눈 소인에게 힌트를 좀 주시면....</div> <div> </div> <div>****사진에 붙은 문제 번호는 무시해주세요...</div> <div> </div> <div> </div> <div><img width="475" height="167" style="border:;" alt="캡처.PNG" src="http://thimg.todayhumor.co.kr/upfile/201806/15279511322899a615a9f04fc5b917f13c1324ac6b__mn372289__w475__h167__f14897__Ym201806.png" filesize="14897"></div> <div> </div> <div>라는 문제가 있어서</div> <div> </div> <div>강의 교재에 있는 아래 예제 두개를 보고</div> <div><font size="3"> </font></div> <div><font size="3">참고한 예제 1.</font></div> <div align="left" class="MsoNormal" style="margin:0cm 0cm 0pt;text-align:left;"><font size="2"><font face="바탕"><span lang="en-us">-fsolve </span>함수를 이용하여 다음의 문제를 풀어 보자<span lang="en-us">.</span></font></font><p></p></div> <div><font size="3"> </font></div> <div align="left" class="MsoNormal" style="margin:0cm 0cm 0pt;text-align:left;"><span style="color:#000000;"><font size="2"><font face="바탕">먼저 다음과 같은<span lang="en-us"> m-file(</span>이름<span lang="en-us">: test.m)</span>을 만들자<span lang="en-us">.</span></font></font></span><p></p></div> <div><font size="3"></font></div> <div align="left" class="MsoNormal" style="margin:0cm 0cm 0pt;text-align:left;"><span lang="en-us" style="color:#000000;"><font size="2"><font face="바탕">function [f]=test(x)</font></font></span><p></p></div> <div><font size="3"> </font></div> <div align="left" class="MsoNormal" style="margin:0cm 0cm 0pt;text-align:left;"><span lang="en-us" style="color:#000000;"><font size="2"><font face="바탕">f(1,1)=x(1)^2+x(2)^2-2;</font></font></span><p></p></div> <div><font size="3"> </font></div> <div align="left" class="MsoNormal" style="margin:0cm 0cm 0pt;text-align:left;"><span lang="en-us" style="color:#000000;"><font size="2"><font face="바탕">f(2,1)=x(1)*x(2)+x(1)+2*x(2)-4;</font></font></span><p></p></div> <div><font size="3"> </font></div> <div align="left" class="MsoNormal" style="margin:0cm 0cm 0pt;text-align:left;"><span lang="en-us" style="color:#000000;"><font size="2"><font face="바탕">return</font></font></span><p></p></div> <div><font size="3"> </font></div> <div align="left" class="MsoNormal" style="margin:0cm 0cm 0pt;text-align:left;"><span lang="en-us" style="color:#000000;"></span><p><font face="바탕" size="2"></font></p><font face="바탕" size="2"></font></div> <div><font face="바탕" size="2"> </font></div><font face="바탕" size="2"><span style="color:#000000;"><font size="2"><font face="바탕">다음으로 아래를 실행하면 근을 찾아 준다<span lang="en-us">.</span></font></font></span></font><p></p> <div></div> <div><font size="3"> </font></div> <div align="left" class="MsoNormal" style="margin:0cm 0cm 0pt;text-align:left;"><span lang="en-us"><font size="2"><font face="바탕">>> [x fval] = fsolve(@test,[2 0])</font></font></span><p></p></div> <div><font size="3"> </font></div> <div align="left" class="MsoNormal" style="margin:0cm 0cm 0pt;text-align:left;"><span lang="en-us"><font size="2"><font face="바탕">Optimization terminated: first-order optimality is less than options.TolFun.</font></font></span><p></p></div> <div><font size="3"> </font></div> <div align="left" class="MsoNormal" style="margin:0cm 0cm 0pt;text-align:left;"><span lang="en-us"><font size="2"><font face="바탕">x =</font></font></span></div> <div><font size="2"><font face="바탕"> </font></font></div><font size="2"><font face="바탕"><span lang="en-us"><font size="2"><font face="바탕"><span>    </span>1.0000<span>    </span>1.0000 </font></font></span></font></font><div> </div><span lang="en-us"><font size="2"><font face="바탕">fval =</font></font></span> <div></div> <div align="left" class="MsoNormal" style="margin:0cm 0cm 0pt;text-align:left;"><span lang="en-us"><font size="2"><font face="바탕"><span lang="en-us"><font size="2"><font face="바탕"><span>  </span>1.0e-007 *</font></font></span></font></font></span></div> <div align="left" class="MsoNormal" style="margin:0cm 0cm 0pt;text-align:left;"><span lang="en-us"><font size="2"><font face="바탕"><span lang="en-us"><font size="2"><font face="바탕"><span>    </span>0.1893</font></font></span></font></font></span></div> <div align="left" class="MsoNormal" style="margin:0cm 0cm 0pt;text-align:left;"><span lang="en-us"><font size="2"><font face="바탕"><span lang="en-us"><font size="2"><font face="바탕"><span lang="en-us"><font size="2"><font face="바탕"><span>   </span>-0.0908</font></font></span></font></font></span></font></font></span><p></p></div> <div><font size="3"> </font></div> <div><font size="3"></font> </div> <div>참고한 예제 2</div> <div style="text-align:left;"><img width="742" height="263" style="border:;" alt="캡처.PNG" src="http://thimg.todayhumor.co.kr/upfile/201806/15279513096c46834c49c64a5996612a5a315e0061__mn372289__w742__h263__f28888__Ym201806.png" filesize="28888"></div> <div> </div> <div>function [F]=myfunc(x)<br>    Ca=x(1);<br>    Cb=x(2);<br>    Cc=x(3);<br>    Cd=x(4);<br>    <br>    k1=1;<br>    k2=0.2;<br>    k3=0.05;<br>    k4=0.4;<br>    Vr=100;<br>    Q=50;<br>    Ca0=1;<br>    <br>    F(1)=-Ca+Ca0+Vr*(-k1*Ca-k2*Ca^1.5+k3*Cc^2)/Q;<br>    F(2)=-Cb+Vr*(2*k1*Ca-k4*Cb^2)/Q;<br>    F(3)=-Cc+Vr*(k2*Ca^1.5-k3*Cc^2+k4*Cb^2)/Q;<br>    F(4)=-Cd+Vr*(k4*Cb^2)/Q;<br>return</div> <div>다음으로 fsolve함수를 이용하여 근을 찾는다.</div> <div>>> [x fval]=fsolve(@myfunc,[1.0 0.0 0.0 0.0])<br>Optimization terminated: first-order optimality is less than options.TolFun.</div> <div>x =<br>    0.3189    0.7839    0.5350    0.4916<br>fval =<br>  1.0e-009 *<br>    0.0696   -0.7838    0.7144    0.7841</div> <div> </div> <div> </div> <div>-----------------------------------------------------------------------------------------------------------------------------</div> <div> </div> <div><img width="512" height="174" style="border:;" alt="캡처.PNG" src="http://thimg.todayhumor.co.kr/upfile/201806/15279514218ed77819344449dcbd6d791655a88252__mn372289__w512__h174__f15013__Ym201806.png" filesize="15013"></div> <div>이 문제를 풀기 위해서 위 두개를 보고 </div> <div>1번 문제를 보고 fsolve를 사용하면 되겠다는 생각이 들어서 </div> <div>fsolve를 쓰는 다른 예제인 2번 예제에서 미지수를</div> <div> Ca=x(1);<br>    Cb=x(2);<br>    Cc=x(3);<br>    Cd=x(4);</div> <div>이런식으로 뒀으니까 저도 식에 있는 미지수를</div> <div>x=x(1);<br>y=x(2);<br>z=x(3);</div> <div style="text-align:left;"> </div> <div>이렇게 둬서</div> <div>솔직히 저거나 이거나 다를거 업잖아요</div> <div>그래서 결국에</div> <div> </div> <div>이렇게 짜 봤어요</div> <div>function [F]=mat04(x)<br>x=x(1);<br>y=x(2);<br>z=x(3);</div> <div>F(1,1)=2*x-y-z;<br>F(2,1)=X^4-0.5*y^2-0.5*z^2;<br>F(3,1)=2*exp(x)-exp(y)*(1+z);</div> <div>return</div> <div> </div> <div>그다음 예제에 나온거랑 비슷하게</div> <div><br>[x fval]=fsolve(@mat04,[1.0 0.9 1.2])</div> <div>이렇게 돌렸어요</div> <div>초기값이 x=1</div> <div>y=0.9</div> <div>z=1.2 </div> <div>로 두라고 해서요</div> <div> </div> <div>2번이랑 문제 형식도 비슷하다고 생각하고 결과값 돌리는거도</div> <div>[x fval]=fsolve(@myfunc,[1.0 0.0 0.0 0.0])<br>[x fval]=fsolve(@mat04,[1.0 0.9 1.2])</div> <div>다를거 하나도 없는데 </div> <div>이런 메세지가 떠서요</div> <div> </div> <div> [x fval]=fsolve(@mat04,[1.0 0.9 1.2])<br>인덱스가 배열 경계를 초과합니다.</div> <div>오류 발생: mat04 (line 3)<br>y=x(2);</div> <div>오류 발생: fsolve (line 242)<br>            fuser = feval(funfcn{3},x,varargin{:});</div> <div>원인:<br>    Failure in initial objective function evaluation. FSOLVE cannot continue.<br> <br>근데 제가 1번의 x(1)*x(2) 이렇게 되어있던걸</div> <div>x=x(1)</div> <div>y=x(2) 이런식으로 고쳐서 돌렸을때도 위와 비슷한 메세지가 떴었어요 </div> <div>대체 뭐가 문제인건가요....</div> <div>아 대체왜안대 ㅈㅇ매ㅓ야ㅜㅁㅈ우ㅕㅁㅇㄿㅁ농ㅎ'메대럼ㄷ'랴ㅓㄴ애푸내ㅕㅎㅁㄷ'ㅑㅗㅁ'게ㅑ런해ㅜㄹㄴㅍ;ㅐ눙ㅎ['ㅁㅔㅑ]댷ㄷ메해ㅏㅁㄷ]ㅎ</div> <div>ㄴ마ㅙㄴㅇ이거말고도할거만은데ㅔㅇㅁ-아므라췅ㅍ'ㅔ냗ㅎㅁ뗴ㅑ헌'야ㅔ퍼ㅑㅑㅁ렂ㅁ</div> <div> </div> <div> </div> <div>소인을좀 구제해 주시옵소서 대감..</div>