How can i write the eq. in attached file in matlab ? please help.....

3 次查看(过去 30 天)
this equation is obtained by finite-difference scheme method.
  4 个评论
Walter Roberson
Walter Roberson 2015-12-8
You need boundary conditions. The form of the boundary conditions will give me more information about how to approach it.
Prince
Prince 2015-12-8
编辑:Prince 2015-12-8
A two-point boundary can be posed the orientation director versus the z-coordinate, and it’s controlled by the equation. Dirichlet type used the two boundary conditions for strong anchoring:- θ(z=0)=90 , θ(z=L)=90 where: L =0.2

请先登录,再进行评论。

采纳的回答

Walter Roberson
Walter Roberson 2015-12-8
编辑:Walter Roberson 2015-12-8
Based upon the information you have given, the equation is
theta(k,i) = 0
for all k and i.
  3 个评论
Walter Roberson
Walter Roberson 2015-12-8
In email, you asked for the expression anyhow, saying that you will find a way to make the rest work. The expression is
theta(k+1, i) = (1/3) * (-(theta(k, i+1) + theta(k, i) + theta(k, i-1)) + 1/6*(V^2/Vth^2)*pi^2*h^2*sin(2*theta(k, i)))
However, if Vth is intended as having a relationship to V then this might need to change a bit.
The validity of the expression depends upon how the boundary conditions are expressed, whether they start from theta(k,0) or from theta(k,1) or theta(k,2); the expression would also need to change if the first k was 0 instead of 1.
You expressed a boundary condition on theta(z=0) without giving any way to interpret that in terms of theta(k,i) . If we guess that maybe you want to start numbering i from 0, then we run into the problem that your other boundary condition is theta(z=0.2) and certainly 0.2 is not a valid array index. Is the idea that you divided the length L up into pieces indexed at i? Or indexed at k? We do not know, but you are going to have to know if you want to make use of the expression you asked for.

请先登录,再进行评论。

更多回答(0 个)

类别

Help CenterFile Exchange 中查找有关 Matrix Indexing 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by