matlab - Error in Ordinary Differential Equation representation -

-

update trying find lyapunov exponents given in link le. trying figure out , understand taking following eqs case. these set of ordinary differential equations (these testing how work cos , sin ode)

f(1)=alpha*(y-x); f(2)=x*(r-z)-y; f(3) = 10*cos(x);

and x=x(1); y=x(2); cos(y)=x(3);

f1 means dx/dt;f2 dy/dt , f3 in case -10sinx. however,when expressing x=x(1);y=x(2);i unsure how express cos.this trial example doing know how work equations have cos,sin etc terms function of variable.

when using ode45 solve these eqs 

[t,res]=sol(3,@test_eq,@ode45,0,0.01,20,[7 2 100 ],10);

it throws following error

??? attempted access (2); index must positive integer or logical.

error in ==> eq @ 19 x=x(1); y=x(2); cos(x)=x(3);
  1. is representation x=x(1); y=x(2); cos(y)=x(3); alright?
  2. how resolve error? thank you

no representation wrong. can't possibly set values in way! start, trying assign value x(3) function.

first not sure understand difference between

x=4

and

4=x

which different meanings. if understand this, you'll see can't possibly assign using cos(x)=x(3).

second: function sol() calling? have defined it?

third, solve or evaluate odes should using deval or solve functions in matlab. files should give examples.


Comments

Post a Comment

Popular posts from this blog

c# - SVN Error : "svnadmin: E205000: Too many arguments" -

-
i trying reposotries using c# code process svncommand = null; var psi = new processstartinfo("svnadmin"); psi.redirectstandardoutput = true; psi.redirectstandarderror = true; psi.useshellexecute = false; psi.arguments= @"dump c:\repositories\myrepo > c:\temp\myrepodumpfile.dump"; using (svncommand = process.start(psi)) { var myoutput = svncommand.standardoutput; var myerror = svncommand.standarderror; debug.write("output :" + environment.newline + environment.newline + myoutput.readtoend() + environment.newline + "error :" + environment.newline + environment.newline + myerror.readtoend()); svncommand.close(); } when use dump command commandline svnadmin dump c:\repositories\myrepo > c:\temp\myrepodumpfile.dump it works fine when try ...
Read more

c++ - Using OpenSSL in a multi-threaded application -

-
i have been writing soap client application in c++ on ubuntu using openssl https transport , pthreads threading. have number of threads - 1 central data acquisition thread periodically gets worker threads make soap requests via shared mutex protected queues. reading documentation openssl found is openssl thread-safe? in openssl faq describes mechanisms required ensure thread safety when using openssl. implemented , works fine. the reason question conceptual difficulty really. thinking implementing same functionality application has already, instead of using threads, create 2 seperate applications : 1 worker thread (of multiple copies running) , main data acquisition thread. use tcp sockets communicate between 2 rather mutex protected queues. may bad idea not important - confusing me would have implement same functions required ensure openssl thread safety in second approach or not? my guess not have , treated independent (indeed surely must many applications use openssl) re...
Read more

All overlapping substrings matching a java regex -

-
is there api method returns (possibly overlapping) substrings match regular expression? for example, have text string: string t = 04/31 412-555-1235; , , have pattern: pattern p = new pattern("\\d\\d+"); matches strings of 2 or more characters. the matches are: 04, 31, 412, 555, 1235. how overlapping matches? i want code return: 04, 31, 41, 412, 12, 55, 555, 55, 12, 123, 1235, 23, 235, 35. theoretically should possible -- there obvious o(n^2) algorithm enumerates , checks substrings against pattern. edit rather enumerating substrings, safer use region(int start, int end) method in matcher . checking pattern against separate, extracted substring might change result of match (e.g. if there non-capturing group or word boundary check @ start/end of pattern). edit 2 actually, it's unclear whether region() expect zero-width matches. specification vague, , experiments yield disappointing results. for example: string line = "xx90xx"; strin...
Read more