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# - Copy ObservableCollection to another ObservableCollection -

-
how copy observablecollection item observablecollection without reference of first collection? here observablecollection item value changes affecting both collections. code private observablecollection<ratemodel> _allmetalrate = new observablecollection<ratemodel>(); private observablecollection<ratemodel> _metalrateondate = new observablecollection<ratemodel>(); public observablecollection<ratemodel> allmetalrate { { return this._allmetalrate; } set { this._allmetalrate = value; notifypropertychanged("metalrate"); } } public observablecollection<ratemodel> metalrateondate { { return this._metalrateondate; } set { this._metalrateondate = value; notifypropertychanged("metalrateondate"); } } foreach (var item in metalrateondate) allmetalrate.add(item); what causing , how can solve it? you need clone object referenced item before it's a...
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