binbinyouli0206 2024-09-05 14:15
浏览 3

PINN解决频域helmholtz问题

PINN解决频域helmholtz方程,其中一个边界条件是dp/dx-iZnCUnK*ρ=0,Zn=1,ρ=1.259,C=340,K=10,Un=10,边界条件包括虚数i,这该怎么处理?我写的代码这左边界条件损失函数一直维持在42801附近?该怎么解决?

def equation2(net):  # 左边界
    x, y, cond = sample2(size1)
    preal = net(torch.cat([x, y], dim=1))[:, 0].unsqueeze(-1)
    pimage = net(torch.cat([x, y], dim=1))[:, 1].unsqueeze(-1)

    preal_x = gradients(preal, x, 1)
    pimage_x = gradients(pimage, x, 1)
    p_x = preal_x + 1j * pimage_x
    ke = torch.tensor(10, dtype=torch.float32)
    c = torch.tensor(340, dtype=torch.float32)
    rou = torch.tensor(1.259, dtype=torch.float32)
    z = torch.tensor(1, dtype=torch.float32)
    uba = torch.tensor(10, dtype=torch.float32)
    left_fun = p_x - 1j * ke * z * rou * c * uba

    loss_left = torch.mean(torch.abs(left_fun))
    
    return loss_left

  • 写回答

0条回答 默认 最新

    报告相同问题?

    问题事件

    • 修改了问题 9月5日
    • 创建了问题 9月5日

    悬赏问题

    • ¥60 如何批量获取json的url
    • ¥15 对法兰连接元件所承受的表面载荷等效转化为法兰开孔接触面上的等效表面载荷?
    • ¥15 comsol仿真压阻传感器
    • ¥15 Python线性规划函数optimize.linprog求解为整数
    • ¥15 llama3中文版微调
    • ¥15 pg数据库导入数据序列重复
    • ¥15 三分类机器学习模型可视化分析
    • ¥15 本地测试网站127.0.0.1 已拒绝连接,如何解决?(标签-ubuntu)
    • ¥50 Qt在release捕获异常并跟踪堆栈(有Demo,跑一下环境再回答)
    • ¥30 python,LLM 文本提炼