%%稀疏矩陣左除bug
L=diag(3*ones(5,1))+diag(2*ones(4,1),-1)+diag(2*ones(4,1),1);
A=blkdiag(L,L,L,L,L)+diag(20*ones(20,1),5)+diag(20*ones(20,1),-5)
b=ones(25,1);
B=sparse(A);
x=B\b
這個求解稀疏矩陣的例子,在不使用sparse函數(shù)的時候可以得到正常結(jié)果,取sparse就會報錯。
在命令行窗口使用enter或者shift+enter均無法換行
L=diag(3*ones(5,1))+diag(2*ones(4,1),-1)+diag(2*ones(4,1),1);
A=blkdiag(L,L,L,L,L)+diag(20*ones(20,1),5)+diag(20*ones(20,1),-5);
B=sparse(A);
log(abs(B(1)))
p = struct();
p.var(1).name=1;
disp(p.var(1).name)
直接運行這段代碼不會得到結(jié)果也不會報錯;調(diào)試這段代碼會卡在第二行,無法繼續(xù)運行
%%稀疏矩陣右除bug
L=diag(3*ones(5,1))+diag(2*ones(4,1),-1)+diag(2*ones(4,1),1);
A=blkdiag(L,L,L,L,L)+diag(20*ones(20,1),5)+diag(20*ones(20,1),-5);
B=sparse(A);
x=B/5
報錯信息顯示計算未支持,請將稠密矩陣轉(zhuǎn)為稀疏矩陣。但是,將B轉(zhuǎn)為full(B)計算可得到正確結(jié)果。
L=diag(3*ones(5,1))+diag(2*ones(4,1),-1)+diag(2*ones(4,1),1);
A=blkdiag(L,L,L,L,L)+diag(20*ones(20,1),5)+diag(20*ones(20,1),-5);
B=sparse(A);
power(abs(B(1)),5)
n = 10;
e = ones(n,1);
A = spdiags([-1*e 2*e -1*e], -1:1, n, n)
B = (1:n)' * (1:n);
X = A \ B
在這個例子中,直接計算不會報錯但會得到錯誤的結(jié)果。使用full(A)能得到正確結(jié)果
當(dāng)前版本暫未支持vpa函數(shù)
給一個沒有預(yù)定義的五維數(shù)組賦值會提示左側(cè)的索引與右側(cè)的大小不兼容。預(yù)定義之后能得到正確結(jié)果。
ii_vg = 1;
ii_vd = 1;
Mu_sub = rand(2,3,4);
Mu_sub_body(:,:,:,ii_vg,ii_vd) = Mu_sub;
disp(size(Mu_sub_body))
將圖窗保存到文件的print函數(shù)未支持
save函數(shù)目前僅支持-struct格式,不支持-ascii格式
目前北太無法打開IN格式的輸入文件