24小时热门版块排行榜    

查看: 559  |  回复: 9
当前主题已经存档。

dreamred8087

金虫 (正式写手)

[交流] 【求助】求助大矩阵生成方法

求助大矩阵生成方法:由于规模太大,我的方法太慢了,求高手指点,谢谢了啊!

大矩阵 A  每行为: A(1,: ), A(2,: ),...,A(m,: )
先生成矩阵 C(i) = [ A(i,: )   0       0  ...   0;
                               0        A(i,: )  0  ...  0;
                              ..................................
                               0        0        0  ....  A(i,: )]
最后生成 C = [C(1);C(2);...;C(m)]

[len_m,len_n] = size(A);  

C = [];
for i=1:len_m
    a_sub = A(i,: );
        
    C_i = [];
     for j=1:len_n
        C_i_sub = sparse(1,(j-1)*len_n+1:j*len_n,a_sub,1,len_n*len_n);
        C_i     = [C_i;C_i_sub];
    end
end
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangen994

荣誉版主 (著名写手)

将军


dreamred8087(金币+1,VIP+0):谢谢支持! 7-5 17:04
占个沙发
2楼2009-07-05 08:32:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sitonyl

铁杆木虫 (文坛精英)

★ ★ ★
wangmonk(金币+1,VIP+0):能否有更多的指教呢? 7-16 21:57
sunxiao(金币+1,VIP+0):谢谢参与 7-19 03:15
dreamred8087(金币+1,VIP+0):谢谢! 7-19 11:40
应该用稀疏矩阵
zz
3楼2009-07-16 21:49:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

woshilsh

荣誉版主 (职业作家)

优秀版主


dreamred8087(金币+1,VIP+0):共同期待 7-19 11:41
期待高手指导!
[center][url=http://www.91cool.net/][img]http://id.91cool.net/sign/?name=小木虫印&say=各位版主辛苦了![/img][/url][/center]
4楼2009-07-16 22:27:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xuqingyuan1972

金虫 (正式写手)

★ ★ ★
sunxiao(金币+2,VIP+0):此方法可行,能快速生成矩阵 7-19 03:15
dreamred8087(金币+1,VIP+0):谢谢! 7-19 11:41
先生成数据文件,再用load命令,最后直接生成矩阵!
5楼2009-07-18 09:32:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

羽寒月

银虫 (小有名气)

★ ★ ★
sunxiao(金币+2,VIP+0):解答不错,谢谢,欢迎常来 7-19 03:15
dreamred8087(金币+1,VIP+0):谢谢! 7-19 11:41
矩阵0太多了,还是用稀疏矩阵吧
sparse(i, j, s); %其中i, j, s 都是列矢量, C(i, j) = s
6楼2009-07-18 19:24:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunxiao

荣誉版主 (著名写手)

小木虫钓鱼协会会长


dreamred8087(金币+1,VIP+0):谢谢! 7-19 11:41
补充楼上的,如果i,j相同,那么所对应的s值会相加
西南交大隧道方向招硕、博士
7楼2009-07-19 03:16:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wen5920708

木虫 (正式写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
我一般都等着。。。确实很慢,期待指点。。。
8楼2009-07-19 12:23:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hitzhang

木虫 (正式写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
wangen994(金币+1,VIP+0):感谢参与讨论 8-12 09:55
循环太多了,另外矩阵最好事先成赋0值,最好是稀疏矩阵。试试blkdiag,vertcat,horzcat这几个函数可能会快一些。
9楼2009-08-11 19:43:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dreamred8087

金虫 (正式写手)

非常感谢楼上的意见,我会按照楼上的意见试试的!
10楼2009-08-12 13:34:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 dreamred8087 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见