function poisson(field, prm)% Calculate Ex from Poisson equationf = field; % reference to the Field objf.ex(prm.X2) = f.ex(prm.X2-1) + f.rho(prm.X2);ex0 = sum(f.ex(prm.X2))/prm.nx %平均電場強(qiáng)度f.ex(prm.X2) = f.ex(prm.X2) - ex0; %？f.ex(1) = f.ex(prm.nxp1); %nxp1=nx+1f.ex(prm.nxp2) = f.ex(2); end >> prm.X2ans =Columns 1 through 162 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17Columns 17 through 3218 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33Columns 33 through 4834 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49Columns 49 through 6450 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65Columns 65 through 8066 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81Columns 81 through 9682 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97Columns 97 through 11298 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113Columns 113 through 128114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129>> f.ex Unable to resolve the name f.ex.>> field.exans =0000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000>> field.ex(prm.X2)ans =00000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000>> field.ex(prm.X2)ans =00000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000>> field.rhoans =1.0e-04 *-0.01610.0068-0.0272-0.02320.02960.02860.0298-0.0726-0.04510.0886-0.06980.1098-0.0172-0.05940.0846-0.05150.00560.0295-0.04370.0732-0.0528-0.05210.1384-0.16000.05710.00280.0386-0.01290.0232-0.05710.0765-0.13430.01530.0398-0.05140.01850.0307-0.05800.03200.01500.00780.0365-0.00910.0169-0.07020.0768-0.0393-0.0055-0.0171-0.0179-0.00300.0138-0.04120.08240.0326-0.0084-0.0116-0.08610.04940.00480.0140-0.04960.0690-0.0408-0.0293-0.04690.1150-0.0627-0.0000-0.04810.07430.0403-0.02000.01640.0429-0.10610.00340.0407-0.08320.06320.0992-0.0657-0.0033-0.03440.0497-0.06510.05640.0068-0.07260.0799-0.02490.0238-0.06030.01910.0342-0.08870.0843-0.05530.0564-0.02010.01760.0327-0.13490.0326-0.00020.0892-0.09870.0773-0.04360.0250-0.06130.03740.0374-0.0022-0.04640.0106-0.0062-0.0130-0.00180.04230.0146-0.06310.03090.0845-0.08240.02270.0094-0.0032-0.01610.0068>> field.rho(130)ans =6.7659e-07>> field.rho(129)ans =-1.6066e-06>> field.rho(prm.x2(129)) Unrecognized method, property, or field 'x2' for class 'Parameters'.>> field.rho(prm.X2(129)) Index exceeds the number of array elements (128).>> field.rho(prm.X2(128))ans =-1.6066e-06>> field.rho(prm.X2)ans =1.0e-04 *0.0068-0.0272-0.02320.02960.02860.0298-0.0726-0.04510.0886-0.06980.1098-0.0172-0.05940.0846-0.05150.00560.0295-0.04370.0732-0.0528-0.05210.1384-0.16000.05710.00280.0386-0.01290.0232-0.05710.0765-0.13430.01530.0398-0.05140.01850.0307-0.05800.03200.01500.00780.0365-0.00910.0169-0.07020.0768-0.0393-0.0055-0.0171-0.0179-0.00300.0138-0.04120.08240.0326-0.0084-0.0116-0.08610.04940.00480.0140-0.04960.0690-0.0408-0.0293-0.04690.1150-0.0627-0.0000-0.04810.07430.0403-0.02000.01640.0429-0.10610.00340.0407-0.08320.06320.0992-0.0657-0.0033-0.03440.0497-0.06510.05640.0068-0.07260.0799-0.02490.0238-0.06030.01910.0342-0.08870.0843-0.05530.0564-0.02010.01760.0327-0.13490.0326-0.00020.0892-0.09870.0773-0.04360.0250-0.06130.03740.0374-0.0022-0.04640.0106-0.0062-0.0130-0.00180.04230.0146-0.06310.03090.0845-0.08240.02270.0094-0.0032-0.0161>> field.rho(2:129)ans =1.0e-04 *0.0068-0.0272-0.02320.02960.02860.0298-0.0726-0.04510.0886-0.06980.1098-0.0172-0.05940.0846-0.05150.00560.0295-0.04370.0732-0.0528-0.05210.1384-0.16000.05710.00280.0386-0.01290.0232-0.05710.0765-0.13430.01530.0398-0.05140.01850.0307-0.05800.03200.01500.00780.0365-0.00910.0169-0.07020.0768-0.0393-0.0055-0.0171-0.0179-0.00300.0138-0.04120.08240.0326-0.0084-0.0116-0.08610.04940.00480.0140-0.04960.0690-0.0408-0.0293-0.04690.1150-0.0627-0.0000-0.04810.07430.0403-0.02000.01640.0429-0.10610.00340.0407-0.08320.06320.0992-0.0657-0.0033-0.03440.0497-0.06510.05640.0068-0.07260.0799-0.02490.0238-0.06030.01910.0342-0.08870.0843-0.05530.0564-0.02010.01760.0327-0.13490.0326-0.00020.0892-0.09870.0773-0.04360.0250-0.06130.03740.0374-0.0022-0.04640.0106-0.0062-0.0130-0.00180.04230.0146-0.06310.03090.0845-0.08240.02270.0094-0.0032-0.0161>> field.rho(2:130)ans =1.0e-04 *0.0068-0.0272-0.02320.02960.02860.0298-0.0726-0.04510.0886-0.06980.1098-0.0172-0.05940.0846-0.05150.00560.0295-0.04370.0732-0.0528-0.05210.1384-0.16000.05710.00280.0386-0.01290.0232-0.05710.0765-0.13430.01530.0398-0.05140.01850.0307-0.05800.03200.01500.00780.0365-0.00910.0169-0.07020.0768-0.0393-0.0055-0.0171-0.0179-0.00300.0138-0.04120.08240.0326-0.0084-0.0116-0.08610.04940.00480.0140-0.04960.0690-0.0408-0.0293-0.04690.1150-0.0627-0.0000-0.04810.07430.0403-0.02000.01640.0429-0.10610.00340.0407-0.08320.06320.0992-0.0657-0.0033-0.03440.0497-0.06510.05640.0068-0.07260.0799-0.02490.0238-0.06030.01910.0342-0.08870.0843-0.05530.0564-0.02010.01760.0327-0.13490.0326-0.00020.0892-0.09870.0773-0.04360.0250-0.06130.03740.0374-0.0022-0.04640.0106-0.0062-0.0130-0.00180.04230.0146-0.06310.03090.0845-0.08240.02270.0094-0.0032-0.01610.0068>> field.rho(2 3 129)field.rho(2 3 129)↑ Error: Invalid expression. Check for missing multiplication operator, missing or unbalanced delimiters, or other syntax error. To construct matrices, use brackets instead of parentheses.>> field.rho(2,3, 129) Index in position 2 exceeds array bounds (must not exceed 1).>> field.rho(2,3) Index in position 2 exceeds array bounds (must not exceed 1).>> field.rho(2:3)ans =1.0e-05 *0.0677-0.2716>> field.rho([2 3 129])ans =1.0e-05 *0.0677-0.2716-0.1607>> poisson(field, prm) >> poisson(field, prm)f = Field with properties:ex: [130×1 double]ey: [130×1 double]ez: [130×1 double]by: [130×1 double]bz: [130×1 double]rho: [130×1 double]ajx: [130×1 double]ajy: [130×1 double]ajz: [130×1 double]>> poisson(field, prm)ex0 =-3.5999e-21f = Field with properties:ex: [130×1 double]ey: [130×1 double]ez: [130×1 double]by: [130×1 double]bz: [130×1 double]rho: [130×1 double]ajx: [130×1 double]ajy: [130×1 double]ajz: [130×1 double]>> poisson(field, prm)ex0 =-3.6264e-21>>

總結(jié)

以上是生活随笔為你收集整理的[MATLAB调试笔记]Possion方程求网格点电场强度的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。