math - Solving Differential equations in Matlab, ode45 -

-

i'm trying solve system 3 differential equations function ode45 in matlab. not understand errors getting , use understanding im doing wrong.

the differential equations following:

f1 = -k1y1+k2(y2-y1) f2 = -k2(y2-y1)+k3(y3-y2) f3 = -k3(y3-y2)

and code in matlab this:

function dz = kopplad(t, z) global m1 m2 m3 k1 k2 k3  dz = [z(2) -k1*z(1)/m1 + k2*(z(2)-z(1))/m1 z(4) -k2*(z(2)-z(1))+k3(z(3)-z(2))/m2 z(6) -k3(z(3)-z(2))/m3];  global m1 m2 m3 k1 k2 k3 m1 = 0.75; m2 = 0.40; m3 = 0.65; k1 = 0.85; k2 = 1.1; k3 = 2.7; [t, z] = ode45(@kopplad, [0, 50], [0 -1 0 1 0 0]); plot(t, z(:,1))  hold on plot(t,z(:,3),'--') plot(t,z(:,5),'*')  legend('y_1', 'y_2', 'y_3'); hold off

the errors i'm recieving following:

attempted access k3(1.00002); index must positive integer or logical.

error in kopplad (line 3) dz = [z(2)

error in ode45 (line 262) f(:,2) = feval(odefcn,t+ha(1),y+f*hb(:,1),odeargs{:});

error in diffekv (line 6) [t, z] = ode45(@kopplad, [0, 50], [0 -1 0 1 0 0]);

a while passed since did matlab programming, far remember, should pass variables function, i.e. write @(x,y)kopplad(x,y), if understand code in correct way. if rest (global variables , equations) correct, shall fine.


Comments

Post a Comment

Popular posts from this blog

java - Jmockit String final length method mocking Issue -

-
using jmockit 1.2 jar, trying mock string's length method getting unexpected invocation exception: failed: test java.lang.illegalstateexception: missing invocation mocked type @ point; please make sure such invocations appear after declaration of suitable mock field or parameter @ stringdemo.testa$1.<init>(testa.java:17) @ stringdemo.testa.test(testa.java:13) i using testng: @test public void test() throws exception { new expectations() { @mocked("length") string astring; { astring.length(); result = 2; } }; system.out.println(a.showa("test")); } } actual class a: public class { public static int showa(string str){ int a= str.length(); return a; } } this wrong way of recording expected results. shouldn't mock , record string's length() method, record showa() instead. here solution @suppresswarnings("unused&qu
Read more

asp.net - Razor Page Hosted on IIS 6 Fails Every Morning -

-
i have simple page in intranet uses razor/asp fetch single record table , display it, plus few graphic, dashboard display. meaning, no user intervention other first time open ie internal url. the problem having every morning ie displays "page not found or network error" message. refresh page, , it's stuck , doesn't display anything. try different pc, open ie internal url , it's stuck... ...until following: login server hosting page run inetmgr go web sites, etc right click on page that's giving me problem , select browse at moment, error message: server error in '/' application. -------------------------------------------------------------------------------- type of page not served. description: type of page have requested not served because has been explicitly forbidden. extension '.cshtml' may incorrect. please review url below , make sure spelled correctly. requested url: /application/dashboard.cshtml ---------------
Read more

c++ - wxwidget compiling on windows command prompt -

-
can please guide me on method,as how compile wxwidget using windows command prompt or msys. [i can't seem find on wxwiki/wxwidget official book/anywhere else] i have compiled wxwidget instructions have provided (using msys). from searching on internet seems 1 can (maybe?) through use of makefile (is correct?) if yes: directory & how should link wx libraries in makefile if no: way besides makefile? i know can use ide's code::block , make life simpler prefer compile using command prompt/msys. thanks in advance.
Read more