查看: 81|回复: 2|关注: 0

[已答复] 求大神debug试图沿模糊的维增大数组

[复制链接]

新手

10 麦片

财富积分


050


1

主题

2

帖子

0

最佳答案
发表于 6 天前 | 显示全部楼层 |阅读模式
代码如下:
  1. length=2;
  2. N=100;iterN=0;eps=1.0e-4;deltQ=2;
  3. T2=(1:N);
  4. T2(1)=8;
  5. T2(N)=1;
  6. dx=length/(N-2);
  7. while(deltQ>=eps && iterN<50)
  8.     T=T2;
  9. for i=1:N
  10.     k(i)=3*T(i)^3;
  11.     sc(i)=20*(i-1)/N;
  12.     sp(i)=4;
  13. end
  14. k(N/2+1:N)=40*T(N/2+1:N).^3;
  15. for i=2:N-1
  16.     if i==2
  17.         a(i-1)=-2*k(i)/dx;
  18.         c(i-1)=-2*k(i)*k(i+1)/dx/(k(i)+k(i+1));
  19.         b(i-1)=-a(i-1)-c(i-1)-sp(i)*dx;
  20.         d(i-1)=sc(i)*dx-a(i-1)*T(i-1);
  21.     else
  22.         if i==N-1
  23.             a(i-1)=-2*k(i)*k(i-1)/dx/(k(i)+k(i-1));
  24.             c(i-1)=-2*k(i+1)/dx;
  25.             b(i-1)=-a(i-1)-c(i-1)-sp(i)*dx;
  26.             d(i-1)=sc(i)*dx-c(i-1);
  27.         else
  28.             a(i-1)=-2*k(i)*k(i-1)/dx/(k(i)+k(i-1));
  29.             c(i-1)=-2*k(i)*k(i+1)/dx/(k(i)+k(i+1));
  30.             b(i-1)=-a(i-1)-c(i-1)-sp(i)*dx;
  31.             d(i-1)=sc(i-1)*dx;
  32.         end
  33.     end
  34. end
  35. if abs(b(1))<=abs(c(1)) || abs(b(N-2)<=abs(c(N-2)))
  36.     warndla('不满足追赶法条件','警告');
  37. end
  38. for i=1:N-2
  39.     if abs(b(i))<abs(c(i)+abs(a(i)))
  40.         warndla('不满足追赶法条件','警告');
  41.     end
  42. end
  43. for i=1:N-2
  44.     if i==1
  45.         u(1)=b(1);
  46.     else
  47.         l(i)=a(i)/u(i-1);
  48.         u(i)=b(i)-l(i)*c(i-1);
  49.     end
  50. end
  51. for i=1:N-2
  52.     if i==1
  53.         y(1)=d(1);
  54.     else
  55.         y(i)=d(i)-l(i)*y(i-1);
  56.     end
  57. end
  58. for j=N-2:-1:1
  59.     if j==N-2
  60.         x(N-2)=y(N-2)/u(N-2);        %该句报错‘试图沿模糊的维增大数组’
  61.     else
  62.         x(j)=(y(j)-c(j)*x(j+1))/u(j);
  63.     end
  64. end
  65.     iterN=iterN+1;
  66.     T2(2:N-1)=x(1:N-2);
  67.     deltQ=max(abs(T2-T));
  68. end
  69. for i=1:10
  70.     Q(i,1:10)=T2((1+(i-1)*10):i*10);
  71. end
  72. disp(Q);
复制代码
不知道为什么报错,求解释,如果可以的话,帮忙给出改法吧,谢谢大家
回复主题 已获打赏: 0 积分

举报

论坛优秀回答者

9

主题

1629

帖子

342

最佳答案
  • 关注者: 81
发表于 6 天前 | 显示全部楼层
在2020a运行正常无报错
1.PNG
回复此楼 已获打赏: 0 积分

举报

新手

10 麦片

财富积分


050


1

主题

2

帖子

0

最佳答案
 楼主| 发表于 6 天前 | 显示全部楼层
啊?!我是2017a,这个错误是真的有版本问题啊,谢谢谢谢
回复此楼 已获打赏: 0 积分

举报

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

本版积分规则

关闭

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

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