renbaifen的博客

fvm_1d_swe理论主要是参考：【教学向】从零开始的流体力学（一维浅水方程+有限体积法） - 知乎 (zhihu.com)而且我话的不是动图，而且耗散项系数为0时，也太震荡了。function fvm_test_1dclear all ;clc;L = 1;div = 100;%space numberdx = L/div;x0 = 0;x1 = 1;y00 = 0;y1 = 0.25;t0 = 0.0;Final_time = 5;g = 1.0;x = (1:

clear all;clc;x0 = 1;b = 1; y0 = 1;z0 = 1; t0 = 0;h = 0.1;df = @diff;[k,X,Y] = Adams4t(df,t0,b,x0,y0,z0,h);% temp = exact(X);% plot(X,Y,'r',X,temp,'*');plot(X,Y,'*');legend('x','y','z');% function result = exact(x)% result = 3./(1+x.^3.

算法的原理就是类似于LU分解，然后求解线性方程组。具体可以参考https://www.cnblogs.com/xpvincent/archive/2013/01/25/2877411.html方法一a = [0,1,1,1]';b = [4,4,4,4]';c = [1,1,1]';x = [1,1,1,0]';m = zeros(3,1);n = 4;c(1) = c(1)/b(1);x(1) = x(1)/b(1);for i = 2:n-1 c(i) = c(i)/

这个博客主要写了高等教育出版社《数值方法引论》 第十一章常微分方程的数值解法的MATLAB 代码，是一个课后题page257的参考代码，问题的条件如下：[0,1]区间,步长h = 0.1;以下所有的代码段已经打包上传到我的代码资源中了，有需要的可以下载。https://download.youkuaiyun.com/download/renbaifen/15810690?spm=1001.2014.3001.5501线性多步，内插function y = multistep_in() %线性多

以下是matlab绘制箱线图的代码   n=10;%n是离散样本点的个数，，，这里样本点的个数可以设的大一些k=1.5;%k是阈值%a=rand(1,n)*n;%产生随机数a=[2.0246,7.6122,6.9196,5.7926,5.6371,3.6518,3.1846,1.6223,0.808,8.6951];Q1=2.3146;%上四分位线,,上下四分位数可以用e...

//阻尼牛顿法function [k,x,val] = dampnm(fun,gfun,hess,x0,epsilon)maxk = 5000;beta = 0.5;sigma = 0.4;k = 0;while(k<maxk)    gk = feval(gfun,x0);    Gk = feval(Hesse,x0);    dk = -Gk\gk;%解方程组，计算搜...

function [k,x,val]=bfgs(fun,gfun,x0,varargin)%输入：fun，gfun 分别是目标函数与目标函数的梯度，x1%是初始目标点，varargin是输入的可变参数变量%输出：k迭代次数，x，val近似最优点和最优值%N=500;N=1000;epsilon = 1.e-5;%给定容许误差beta = 0.55;sigma = 0.4;n=le...

fid=fopen('data.txt','w');for i = 1:size(Echo_h1,1) for j = 1:size(Echo_h1,2) fprintf(fid,'%f\t',Echo_h1(i,j)); endendfclose(fid);