查看: 547|回复: 1|关注: 0

[未答复] 利用sparse创建稀疏矩阵时报错,求指正

[复制链接]

新手

7 麦片

财富积分


050


6

主题

11

帖子

0

最佳答案
发表于 2020-1-16 22:55:21 | 显示全部楼层 |阅读模式
想利用sparse函数建立一个上三角形(NS-3)*(NS-1)稀疏矩阵,主对角线上元素是一个函数值向量c1。c2、c3也是函数值向量。
  1. 9    al = PML_model.a(x(2:NS-2));
  2. 10   dal = PML_model.da(x(2:NS-2));
  3. 11   c1 = -1/(h^2)./al-dal./h/(al).^2-al.*k^2;
  4. 12   c2 = 2/(h^2)./al+dal./h/(al).^2;
  5. 13   c3 = -1/(h^2)./al;
  6. 14   A = sparse(1:NS-3,1:NS-3,c1,NS-3,NS-1) + sparse(1:NS-3,2:NS-2,c2,NS-3,NS-1)...
  7. 15   + sparse(1:NS-3,3:NS-1,c3,NS-3,NS-1);
复制代码


                               
登录/注册后可看大图



报错如下
错误使用 sparse
矢量长度必须相同。


出错 PML_FD1d_bvp (line 14)
A = sparse(1:NS-3,1:NS-3,c1,NS-3,NS-1) + sparse(1:NS-3,2:NS-2,c2,NS-3,NS-1)...

回复主题 已获打赏: 0 积分

举报

新手

7 麦片

财富积分


050


6

主题

11

帖子

0

最佳答案
 楼主| 发表于 2020-2-14 12:17:35 | 显示全部楼层
差分出了问题,已经解决了:victory:
回复此楼 已获打赏: 0 积分

举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

站长推荐上一条 /4 下一条

快速回复 返回顶部 返回列表