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

[已答复] matlab数值积分

[复制链接]

新手

10 麦片

财富积分


050


1

主题

3

帖子

0

最佳答案
发表于 2019-8-12 14:09:00 | 显示全部楼层 |阅读模式
本帖最后由 溜溜梅 于 2019-8-12 17:53 编辑

求高手帮我看看这个程序怎么计算出表达式的积分,可能是因为结果太小而只显示0,想用符号变量来做,但是integral2函数又会提示错误,错误是说调用某个函数的时候要求的是数值,但是我用的符号。求大神指教,万分感谢f.m中的代码
  1. function z = f(a,x )
  2. n=80;   
  3. p=8;
  4. K=50;
  5. T=2309;
  6. beta=711;
  7. theta=3.6288e+06;
  8. w=[126;126;96;96;63;63;30;30];
  9. TIJ=[686,810,1658,1790,2185,0,0
  10.     801,2202,2212,2217,0,0,0
  11.     351,1026,1840,2150,0,0,0
  12.     1097,1584,1755,2032,0,0,0
  13.     233,359,701,1867,2110,0,0
  14.     1842,0,0,0,0,0,0
  15.     0,0,0,0,0,0,0
  16.     310,530,918,1096,1582,0,0
  17.     106,118,1713,2056,0,0,0
  18.     102,525,952,1209,1267,1716,2125
  19.     1962,2053,2119,2174,0,0,0
  20.     819,1479,0,0,0,0,0
  21.     1029,1074,1112,1147,1865,1954,0
  22.     0,0,0,0,0,0,0
  23.     1795,1862,1963,0,0,0,0
  24.     1158,1744,0,0,0,0,0
  25.     1665,0,0,0,0,0,0
  26.     1465,0,0,0,0,0,0
  27.     1269,0,0,0,0,0,0
  28.     59,1164,1500,0,0,0,0
  29.     0,0,0,0,0,0,0
  30.     0,0,0,0,0,0,0
  31.     0,0,0,0,0,0,0
  32.     0,0,0,0,0,0,0
  33.     0,0,0,0,0,0,0
  34.     177,334,592,639,820,1792,0
  35.     1461,0,0,0,0,0,0
  36.     914,0,0,0,0,0,0
  37.     0,0,0,0,0,0,0
  38.     0,0,0,0,0,0,0
  39.     1026,0,0,0,0,0,0
  40.     0,0,0,0,0,0,0
  41.     1170,0,0,0,0,0,0
  42.     0,0,0,0,0,0,0
  43.     779,0,0,0,0,0,0
  44.     0,0,0,0,0,0,0
  45.     0,0,0,0,0,0,0
  46.     0,0,0,0,0,0,0
  47.     0,0,0,0,0,0,0
  48.     0,0,0,0,0,0,0
  49.     458,0,0,0,0,0,0
  50.     0,0,0,0,0,0,0
  51.     569,992,0,0,0,0,0
  52.     540,0,0,0,0,0,0
  53.     0,0,0,0,0,0,0
  54.     936,0,0,0,0,0,0
  55.     532,0,0,0,0,0,0
  56.     0,0,0,0,0,0,0
  57.     484,0,0,0,0,0,0
  58.     0,0,0,0,0,0,0];
  59. TIE=[2309;2309;2279;2279;2246;2246;2213;2213;2096;2182;2174;2183;2131;1452;2174;2086
  60. 2138;1938;2058;2138;1543;1248;1603;1804;1523;2175;1573;1196;1445;1145;1235;1270;1266
  61. 1369;788;983;1262;1019;1178;728;1121;966;1130;1029;1180;944;1051;672;781;28];
  62. w=round(w.*9.45);
  63. TIJ=round(TIJ.*9.45);
  64. TIE=round(TIE.*9.45);
  65. T=round(T*9.45);
  66. S1=0;
  67. for i=1:1:p
  68.     S1=w(i).^x+S1;
  69. end
  70. S2=0;
  71. for i=1:1:K
  72.     S2=TIE(i).^x+S2;
  73. end
  74. non0=(TIJ~=0);
  75. ni=sum(non0,2); %计算ni的值
  76. S4=1;
  77. for i=1:1:K
  78.     for j=1:1:ni(i)
  79.         if ni(i)==0
  80.             break;
  81.         end
  82.         S4=((TIJ(i,j)).^(x-1)).*S4;
  83.     end
  84. end

  85. A = S2-S1+theta;
  86. B = S4;
  87. % z = a.^(n+beta).*x.^n.*exp(-(A).*a).*B;        %%最终表达式
  88. end

复制代码
main.m中的代码
  1. x1=0.9;
  2. x2=1.1;
  3. fa=@(a,x)f(a,x);
  4. integral2(fa,0,inf,x1,x2)
复制代码



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

举报

论坛优秀回答者

中级

513 麦片

财富积分


5001500


5

主题

1163

帖子

112

最佳答案
  • 关注者: 24
发表于 2019-9-18 20:12:46 | 显示全部楼层
过去太久,不知道你是不是已经解决了
回复此楼 已获打赏: 0 积分

举报

新手

10 麦片

财富积分


050


1

主题

3

帖子

0

最佳答案
 楼主| 发表于 2019-9-18 20:16:30 | 显示全部楼层
jiangxuzhen 发表于 2019-9-18 20:12
过去太久,不知道你是不是已经解决了

还没有啊,还望指教啊大神:'(:'(
回复此楼 已获打赏: 0 积分

举报

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

本版积分规则

关闭

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

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