%%稀疏矩陣左除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
這個(gè)求解稀疏矩陣的例子,在不使用sparse函數(shù)的時(shí)候可以得到正常結(jié)果,取sparse就會(huì)報(bào)錯(cuò)。
在命令行窗口使用enter或者shift+enter均無(wú)法換行
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)
直接運(yùn)行這段代碼不會(huì)得到結(jié)果也不會(huì)報(bào)錯(cuò);調(diào)試這段代碼會(huì)卡在第二行,無(wú)法繼續(xù)運(yùn)行
%%稀疏矩陣右除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
報(bào)錯(cuò)信息顯示計(jì)算未支持,請(qǐng)將稠密矩陣轉(zhuǎn)為稀疏矩陣。但是,將B轉(zhuǎn)為full(B)計(jì)算可得到正確結(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
在這個(gè)例子中,直接計(jì)算不會(huì)報(bào)錯(cuò)但會(huì)得到錯(cuò)誤的結(jié)果。使用full(A)能得到正確結(jié)果
當(dāng)前版本暫未支持vpa函數(shù)
給一個(gè)沒(méi)有預(yù)定義的五維數(shù)組賦值會(huì)提示左側(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))
運(yùn)行如下代碼程序會(huì)直接閃退,
global A B X;
linspace(A,B,X)
期望得到的結(jié)果是報(bào)錯(cuò)(如輸入必須為標(biāo)量)。
將圖窗保存到文件的print函數(shù)未支持
save函數(shù)目前僅支持-struct格式,不支持-ascii格式
目前北太無(wú)法打開(kāi)IN格式的輸入文件