飞蓝飘雪 2021-10-17 16:23 采纳率: 100%
浏览 49
已结题

微分方程编程MATLAB

img

你好,请问这个微分方程怎么用MATLAB求解啊?谢谢,希望给个解答或者思路。真心感谢

  • 写回答

1条回答 默认 最新

  • joel_1993 2021-10-17 16:54
    关注

    你好,这个边值问题的方程建议用中心差分法来求解,代码供参考

    L = 1;%L大小
    nx = 100;%100段
    x = linspace(0,L,nx+1);
    dx = L/nx;
    afun = @(x) (1+x.^2); % a随x的分布函数
    ffun = @(x) (3-x).^2; % f随x的分布函数
    % d^2u/dx^2 = -f./a
    % (u(i-1)+u(i+1)-2*u(i))/dx^2=d^2u/dx^2|x(i) %离散方法
    n = length(x);% 布点总个数
    a = afun(x);%求x对应的a
    f = ffun(x);%求x对应的f
    K = diag(-2*ones(n-2,1)) + diag(ones(n-3,1),1) + diag(ones(n-3,1),-1);
    F = - f(2:end-1)'./ a(2:end-1)'*dx^2;
    u0 = 0; % x=0处边界条件
    P = 1; uL = P; % x=L处边界条件
    F(1) = F(1) - u0;% 施加边界条件
    F(end) = F(end) - uL;% 施加边界条件
    u = K\F;%求解u(除两个端点之外的u)
    utot = [u0; u; uL];%囊括两个端点之后的u
    plot(x,utot)% 画出u随x变化图
    xlabel('x'); ylabel('u')
    set(gca, 'fontsize',16)
    

    效果图u-x

    img

    有帮助还请给个采纳支持一下答主冲榜哦,万分感谢

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

报告相同问题?

问题事件

  • 系统已结题 10月25日
  • 已采纳回答 10月17日
  • 创建了问题 10月17日

悬赏问题

  • ¥15 MATLAB动图问题
  • ¥15 【提问】基于Invest的水源涵养
  • ¥20 微信网友居然可以通过vx号找到我绑的手机号
  • ¥15 寻一个支付宝扫码远程授权登录的软件助手app
  • ¥15 解riccati方程组
  • ¥15 display:none;样式在嵌套结构中的已设置了display样式的元素上不起作用?
  • ¥15 使用rabbitMQ 消息队列作为url源进行多线程爬取时,总有几个url没有处理的问题。
  • ¥15 Ubuntu在安装序列比对软件STAR时出现报错如何解决
  • ¥50 树莓派安卓APK系统签名
  • ¥65 汇编语言除法溢出问题