24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 1084  |  回复: 4

安的三

新虫 (初入文坛)

[求助] matlab中建立稀疏矩阵组并操作

1.这个稀疏矩阵共有9个,分别从3到11阶,只有紧邻对角线的上下两斜线有元素,请问如何
建立;然后分别求特征值与特征向量;
2.初学matlab,自己写的东西不知道哪里不对,大牛求解.

程序下付:
N=0:1:10;
%N as a vector
for k=2:10;
    A{k}=rand(k+1);
    B{k}=rand(k+1);
    C{k}=rand(k+1);
    x{k}=2:1:k+1;
    y{k}=1:1:k;
    for l=2:k+1;
        s(l)=sqrt((l+1)*(l-1));
    end
    A{k}=sparse(x{k},y{k},s);
  % C(k)=full(A(k));
    p{k}=1:1:k;
    q{k}=2:1:k+1;
    for n=2:k+1;
        r(n)=sqrt((n-1)*(n+1));
    B{k}=sparse(p{k},q{k},r);
   %D(k)=full(B(k));
    %F(k)=C(k)+D(k);
    %[V,E]=eig(F(k));
    C{k}=A{k}+B{k};
    [V,E]=eig(C(k));
end
V E

--
回复此楼
LinnaAnUSTCchemistry
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

csgt0

荣誉版主 (著名写手)

彩色挂图

【答案】应助回帖


感谢参与,应助指数 +1
安的三: 金币+1 2012-06-22 13:29:53
首先,x,y的长度为k,s的长度为k+1(s(1)你没给,就自动赋0了),而sparse需要三个一样长
其次A{k}和C{k}的大小不同,一个是m*n,一个是n*m
第三,C是cell,eig应该用C{k},且C{k}必须是方阵
showmethemoney
2楼2012-06-21 10:34:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

安的三

新虫 (初入文坛)


dubo: 金币+1, 欢迎讨论 2012-06-22 21:59:45
引用回帖:
2楼: Originally posted by csgt0 at 2012-06-21 10:34:51
首先,x,y的长度为k,s的长度为k+1(s(1)你没给,就自动赋0了),而sparse需要三个一样长
其次A{k}和C{k}的大小不同,一个是m*n,一个是n*m
第三,C是cell,eig应该用C{k},且C{k}必须是方阵

soga,可是我在一开始就把A确定成为K+1阶的仿真了啊?
LinnaAnUSTCchemistry
3楼2012-06-21 15:16:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

csgt0

荣誉版主 (著名写手)

彩色挂图

【答案】应助回帖


dubo: 金币+1, 欢迎讨论 2012-06-22 21:59:48
引用回帖:
3楼: Originally posted by 安的三 at 2012-06-21 15:16:51
soga,可是我在一开始就把A确定成为K+1阶的仿真了啊?...

比如这一句A{k}=sparse(x{k},y{k},s);
当k=2时
需要x{2},y{2},s一样长
但此时你的x{2},y{2}的长度为2
而s的长度为3
showmethemoney
4楼2012-06-21 15:42:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

安的三

新虫 (初入文坛)


dubo: 金币+1, 欢迎讨论 2012-06-22 21:59:51
引用回帖:
4楼: Originally posted by csgt0 at 2012-06-21 15:42:33
比如这一句A{k}=sparse(x{k},y{k},s);
当k=2时
需要x{2},y{2},s一样长
但此时你的x{2},y{2}的长度为2
而s的长度为3...

soga ,thanks
nice day!
LinnaAnUSTCchemistry
5楼2012-06-22 13:29:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 安的三 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿北交大材料工程总分358求调剂 +6 cs0106 2026-04-05 6/300 2026-04-06 08:48 by jj987
[考研] 一志愿211,化学学硕,310分,本科重点双非,求调剂 +13 努力奋斗112 2026-04-04 13/650 2026-04-06 07:13 by jj987
[考研] 308求调剂 +11 倘若起风了呢 2026-04-05 11/550 2026-04-05 23:21 by 来看流星雨10
[考研] (调剂)一志愿报考哈尔滨工业大学0857资源与环境专业378分考生 +7 狠狠加油 2026-04-05 7/350 2026-04-05 22:31 by dongzh2009
[考研] 材料专硕 调剂 +14 CXN123456 2026-04-03 14/700 2026-04-05 17:18 by Hdyxbekcb
[考研] 0703总分331求调剂 +10 ZY-05 2026-04-04 13/650 2026-04-05 11:03 by xiayan13521
[考研] 323分(计算机视觉和大模型项目)能直接上手 +3 chaoxiicy 2026-04-01 3/150 2026-04-05 00:50 by chongya
[考研] 一志愿郑大0705求调剂 +3 橘十一 2026-04-02 4/200 2026-04-05 00:05 by chongya
[考研] 材料与化工306分找调剂 +12 沧海轻舟e 2026-04-03 13/650 2026-04-04 23:45 by lqwchd
[考研] 324求调剂 +14 想上学求调 2026-04-02 15/750 2026-04-04 20:31 by 无际的草原
[考研] +5 雾与海 2026-04-02 6/300 2026-04-04 19:53 by 蓝云思雨
[考研] 294求调剂 +6 Grey_Ey 2026-04-03 6/300 2026-04-03 20:46 by 欣喜777
[考研] 301求调剂 +14 A_JiXing 2026-04-01 14/700 2026-04-03 18:31 by ls刘帅
[考研] 求材料调剂 一志愿南昌大学 328分 +5 yyy..... 2026-04-03 5/250 2026-04-03 13:46 by 百灵童888
[考研] 319求调剂 +18 太容易1018 2026-04-01 18/900 2026-04-03 11:18 by linyelide
[考研] 一志愿a区211,085601-307分求调剂 +13 党嘉豪 2026-03-31 26/1300 2026-04-03 08:33 by 495374996
[考研] 初试301,代码085701环境工程,本硕一致,四六级已过,有二区一作,共发表5篇论文 +6 axibli 2026-04-01 6/300 2026-04-02 13:42 by Ecowxq666!
[考研] 377求调剂 +3 RASKIN 2026-04-02 3/150 2026-04-02 09:45 by zzchen2000
[考研] 求调剂,一志愿北林食品与营养095500,301分,已过六级,有科研经历 +4 快乐储蓄罐 2026-03-31 4/200 2026-04-01 09:26 by JourneyLucky
[考研] 335求调剂 +3 321* 2026-03-31 4/200 2026-04-01 00:00 by 321*
信息提示
请填处理意见