作者edcamel (大坏蛋)
看板MATLAB
标题[问题] sparse与偏微方
时间Thu Jan 26 22:49:22 2012
我用Octave
看了help和GNU上的教学文件後，依然不懂sparse怎麽用
假设今天K是一个N*N矩阵，然後
K(i,i)=2 for i=1:N; K(i-1,i)=K(i,i-1)=-1 for i=2:N;
我要造一个(N^2)*(N^2)的矩阵，叫K2D： (I=eye(N))
K2D=kron(K,I)+kron(I,K);
请问，我要怎麽告诉Octave K2D是一个sparse矩阵?
因为我试了先造K2D，再使用sparse(K2D)，耗时有点久，很怀疑自己搞错
然後我要解
K2D*u=F， F和u是N*1的vector, F已知，要解u
我查了一些资料，说是如果在sparse的情况下，直接backslash F\K2D几乎会是最快的
方法
但我实在不懂怎麽利用sparse啊...
先谢谢大家的指教!
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 61.56.129.110
※ 编辑: edcamel 来自: 61.56.129.110 (01/26 22:52)
1F：推 sunev:octave没有针对sparse做处理吗? matlab是没问题的 01/27 02:26
2F：→ sunev:你可以参考matlab的kron的写法，再转成octave 01/27 02:27