cyberck的博客

简单体会数值求解一维扩散方程

cyberck — Wed, 25 Aug 2021 21:12:43 +0800

借用这个题目然后用比较简单的代码来求解 %空间0到1，切了5份，时间0到1，切了十份。 clc,clear v = 1; dx = 0.2; x = 0:dx:1; dt = 0.1; nx = 6; nt = 11; t = 0:dt:1; b = v/(2*dx*dx); c = b; a = 1/dt+b+c; Uold = exp(x); Unew = ones(1,nx); Unew(1) = exp(t(2)); Unew(nx) = exp(1+t(2)); %其实最关键的就是得

Matlab:Crank Nicolson方法求解线性抛物方程

cyberck — Wed, 25 Aug 2021 20:50:59 +0800

tic; clear clc M=[10,20,40,80,160,320,640];%x的步数 K=M; %时间t的步数 for p=1:length(M) hx=1/M(p); ht=1/K(p); r=ht/hx^2; %网格比 x=0:hx:1; t=0:ht:1; numerical=zeros(M(p)+1,K(p)+1); numerical(:,1)=exp(x); %初始值 numerical(1,:)=exp(t); %边值 numerical(M(p)+1,:)=exp(t+1); .

【转载】crank-nicolson求解一维扩散方程

cyberck — Wed, 25 Aug 2021 20:43:56 +0800

OpenFOAM中的PISO算法

cyberck — Mon, 28 Dec 2020 16:54:09 +0800

流体力学的控制方程是耦合方程组，形式上体现为连续方程和运动方程的耦合，变量上体现为速度和压强的耦合。在数值计算过程中，同时求解速度和压强的方法称为“同时算法”，这通常比较困难繁琐，为此“分离算法”得到了广泛的应用。分离算法的思路是建立速度和压强之间的迭代关系，在迭代中去寻求稳定解。有一点促使了分离算法的应用：连续方程和运动方程在离散后均得到速度-压强关系。这使得在迭代过程中只需要在两个方程间进行切换，而无需想办法构造一个不动点迭代。 Jasak的论文中对OpenFOAM所采用的PISO算法给出了描述，但是

Ubuntu终端常用的快捷键

cyberck — Fri, 25 Dec 2020 22:08:24 +0800

Ubuntu中的许多操作在终端（Terminal）中十分的快捷，记住一些快捷键的操作更得心应手。在Ubuntu中打开终端的快捷键是Ctrl+Alt+T。其他的一些常用的快捷键如下：快捷键功能 Tab 自动补全 Ctrl+a 光标移动到开始位置 Ctrl+e 光标移动到最末尾 Ctrl+k 删除此处至末尾的所有内容 Ctrl+u 删除此处至开始的所有内容 Ctrl+d 删除当前字符 Ctrl+h 删除当前字符前一个字符 Ctrl+w 删除此处到左边的单词 Ctrl+y

matlab中的meshgrid

cyberck — Fri, 25 Dec 2020 19:37:08 +0800

x=[-1:0.5:1]; y=[-1:0.5:1]'; [X,Y]=meshgrid(x,y); 生成的Y是 Y = -1.0000 -1.0000 -1.0000 -1.0000 -1.0000 -0.5000 -0.5000 -0.5000 -0.5000 -0.5000 0 0 0 0 0 0.5000 0.5000 0.5000 0.5000

matlab中gradient函数的参数h

cyberck — Fri, 25 Dec 2020 19:28:49 +0800

help文件中关于gradient： [___] = gradient(F,h) 使用 h 作为每个方向上的点之间的均匀间距。您可以指定上述语法中的任何输出参数。 [___] = gradient(F,hx,hy,…,hN) 为 F 的每个维度上的间距指定 N 个间距参数。其中h的意义一开始并不明白，后来测试后才知道其实就是生成x与y的一维数组的时候的间距，有的时候是0.1 有的时候是1，取决于计算的需要比如 x=[-1:0.1:1]; y=[-1:0.1:1]'; [X,Y]=meshgrid(x,y

openfoam3.0.1编译安装

cyberck — Wed, 23 Dec 2020 15:39:29 +0800

ubuntu14.04下编译安装openfoam3.0.1

【转载】[学习笔记]共轭转置矩阵与伴随矩阵都用A*表示合理吗？

