查看: 92|回复: 0|关注: 0

[未答复] 如何求叠加和积分的混合解

[复制链接]

新手

5 麦片

财富积分


050


1

主题

1

帖子

0

最佳答案
发表于 2019-10-10 16:27:35 | 显示全部楼层 |阅读模式
本帖最后由 demonseasons 于 2019-10-10 16:32 编辑

目前我编的代码如下

  1. for m1=1:Com
  2.     for j=1:size(tDe,2)  %时间循环
  3.         f=log(2)/tDe(j);
  4.         for i=1:N
  5.             PwD(m1,j)=0;
  6.             z=i*f;
  7. %            u=z;
  8.               syms arph n   
  9.               beta=n*pi;
  10.               kesai0=sqrt(z*hD^2*LD^2+beta*LD^2);
  11.               kesai1=besselk(0,sqrt((xD-arph)^2)*kesai0);
  12.               int2=int(kesai1,arph,-1,1);%积分部分
  13.               int2=vpa(int2);
  14.               int3=int2*cos(beta*zD)*cos(beta*zwD);
  15.               int1=symsum(int3,n,1,inf);
  16.                 pD(m1,j)=0.5/z*int2+2*int1;
  17.                 PwD(m1,j)=PwD(m1,j)+V(i)*pD(m1,j);            %报错位置
  18.         end
  19.             PwD(m1,j)=PwD(m1,j)*f;
  20.     end
  21. end   
复制代码

错误展示:
  1. 错误使用 symengine
  2. Unable to convert expression into double array.

  3. 出错 sym/double (line 661)
  4.         Xstr = mupadmex('symobj::double', S.s, 0);

  5. 出错 horizontal (line 71)
  6.                 PwD(m1,j)=PwD(m1,j)+V(i)*pD(m1,j);
复制代码


本人编程小白,在网上查到可能是由于数值解和解析解混淆了,但不知道怎么修改,希望得到各位菊苣的指点!


公式

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

本版积分规则

关闭

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

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