wjp314 2021-09-01 16:20 采纳率: 33.3%
浏览 73
已结题

matlab中pdepe帮助文件里的矩阵问题

当pdepe函数中c矩阵不是对角矩阵时,应该怎样编写?我尝试了一下,会报错。详情如下。
帮助文档里的示例是这个样子的。

img

img
我改成非对角矩阵,发现无法计算了

img

  • 写回答

1条回答 默认 最新

  • joel_1993 2021-09-01 17:17
    关注

    人家都规定了,仅对角,如果是耦合的话,matlab求解器根本没法算除非你解耦计算。
    你可以左右同时左乘这个矩阵的逆,然后化成新的pde

    function [c,f,s] = pdefun(x,t,u,dudx) % Equation to solve
    global cmat %这里的cmat就是你的c矩阵,global是为了方便主函数调用
    c = [1; 1];
    f = cmat\([0.024; 0.17] .* dudx); % 相当于是左乘逆
    y = u(1) - u(2);
    F = exp(5.73*y)-exp(-11.47*y);
    s = cmat\[-F; F];
    end
    

    下面是主程序

    global cmat
    cmat = [1,0; 0.5,1];%这里是你需要的c的矩阵形式
    m = 0;
    x = [0 0.005 0.01 0.05 0.1 0.2 0.5 0.7 0.9 0.95 0.99 0.995 1];
    t = [0 0.005 0.01 0.05 0.1 0.5 1 1.5 2];
    sol = pdepe(m,@pdefun,@pdeic,@pdebc,x,t);
    u1 = sol(:,:,1);
    u2 = sol(:,:,2);
    surf(x,t,u1)
    title('u_1(x,t)')
    xlabel('Distance x')
    ylabel('Time t')
    

    其他的函数程序保持不变即可

    本回答被题主选为最佳回答 , 对您是否有帮助呢?
    评论

报告相同问题?

问题事件

  • 系统已结题 9月9日
  • 已采纳回答 9月1日
  • 创建了问题 9月1日

悬赏问题

  • ¥15 学不会递归,理解不了汉诺塔参数变化
  • ¥15 基于图神经网络的COVID-19药物筛选研究
  • ¥30 软件自定义无线电该怎样使用
  • ¥15 R语言mediation包做中介分析,直接效应和间接效应都很小,为什么?
  • ¥15 Jenkins+k8s部署slave节点offline
  • ¥15 如何实现从tello无人机上获取实时传输的视频流,然后将获取的视频通过yolov5进行检测
  • ¥15 WPF使用Canvas绘制矢量图问题
  • ¥15 用三极管设计一个单管共射放大电路
  • ¥15 孟德尔随机化r语言运行问题
  • ¥15 pyinstaller编译的时候出现No module named 'imp'