- 博客(79)
- 收藏
- 关注
RSS订阅原创 传热学一些“数”和意义
Nu与Bi的表达式相同，但是意义是无量纲的h。它们表达式里的长度取值不同，比如同样一个平板，Bi的L是厚度，Nu是长度，因为Bi面向固体，λ为固体的，Nu的λ是气体的。管子的特征长度是直径，板子如果单面受热特征长度就是一个厚度，如果双面受热特征长度就是半个厚度。速度边界层厚度/温度边界层厚度 （粘性扩散能力/热扩散能力）物体单位面积上的导热热阻/单位表面积上的对流换热热阻。
2024-10-08 16:00:00
200
原创 关于直肋的效率
这个的分子的导数小于0，分子的零点明显为0，因此在x>0时该导数为负，因此随mH的减小，效率提升。可以如此理解，由效率表达式上面的分子在x趋近于∞时为1，可以看做一个反比例函数。肋片越矮胖，肋片各个部分的温度越趋近于肋根温度，效率越高。对于截面为圆柱形的，效率随半径的增大而增大。假如固定其他条件，当λ变大，mH减小，对。效率随材料的导热系数的增大而增大。效率随对流换热系数的增大而减小。效率随肋片的增高而减小。进行求导，得到导数为。
2024-10-08 15:50:24
675
原创 k=1.4,R=287.06. 总温和临界声速的关系
T*c_cr100183.003643101183.916385102103104105106107108109110191.93584111112113114115116117118119120121122123124125126127128129130131132133134135136137138139215.758113140141217.304786142143144145146147。
2024-10-08 08:46:18
117
原创 关于复变函数解析与可导
一般情况下，有lm(z)或|z|存在的函数不可导，但是也有例外，比如。一般看f(z)在定义域D内解析只看是否满足柯西-黎曼方程，即。两个式子同时成立，但是还要满足u和v在D内可微。
2024-10-06 18:24:56
992
原创 关于abaqus里一些问题的记录
为了帮助您直观地了解偏差系数的影响，Abaqus/CAE 会显示它将围绕与您输入的设置相对应的圆圈创建的元素数。指定最小尺寸控件有助于防止 Abaqus/CAE 在高曲率区域创建不必要的精细网格（如果这些区域对分析意图并不重要）。最后输出的结果没有力，可能是因为材料的参数给的不对，忘记了压强的单位是MPa，尤其注意密度单位t/mm^2。在进行布种时，会遇到最大偏离因子和最小尺寸因子，在帮助文档里，是这么解释。这里有一个恢复模块，不是回弹，是重新提交之前被中断的explicit任务。输入绝对最小元素大小。
2024-10-04 08:30:07
305
原创 传热学测温套管相关
公式的推导过程见传热学的书籍，假设套管长为H，外径为d，管厚δ。套管伸入大管子，大管子壁上温度为t0，温度计示数为tw，则外温是。，如果给定测量误差α和来流温度tf，则θ=αtf。对于一般的肋片，如果顶端绝热，则顶端的过余温度为。，如果已知m，可以由此解出H，代码如下。
2024-09-30 21:20:13
666
原创 关于特征方程的想法
解常微分方程所用到的特征方程或许和拉氏变换有关，比如解。，这分母正是特征方程，可以解出通解为。，这和由特征方程解出的一样，其中。
2024-09-18 14:24:13
954
原创 matlab实现简单的保角变换
用虚数的思想可以在虚坐标系内绘制圆，同样可以用虚数的想法将这个圆进行变换。用MATLAB绘制一个坐标在(1,1)的圆代码如下。
2024-09-04 14:21:10
400
原创 古典显示格式解一偏微分方程并绘制结果的彩色图
最后画图时如果用mesh(x,t,T)，会显示一个三维彩色网格图如下。如果想要实心的格子，就用surf(x,t,T)，效果如下。通过转动也可以显示出x-T或t-T二维图。以上公式的Latex代码。MATLAB代码如下。
2024-09-02 17:00:35
240
原创 离散的点进行积分
这里的代码表明，拟合后的结果可以直接转化为句柄，但并不是标准的形式，因此不能用integral直接积分。2，可以先进行多项式拟合，再积分。1，可以直接用梯形公式。
2024-09-01 09:15:47
272
原创 穿越沙漠小游戏的Python代码
本游戏以命令行的形式进行，玩家可以通过改变代码的方式改相关的参数，游戏背景请参考。这里假设在一条线上移动，天气随机决定。
2024-08-28 10:35:45
143
原创 logistic方程求解和模拟
其中N(t)为t时刻种群的数量，N为环境最大容纳的的种群数。k的取值通常与N有关，大约在1/N和4/N之间。logistic规律可以用以下式子表示。假设一些参数后的代码如下。
2024-08-25 13:53:23
246
原创 采用不高于3次的勒让德多项式拟合原函数
利用勒让德多项式进行拟合的区域是[-1,1]，如果不是这个区域，比如是[a,b]，利用。integrate为求积分的函数，见下。如果用MATLAB，则可以使用函数句柄。转化到[-1,1]。参考以下例题计算系数。
2024-08-24 14:41:11
145
原创 一个可以搜索页面内超链接的HTML页面
页面效果如上，含有标题，搜索框和一些超链接，在搜索框输入一些文字可以在下方显示含有这些文字的超链接。这里的搜索不仅可以中文，还可以是英文，数字。
2024-08-24 12:32:19
335
原创 用C#写一个随机音乐播放器
链接：https://pan.baidu.com/s/105w5fij6kkrUoWcjMuiBWA?这里的文件夹地址是直接Ctrl+Shift+C复制过来的，注意要去掉引号。实际是随机打开一个文件，鉴于VS内置的库不能播放除了wav以外的音频。form1中namespce里的代码如下。
2024-08-19 13:34:15
635
原创 两种用MATLAB绘制色块的方法
这种方式本质上是创造一个图片，因此显示时没有坐标轴，不过由于仅仅填充像素，当尺寸设置过小，图像也过小。绘制色块首先可以想到用填充像素的方式。另外可以想到填充用坐标线围成的网格。
2024-08-19 12:58:54
366
原创 y=λsin(πx)分岔的研究
加密步长，取r区间为[0.4,0.9]，绘制feigenbaum图。用如下代码对收敛的λ的值进行探究，这里的r代表λ。使用如下的迭代格式，λ为可变的参数。可以观察到从0.9开始就开始分岔。
2024-07-25 22:21:27
224
原创 简单的差分格式解一个一元二阶常微分方程
将a到b的区域平均分成m份，解以上的式子可以看做解一个三对角方程组。这里也包含了追赶法求三对角矩阵的函数。对于一个如下的一元二阶常微分方程。其中二阶导数可以换成差分格式。之后展开成的矩阵形式，取。这里是左右两边同时乘以。
2024-07-06 11:46:41
749
1
原创 C++创建一个可以在不同源文件之间访问的二维数组
这里定义了数组的大小，定义二维数组时要在前面加extern。首先需要创建一个源文件（如Moudle.h），代码如下。之后，可以在一个源文件里给u赋值。注意这里u的声明要放函数外面。在主文件里，代码如下。
2024-07-05 13:53:05
116
原创 MATLAB将两个折线图画在一个图里
输入行数和列数，点击开始填入数据，其中第一列为x值，后面几列，每一列都是y坐标值，填好后点击画在同一张图里即可。点击置零就把所有数变成0，另外也可以选择节点样式。.mlapp格式的文件如下。
2024-07-02 16:12:57
600
1
原创 关于一维，二维正态分布的绘制
完全正相关或完全负相关，即一个变量的增加（减少）总是伴随着另一个变量按确定比例的增加（减少）。，可以在保持各自分布形态不变的情况下，描述两个随机变量之间不同的依赖关系。二维正态分布的概率密度函数（PDF）可以用两个一维正态分布的随机变量。是独立的，其联合分布可以分解为各自边际分布的乘积；和两个随机变量标准差的乘积的比值，即。之间的线性关系强度和方向，由相关系数。对于二维正态分布，其联合概率密度函数。若要算某一点的概率密度，代码如下。，以及它们之间的相关系数。绘制二维正态分布代码如下。本文采用AI辅助创作。
2024-07-02 08:04:53
1014
原创 强度理论的相当应力在坐标图里的表示
以下的图均只考虑二维，且没有切应力，有切应力的情况可以转换成只有两个正应力的情况，另外，这里画出的图均认为相当应力就是许用应力。两种强度理论的图均明显有一个转角，这是因为主应力的转变，当。从图形上看，在极坐标系下的和在直角坐标系下的相同。这里绘制畸变能理论的图，公式为。第一，二强度理论用于脆性材料。绘制最大拉应变理论的代码如下。第三，四强度理论用于脆性材料。
2024-07-01 13:00:04
181
原创 学习数值方法解常微分方程的笔记
比如y''' + 2y'' + 2y' + y = 80，化为y1' = y2, y2' = y3, y3' = 80 - 2y3 - 2y2 - y1，这个80 - 2y3 - 2y2 - y1相当于一阶里的f(t,y)工程上在满足精度要求的前提下，多用4阶龙格库塔法解常微分方程。取步长0.001，步数10，C语言代码如下。对于多阶常微分方程，可以逐级化为一阶的。与精确值的误差几乎没有。
2024-04-27 11:30:47
276
1
原创 单片机小项目——直流电机+按键
不过在未按下按键时，电动机仍然有一定的运转，而且在按下前15个按键的瞬间电动机仍要加速转动，猜想是矩阵键盘的扫描与电机共用一个IO口，互相影响了电平。使用4*4键盘控制，第一个到第15按键分别对应数码管显示1到F，计划在按下第16的按键时电动机运转数码管显示的秒数，代码如下。利用普中单片机的代码实现在按下第k个独立按键时，直流电机运作k秒。
2024-04-20 22:32:46
616
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人