我怎么解决”指数在位置2超出数组边界(不得超过1)。”

24日视图(30天)

显示旧的评论

撒母耳Suakye 2023年7月10日

0
链接

这个问题直接联系

//www.tianjin-qmedu.com/matlabcentral/answers/1994123-how-do-i-fix-index-in-position-2-exceeds-array-bounds-must-not-exceed-1

评论道: 撒母耳Suakye2023年7月11日17:47

                          clc;
                         
                          清晰;
                         
                          关闭所有;
                         
                          % %数据
                         
                          %定义参数
                         
                          量化宽松政策= 1;
                         
                          n = 5 e16天;
                         
                          p = 0.525;
                         
                          E_0 = 10 e3;
                         
                          v_s = 3.5 e3;
                         
                          epsilon_d = 50;
                         
                          问题= 1;
                         
                          m0 = 9.1094 e-31;
                         
                          m = 0.44 * m0;
                         
                          epsilon_0 = 1;
                         
                          h = 4.1356 e15汽油;
                         
                          百巴= h /(2 *π);
                         
                          b = 0.142 e-9;
                         
                          = 3 * b /(2 *百巴);
                         
                          tt = 2.6;
                         
                          δ= 5.7 / 2;
                         
                          G =√△^ 2 + 5 * tt ^ 2);
                         
                          %为ω_p定义范围和B ^ (1)
                         
                          omega_p = linspace (0, 150);
                         
                          B_inv = linspace (0.04, 150);
                         
                          k_q = omega_p / v_s;
                         
                          v = (2 * tt ^ 2 * * sin (p *) / G);
                         
                          τ=问题* p / v (qe *);
                         
                          tau_1 =τ/ 2;
                         
                          omega_c = qe * B_inv‘* v / p;
                         
                          γ= 1 + 1我* omega_c * tau_1;
                         
                          sigma_g = qe ^ 2 * n * v *τ/ p;
                         
                          6 = p ^ 2 /(2 *π*百巴^ 2);
                         
                          里克斯= 2 * p / m / v;
                         
                          b_g = 2 *π*百巴^ 2 * epsilon_0 * v / (qe ^ 2 * p);
                         
                          beta_g = (k_q”。* v) / omega_c;
                         
                          %初始化结果矩阵
                         
                          j_x = 0(元素个数(omega_p),元素个数(B_inv));
                         
                          %计算sigma_xx
                         
                          sigma_xx = 0(元素个数(omega_p),元素个数(B_inv));
                         
                          为r =三3
                         
                          J_r = besselj (r, beta_g);
                         
                          sigma_xx = sigma_xx + (2 * sigma_g ' * J_r) / (1 - 1 * (omega_p - r * omega_c) *τ);
                         
                          结束
                         
                          %计算R_x
                         
                          R_x = 0(元素个数(omega_p),元素个数(B_inv));
                         
                          为r =三3
                         
                          J_r = besselj (r, beta_g);
                         
                          R_x = R_x + (omega_c。/ k_q)。* (J_r ^ 2)。/ (1 - 1 * (omega_p - r * omega_c) *τ);
                         
                          结束
                         
                          %计算g (omega_q k_q)
                         
                          g_omega_q_k_q = 1 + 1 / (epsilon_0 * (epsilon_d + 1) * (sigma_xx。/ (v_s - R_x));
                         
                          %计算方程为每个omega_p和B ^ 1
                         
                          为i = 1:元素个数(omega_p)
                         
                          为j = 1:元素个数(B_inv)
                         
                          sum_term = 0;
                         
                          为r = 1:1
                         
                          J_r = besselj (r, beta_g(我));
                         
                          sum_term = sum_term + J_r / (1 - 1 * (omega_p (i) - r * omega_c (j)) *τ);
                         
                          结束
                         
                          j_x (i, j) = (1 / (2 * qe * n * v (i))) *((τ* omega_c (j) * p ^ 2) /(8 *π*(百巴^ 2)))* (abs ((sigma_g (i, j) * E_0) / g_omega_q_k_q (i, j))) ^ 2…
                         
                          * ((1 / beta_g(我))/ (1 + (omega_c (j) ^ 2 * tau_1 ^ 2))) ^ 2 * (sum_term)…
                         
                          *(我*γ(i, j) ^ 2 * (sum_term + 1) * besselj (sum_term + 1, beta_g(我))+ 1我*连词(γ(i, j)) ^ 2 * (sum_term - 1) * besselj (sum_term - 1, beta_g(我)));
                         
                          结束
                         
                          结束
                         
                             指数在位置2超出数组边界。指标不得超过1。

                          %绘制结果
                         
                          [Omega_P, B_inv] = meshgrid (Omega_P B_inv);
                         
                          冲浪(Omega_P B_inv j_x);
                         
                          包含(“\ omega_p”);
                         
                          ylabel (" B ^ {1}”);
                         
                          zlabel (“j_x”);
                         
                          标题(的情节j_x反对\ omega_p和B ^ {1}”);
                         
                          标题(的情节j_x反对\ omega_p和B ^ {1}”);

0评论
显示1年长的评论藏1年长的评论

登录置评。

在回答这个问题。

答案(1)

Torsten 2023年7月10日17:57

0
链接

直接链接到这个答案

//www.tianjin-qmedu.com/matlabcentral/answers/1994123-how-do-i-fix-index-in-position-2-exceeds-array-bounds-must-not-exceed-1 answer_1270628

编辑:Torsten 2023年7月10日在21:16

在进入循环后,弹出错误消息类型

                              大小(γ)
                             
                              大小(sigma_g)
                             
                              大小(g_omega_q_k_q)

并检查是否这些数组的第二个维度是> =元素个数(B_inv)是必需的,因为下面的循环。

9日评论
显示8年长的评论隐藏8年长的评论

Torsten 2023年7月11日在12:54

编辑:Torsten 2023年7月11日12:56

像这样的吗?“伽马”仍有一个错误的二维度和sigma_g错第一维度。

                                   clc;
                                  
                                   清晰;
                                  
                                   关闭所有;
                                  
                                   % %数据
                                  
                                   %定义参数
                                  
                                   量化宽松政策= 1;
                                  
                                   n = 5 e16天;
                                  
                                   p = 0.525;
                                  
                                   E_0 = 10 e3;
                                  
                                   v_s = 3.5 e3;
                                  
                                   epsilon_d = 50;
                                  
                                   问题= 1;
                                  
                                   m0 = 9.1094 e-31;
                                  
                                   m = 0.44 * m0;
                                  
                                   epsilon_0 = 1;
                                  
                                   h = 4.1356 e15汽油;
                                  
                                   百巴= h /(2 *π);
                                  
                                   b = 0.142 e-9;
                                  
                                   = 3 * b /(2 *百巴);
                                  
                                   tt = 2.6;
                                  
                                   δ= 5.7 / 2;
                                  
                                   G =√△^ 2 + 5 * tt ^ 2);
                                  
                                   %为ω_p定义范围和B ^ (1)
                                  
                                   omega_p = linspace (0, 150);
                                  
                                   B_inv = linspace (0.04, 150);
                                  
                                   k_q = omega_p / v_s;
                                  
                                   v = (2 * tt ^ 2 * * sin (p *) / G);
                                  
                                   τ=问题* p / v (qe *);
                                  
                                   tau_1 =τ/ 2;
                                  
                                   omega_c = qe * B_inv‘* v / p;
                                  
                                   γ= 1 + 1我* omega_c * tau_1;
                                  
                                   sigma_g = qe ^ 2 * n * v *τ/ p;
                                  
                                   6 = p ^ 2 /(2 *π*百巴^ 2);
                                  
                                   里克斯= 2 * p / m / v;
                                  
                                   b_g = 2 *π*百巴^ 2 * epsilon_0 * v / (qe ^ 2 * p);
                                  
                                   beta_g = (k_q”。* v) / omega_c;
                                  
                                   %初始化结果矩阵
                                  
                                   j_x = 0(元素个数(omega_p),元素个数(B_inv));
                                  
                                   %计算sigma_xx
                                  
                                   sigma_xx = 0(元素个数(omega_p),元素个数(B_inv));
                                  
                                   为r =三3
                                  
                                   J_r = besselj (r, beta_g);
                                  
                                   sigma_xx = sigma_xx + (2 * sigma_g ' * J_r) / (1 - 1 * (omega_p - r * omega_c) *τ);
                                  
                                   结束
                                  
                                   %计算R_x
                                  
                                   R_x = 0(元素个数(omega_p),元素个数(B_inv));
                                  
                                   为r =三3
                                  
                                   J_r = besselj (r, beta_g);
                                  
                                   R_x = R_x + (omega_c。/ k_q)。* (J_r ^ 2)。/ (1 - 1 * (omega_p - r * omega_c) *τ);
                                  
                                   结束
                                  
                                   %计算g (omega_q k_q)
                                  
                                   g_omega_q_k_q = 1 + 1 / (epsilon_0 * (epsilon_d + 1) * (sigma_xx。/ (v_s - R_x));
                                  
                                   %调整伽马和sigma_g维度
                                  
                                   γ= repmat(1 + 1我* omega_c * tau_1元素个数(omega_p), 1);
                                  
                                   sigma_g = repmat (qe ^ 2 * n * v * t / p, 1,元素个数(B_inv));
                                  
                                   %计算方程为每个omega_p和B ^ 1
                                  
                                   大小(γ)
                                  
                                      ans =
                                      1×2
                                     
                                       22500年1

                                   大小(sigma_g)
                                  
                                      ans =
                                      1×2
                                     
                                       1 150

                                   大小(g_omega_q_k_q)
                                  
                                      ans =
                                      1×2
                                     
                                       150 150

                                   元素个数(B_inv)
                                  
                                     ans = 150

                                   元素个数(omega_p)
                                  
                                     ans = 150

                                   为i = 1:元素个数(omega_p)
                                  
                                   为j = 1:元素个数(B_inv)
                                  
                                   sum_term = 0;
                                  
                                   为r = 1:1
                                  
                                   J_r = besselj (r, beta_g(我));
                                  
                                   sum_term = sum_term + J_r / (1 - 1 * (omega_p (i) - r * omega_c (j)) *τ);
                                  
                                   结束
                                  
                                   j_x (i, j) = (1 / (2 * qe * n * v (i))) *((τ* omega_c (j) * p ^ 2) /(8 *π*(百巴^ 2)))* (abs ((sigma_g (i, j) * E_0) / g_omega_q_k_q (i, j))) ^ 2…
                                  
                                   * ((1 / beta_g(我))/ (1 + (omega_c (j) ^ 2 * tau_1 ^ 2))) ^ 2 * (sum_term)…
                                  
                                   *(我*γ(i, j) ^ 2 * (sum_term + 1) * besselj (sum_term + 1, beta_g(我))+ 1我*连词(γ(i, j)) ^ 2 * (sum_term - 1) * besselj (sum_term - 1, beta_g(我)));
                                  
                                   结束
                                  
                                   结束
                                  
                                      指数在位置2超出数组边界。指标不得超过1。

                                   指数在位置2超出数组边界。指标不得超过1。
                                  
                                   %绘制结果
                                  
                                   [Omega_P, B_inv] = meshgrid (Omega_P B_inv);
                                  
                                   冲浪(Omega_P B_inv j_x);
                                  
                                   包含(“\ omega_p”);
                                  
                                   ylabel (" B ^ {1}”);
                                  
                                   zlabel (“j_x”);
                                  
                                   标题(的情节j_x反对\ omega_p和B ^ {1}”);
                                  
                                   标题(的情节j_x反对\ omega_p和B ^ {1}”);

撒母耳Suakye 2023年7月11日17:47

我想计算j_x和策划一个3 d图形

登录置评。

在回答这个问题。

类别

人工智能,数据科学和统计数据统计和机器学习工具概率分布连续分布均匀分布(连续)

找到更多的在均匀分布(连续)在帮助中心和文件交换

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!

我怎么解决”指数在位置2超出数组边界(不得超过1)。”

0评论
显示1年长的评论藏1年长的评论

答案(1)

9日评论
显示8年长的评论隐藏8年长的评论

另请参阅

类别

标签

社区寻宝

我怎么解决”指数在位置2超出数组边界(不得超过1)。”

0评论 显示1年长的评论藏1年长的评论

答案(1)

9日评论 显示8年长的评论隐藏8年长的评论

另请参阅

类别

标签

社区寻宝

0评论
显示1年长的评论藏1年长的评论

9日评论
显示8年长的评论隐藏8年长的评论