cyberck的博客

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

借用这个题目然后用比较简单的代码来求解%空间0到1，切了5份，时间0到1，切了十份。clc,clearv = 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方法求解线性抛物方程

tic;clearclcM=[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求解一维扩散方程

转载 OpenFOAM中的PISO算法

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

转载 Ubuntu终端常用的快捷键

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

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

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编译安装

ubuntu14.04下编译安装openfoam3.0.1

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

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

原创 7-4 List Leaves (25分)

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分)

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

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

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

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

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

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

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

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

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

原创 第七周作业

#include<stdio.h>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); ...

原创 第六周作业

#include<stdio.h>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; ...