S**********r
发帖数: 28 | 1 小女子打算求一个ODE的特征值：
(d^4/dt^4) y(t)=lambda* y(t),t in [0,1]
边值条件是y'(0)=y'(1)=y''(0)=y''(1)=0
据说有个很神奇的软件叫chebfun()可以做，小女子于是编了下面这个code
trunc=30;
tic
L = chebop(@(x,u) diff(u,4), [0,1]);
L.lbc = @(u) diff(u,1);
L.lbc =@(u) diff(u,2);
L.rbc = @(u) diff(u,1);
L.rbc =@(u) diff(u,2);
[V,D] = eigs(L,trunc);
disp(diag(D)), toc
可是有报错说
Error using eig
Matrix must be square.
各位计算数学大侠请帮忙看看小女子哪里错了？谢谢！ | w**a
发帖数: 1024 | | N***m
发帖数: 4460 | 3 找个男人替你去做
【在 S**********r 的大作中提到】
: 小女子打算求一个ODE的特征值：
: (d^4/dt^4) y(t)=lambda* y(t),t in [0,1]
: 边值条件是y'(0)=y'(1)=y''(0)=y''(1)=0
: 据说有个很神奇的软件叫chebfun()可以做，小女子于是编了下面这个code
: trunc=30;
: tic
: L = chebop(@(x,u) diff(u,4), [0,1]);
: L.lbc = @(u) diff(u,1);
: L.lbc =@(u) diff(u,2);
: L.rbc = @(u) diff(u,1);
| E*****T
发帖数: 1193 | | d********8
发帖数: 3019 | 5 我刚好认识一个人刚做过这个这个方程，同样的特征值方程和变边界条件，不过他是个
男生。如果需要他的联系方式，请pm给我。
----
发自 MITBBS快读 for Android Tablet/Phone v3.2.3
https://market.android.com/details?id=com.mitbbsreader.app
【在 S**********r 的大作中提到】
: 小女子打算求一个ODE的特征值：
: (d^4/dt^4) y(t)=lambda* y(t),t in [0,1]
: 边值条件是y'(0)=y'(1)=y''(0)=y''(1)=0
: 据说有个很神奇的软件叫chebfun()可以做，小女子于是编了下面这个code
: trunc=30;
: tic
: L = chebop(@(x,u) diff(u,4), [0,1]);
: L.lbc = @(u) diff(u,1);
: L.lbc =@(u) diff(u,2);
: L.rbc = @(u) diff(u,1);
| R*******T
发帖数: 1171 | 6 小女子打算求一个ODE的特征值：
(d^4/dt^4) y(t)=lambda* y(t),t in [0,1]
边值条件是y'(0)=y'(1)=y''(0)=y''(1)=0
据说有个很神奇的软件叫chebfun()可以做，小女子于是编了下面这个code
trunc=30;
tic
L = chebop(@(x,u) diff(u,4), [0,1]);
L.lbc = @(u)[ diff(u,1), diff(u,2)];
L.rbc = @(u)[ diff(u,1), diff(u,2)]
[V,D] = eigs(L,trunc);
disp(diag(D)), toc
试试这个，不知道对不对。 |
|