cyberck — Tue, 18 Feb 2020 15:57:32 +0800

最近在复习线代的时候被矩阵的厄米共轭弄得有点晕找到一篇文章感觉不错以下内容转载自知乎裴以鹏链接如下： https://zhuanlan.zhihu.com/p/87330558 这里是引用 [学习笔记]共轭转置矩阵与伴随矩阵都用A*表示合理吗？量子力学中，通常用上标一个dagger（匕首）来表示 ...

7-4 List Leaves (25分)

cyberck — Sun, 22 Dec 2019 02:17:13 +0800

Given a tree, you are supposed to list all the leaves in the order of top down, and left to right. Input Specification: Each input file contains one test case. For each case, the first line gives a po...

03-树1 树的同构 (25分)

cyberck — Sat, 21 Dec 2019 15:50:08 +0800

给定两棵树T1和T2。如果T1可以通过若干次左右孩子互换就变成T2，则我们称两棵树是“同构”的。例如图1给出的两棵树就是同构的，因为我们把其中一棵树的结点A、B、G的左右孩子互换后，就得到另外一棵树。而图2就不是同构的。图1 图2 现给定两棵树，请你判断它们是否是同构的。输入格式: 输入给出2棵二叉树树的信息。对于每棵树，首先在一行中给出一个非负整数N (≤10)，即该树的结点数（此时假设...

02-线性结构1 两个有序链表序列的合并 (15分)

cyberck — Fri, 20 Dec 2019 21:04:47 +0800

本题要求实现一个函数，将两个链表表示的递增整数序列合并为一个非递减的整数序列。函数接口定义： List Merge( List L1, List L2 ); 其中List结构定义如下： typedef struct Node *PtrToNode; struct Node { ElementType Data; /* 存储结点数据 */ PtrToNode Next; /*...

02-线性结构2 一元多项式的乘法与加法运算

cyberck — Fri, 20 Dec 2019 21:01:19 +0800

设计函数分别求两个一元多项式的乘积与和。输入格式: 输入分2行，每行分别先给出多项式非零项的个数，再以指数递降方式输入一个多项式非零项系数和指数（绝对值均为不超过1000的整数）。数字间以空格分隔。输出格式: 输出分2行，分别以指数递降方式输出乘积多项式以及和多项式非零项的系数和指数。数字间以空格分隔，但结尾不能有多余空格。零多项式应输出0 0。输入样例: 4 3 4 -5 2 6 1 ...

C语言程序设计进阶-第4周编程题（链表）

cyberck — Mon, 16 Dec 2019 18:34:42 +0800

1 逆序输出（10分）题目内容：你的程序会读入一系列的正整数，预先不知道正整数的数量，一旦读到-1，就表示输入结束。然后，按照和输入相反的顺序输出所读到的数字，不包括最后标识结束的-1。输入格式: 一系列正整数，输入-1表示结束，-1不是输入的数据的一部分。输出格式：按照与输入相反的顺序输出所有的整数，每个整数后面跟一个空格以与后面的整数区分，最后的整数后面也有空格。输入样例： 1 2...

C语言程序设计进阶-第1周编程题（指针与字符串）

cyberck — Sat, 14 Dec 2019 11:48:03 +0800

1 字符串比对（10分）题目内容：题目说起来很简单，你会读到两个字符串，每个字符串占据一行，每个字符串的长度均小于10000字符，而且第一个字符串的长度小于第二个字符串的。你的程序要找出第一个字符串在第二个字符串中出现的位置，输出这些位置，如果找不到，则输出-1。注意，第一个字符的位置是0。注意，第一个字符串在第二个字符串中的位置可能不止一处。注意，字符串中可能含有空格。注意，两个字符...

第七周作业

cyberck — Thu, 28 Nov 2019 21:17:35 +0800

#include int main() { int a[101]={0}; int j=0; int i=0; int d=0; do{ scanf("%d",&j); scanf("%d",&a[j]); }while(j!=0); ...

第六周作业

cyberck — Thu, 28 Nov 2019 18:15:47 +0800

#include int main() { int a,b; scanf("%d/%d",&a,&b); int i=0; int j=a; int count=0; printf("0."); do{ i=j*10/b; ...