- 博客(85)
- 收藏
- 关注
原创 COMSOL求解器教程
PDE到FEM弱形式的好处:①把解空间做了一定弱化,方便数值求解。②降低要求解的PDE的阶次。有限元方法对于更一般的、包含了时间项的问题,要对时间进行离散。Jacobian矩阵代表的就是线性下降速度最快的方向。一些出错的可能性①给定的初始值不合理②边界条件设置不合理。
2024-12-25 15:56:28
740
原创 COMSOL异常处理
不收敛的原因很多,主要还是要对研究的现象有充分的了解。来解决由于参数值和初始值差异较大导致的瞬态不收敛问题。只支持从“有细节的一面”向“无细节的一面”进行扫掠。默认式子的另一侧为0。COMSOL中的表达式。
2024-12-23 14:02:02
317
原创 COSMOL中的状态变量
在的下方,点击,并勾选。在Component中的中定义状态变量,输入初始值和更新表达式。状态变量的更新指的是 用initial value作为第一个输入的值,计算第一次的结果。After step指的是 把计算出来的结果作为输入的值,结果作为第一次的结果。
2024-12-18 16:25:01
682
原创 冷热水混合器的设计与优化和事件的应用
当出水口温度低于50[degC]时需要开启水泵进行加热;当出口水温高于55[degC]时停止水泵加热。降低入口水的流速能够使得(有水泵加热的)入口水温升高,最终使整体水温升高。很多时候加热是功率加热,并没有给定恒定的温度。隐式事件1:温度低于50[degC],则开关打开。将入口水温由80[degC]更改为20[degC]隐式事件2:温度高于55[degC],则开关关闭。将功率加热的功率设置为。定义开关开启和关闭的。出口平均温度aveop1(T)在初始时刻,开关是开启的。,设置字母,并指定为。
2024-04-02 11:08:11
663
原创 COMSOL热应力仿真
在参数部分,设定体积参考温度Tref,假定在25[degC]模型无热应变。APP开发器-新建表单-选择表单模板-挑选参数-完成。假如直接知道热膨胀大小,可以直接对。对于自由膨胀,可以添加抑制刚体运动。
2024-03-09 14:51:56
2326
原创 COMSOL中使用自定义函数
下方if函数的意思是,如果u<0,则u=0;如果u>0,则u=u。基于RGB值自动生成二维函数。利用三角函数定义周期性矩形波。物理场参数的跳变情况。
2024-03-08 15:14:01
5409
原创 COMSOL传热建模
既有固体,又有流体的情况,包括固体和流体之间的传热。基于显微镜下模型积分,得到宏观模型。(2)在参与介质中的辐射。(1)固体传热(推荐)(3)薄壳传热(推荐)(1)表面对表面辐射。
2024-03-02 11:13:37
2196
原创 Chapter20 音乐
每个键的频率比相邻的键高或低半音程。每一个黑键可以理解成上一个白键的升音(sharp,记为#),或下一个白键的降音(flat,b)和弦是两个或两个以上的音符同时发声的演奏方式。在pianoex程序界面上可以选中1到12的双态开关来产生和弦,标注为0的开关总处于选中的状态。在微型键盘上用for循环演奏出两个八音度的25个琴键半音音阶。因为12个半音程构成了二倍的八音度频率范围,所以一个半音程的频率倍数为。不同的调式在不同的特征频率或特征值下振动。的幂次定义,纯律是由一系列分数定义的。八音度(octave)
2023-08-04 23:43:44
336
原创 matlab编程实践18、19
waterwave求解的区域为正方形区域,有反射的边界条件。在初始时刻,在整个区域都有h=1,u=0,v=0,这样解是静态的。然后在连续几个时间步长内,二维高斯型峰值添加到h处,用来模拟水滴滴到水面上的冲量扰动作用。从根节点开始向左移动一个链接表示一个“点”,向右表示一个“划”。,这种递归调用的结果把C表格和含有“点”和“划”的dd字符串合并。'...---...'表示SOS,‘...--...’表示SMS。栈S为单元数组,只要栈是非空的,while循环就一直进行下去。在这样的记号下,浅水方程为。
2023-08-02 23:19:39
863
原创 matlab编程实践16、17
单物种logistic模型是有求解公式的,但对捕食者与猎物模型来说,不能得出包括指数函数、三角函数或其他基本函数的解析解,这时。捕食者和猎物模型为两个方程,这设计两个竞争物种y1(t)和y2(t)的变化情况。在使用ode45求解时,@符号和@(t,y)可以定义出t和y的函数。该简单模型在人口增长的初始阶段是正确的,因为初始阶段对人口没有限制。在初始时刻小球上抛之后,地球的引力是的速度每步都按照固定的比率g减少。,当y(t)接近于μ,增长率接近于0,人口增长逐渐停止。在人口增长或衰减的最简单模型中,增长。
2023-08-01 21:17:57
1397
原创 matlab编程实践14、15
数独网格的一些运算可能改变其在图形用户界面中的显示,但不能改变其基本特性。将单选数填在网格内,如果没有单选数,则采用递归回溯法来求解。matlab提供了很多求给定常微分方程数组近似解的函数,这一常微分方程组数值解的函数包括。比如ode23算法比较二阶和三阶算法来估计计算步长,而ode45比较的是四阶和五阶算法。以上的数独是图14-2的两个矩阵的和,左侧的矩阵可以由。函数名中的数字表示所用算法的阶次,阶次和算法的复杂程度和精度有关。所选择的阶次越高,每一步计算量越大,但是所用的步长也越大。
2023-07-31 22:53:10
965
原创 MATLAB编程实践12、13
如果一个存活的网格有两个(或三个)存活的邻居,或者一个死亡的网格周围有三个存活的邻居,则它在下一步中也是存活的。滑翔机五网格群体初始结构,每一步进化过程中,都有两个网格死亡,并有两个新网格出生。的离散时间步距进化。是一个包含z0值的复平面区域,其轨迹定义为。是无限可扩展的二维矩形网格,
2023-07-30 22:55:44
416
原创 幻方问题(Magic Squares)
n×n幻方是含有1到n^2的整数数组,排列后是的每一行、每一列、正反两条主对角线上数字的和都是相同的。对于每个n>2都有很多不同的n×n幻方,函数magic(n)只产生一个特定的幻方。rot90将数组旋转 90 度sum(A)sum(A')'sum(1:9)/3rot90(A,k)end。
2023-07-27 18:56:29
1028
原创 T形积木(T puzzle)
函数负责处理鼠标释放事件,即在将T形积木拖动到其他块体附近并释放鼠标键时,程序会检测是否有其他块体的顶点靠近被释放的积木,如果有足够近的块体,则将积木自动吸附到靠近的块体上。这会导致在旋转过程中可能出现明显的跳跃,特别是在慢速旋转时,可能感觉不够平滑。(3)T形积木问题利用了复数乘以。将该复数值旋转b角的事实。查找具有特定属性的图形对象。
2023-07-26 23:32:36
602
原创 指数函数exp
对指数函数进行二项式展开,求和式的每一项都通过前一项简单计算得出。该循环在r==s时结束,即两个相连的部分浮点数相同时结束。(根据比值=1推导出e的值)年利、月利、连续利息。
2023-07-25 21:59:09
516
原创 蕨型叶分形
蕨型叶是通过一个点的反复变换产生的,假设x是一个含有两个分量的向量,可以用来表示平面内的一个点,则可以用Ax+b的形式对其进行变换。语句darkgreen=[0 2/3 0]设置深绿色的颜色变量,在matlab中用红绿蓝分量来表示颜色,这里将绿色分量设置为2/3,其余为0。在RGB表示法中,每个颜色的红、绿、蓝三个分量的取值范围为0到1,表示颜色的深浅程度。在MATLAB中,颜色可以使用多种方式来表示,其中一种常用的方式是使用。来绘制蕨型叶,h为蕨型叶图形的句柄。蕨型叶的起始点均为(0,0)。
2023-07-24 13:38:01
1167
2
原创 矩阵Matrices
返回一个主对角线元素为 1 且其他位置元素为 0 的数组。返回一个主对角线元素为 1 且其他位置元素为 0 的 2×3 数组。返回一个主对角线元素为 1 且其他位置元素为 0 的。返回一个主对角线元素为 1 且其他位置元素为 0 的。相同的数据类型、稀疏性和复/实性(实数或复数)。返回由 8 位整数组成的 5×5 单位矩阵。房子图形的矩阵变换(wiggle动态变换)您指定的输出参数的数目决定输出的维度。ndgrid-N 维空间中的矩形网格。指定对所有维度使用单一网格向量。在上述语法的基础上进一步指定。
2023-07-07 23:32:47
725
原创 日历与时钟
此 MATLAB 函数 将输入数组 t 中的 datetime 或 duration 值转换为日期序列值。删除当前图窗的所有子级,不管其句柄可见性如何。在公历中,当年份为4的整数倍,但不是100的整数倍时,会出现闰年的现象。此 MATLAB 函数 返回表示 D 中每个元素的星期几的数字。使当前图窗可见,并将其置于屏幕上所有其他图窗的前面。以400年为周期的历法下,某个月的13日正好是星期五的概率。体力周期为23天,感情周期为28天,智力周期为33天。删除指定图窗的所有子级并重置其属性。
2023-07-04 15:17:56
657
2
原创 斐波那契数列
此 MATLAB 函数 在 x 轴上使用线性刻度、在 y 轴上使用以 10 为底的对数刻度来绘制 x 和y 坐标。MATLAB 读取执行 toc函数时的内部时间,并显示最近一次不带输出调用 tic 函数以来经过的时间。一个人将一对兔子放到一个封闭的围墙内,并假设每对兔子每个月都繁殖出一对兔子,且新生兔子从第二个月开始有繁殖能力,那么一年以后这个封闭的围墙内有多少只兔子。本月月底的兔子的对数等于月初兔子对数加上刚成熟的新生兔子的对数。是一个结构体,其字段名称是对象属性名称,字段值是对应的属性值。
2023-07-03 16:58:34
515
原创 MATLAB Onramp
目录任务音频频率以表的形式导入数据循环恒星运动 知识点摘要 向量 现在包含 个点。要将这些点转换为频率,您可以将整个向量乘以采样频率 (),然后除以点数 ()。 将包含从 到 的频率。主频位于 的开头位置。您可以使用 函数放大所关注的区域。以表的形式导入数据 恒星运动 将 x = 、y = 处的点绘制成一个标记大小 () 为 的红色方框 (),添加到现有图中。知识点摘要MATLAB Course
2023-07-01 18:34:41
590
原创 MATLAB文化算法
此 MATLAB 函数 返回 A 的数组元素的乘积。如果 A 是向量,则 prod(A) 返回元素的乘积。如果 A 为非空矩阵,则 prod(A) 将 A 的各列视为向量,并返回一个包含每列乘积的行向量。任何满足文化算法结构的算法都能融入种群空间 中,通过相互融合,可以弥补传统算法的不足。此 MATLAB 函数 返回数组 A 中的元素数目 n 等同于 prod(size(A))。生成一个均匀随机数的数组,其中sz1,...,szN表示每个维度的大小。的每个元素四舍五入到指定单位时间的最接近的数。
2023-06-14 11:36:00
1972
1
原创 蚁群优化算法
蚁群优化(ACO)是一套概率元启发法和智能优化算法,其灵感来源于蚂蚁的社会行为。ACO算法也被归类为群集智能方法,因为在这些算法的结构中通过模拟蚂蚁的行为来实现这一范式。中的第一个其大小不等于 1 的数组维度开始返回。可选择性地使用上述任何语法指定方向。指定在上述任意语法的计算中包括还是忽略。的第二个维度返回其中各行的累积和。返回由结果的线性索引组成的列向量。,即查找与非零元素对应的前。中的非零元素对应的最后。中的非零元素对应的前。每列的累积和的矩阵。
2023-06-12 10:28:32
838
原创 模拟退火(SA)算法
在 x 轴上使用线性刻度、在 y 轴上使用以 10 为底的对数刻度来绘制 x 和 y 坐标。在同一组坐标轴上绘制多对 x 和 y 坐标。此语法可替代将坐标指定为矩阵的形式。返回一个数组,该数组在其行维度和列维度包含。要在同一组坐标轴上绘制多组坐标,请将。使用指定的线型、标记和颜色创建绘图。要绘制由线段连接的一组坐标,请将。中的至少一个指定为矩阵。指定为相同长度的向量。生成均匀分布的随机数。生成正态分布的随机数。生成连续统一的随机数。
2023-06-11 21:47:31
752
原创 DBSCAN 集群
对于任何一个相邻点,如果它的ε-相邻点包含预定数量的点,聚类就会被扩展到包含其相邻点。然而,对于未分配的点,如果邻域中的点的数量少于预定的阈值,则该点被认为是一个噪音。以三列数组形式返回 hsv 颜色图,其中包含的行数与当前图窗的颜色图相同。D = pdist2(___,Name,Value)指定了一个额外的选项,除了前面语法中的任何参数外,还使用了一个名-值对参数'Smallest'或'Largest'。大括号,用于cell型的数组(就是前面讲的单元数组)的分配或引用。相同的、由逻辑值“0”构成的数组。
2023-06-10 16:51:02
889
原创 使用 GMDH 进行时间序列预测
sz 必须至少包含 2 个元素,prod(sz) 必须与 numel(A) 相同。此 MATLAB 函数 返回逻辑值 1 (true),否则返回逻辑值 0 (false)。空数组、表或时间表有。的每个元素四舍五入为最近的整数。至少一个长度为 0 的维度,如 0×0 或 0×5。指定为正整数向量,以一次查询多个维度长度。返回由表中的行数和变量数组成的二元素行向量。返回一个数组,该数组在其行维度和列维度包含。还会为上述任意语法返回一个索引向量的集合。的第二个维度和第三个维度的长度。的副本在每个维度中如何排列。
2023-06-08 19:10:58
1043
原创 MATLAB应用
如果二者相同,函数将返回 1(true),否则返回 0 (false)。例如,如果已启用暂停功能,oldState = pause('off') 会在 oldState 中返回 'on' 并禁用暂停。使用智能聚类方法:(a)k均值算法,(b)模糊c均值聚类(FCM)和(c)自组织神经网络,使用 RGB 和 HSV 颜色编码来执行聚类任务。例如,指定余弦距离、使用新初始值重复聚类的次数或使用并行计算的次数。的各列视为向量,并返回一个包含每列乘积的行向量。将返回在对话框中选择的筛选器的索引。
2023-06-07 20:45:42
2047
1
原创 中国人口增长预测
您可以使用上述语法中的任何输入参数。除外)自动确定的,该小写字母字符不属于通过多个字母字符构成的字词的一部分。2005年为t=51,作为终止时刻。中的预测变量的非线性回归的估计系数。) 传递到函数句柄来计算该值的平方,与您将输入参数传递到标准函数一样。2005年为终止年,t=42。如果找到两个字符,则选择在字母表中靠后的字符。指定输入参数的内联函数。③将1980年看为起始年,2005年视为终止年。包括函数的输入参数。,并显式返回单个输出,即大小与包含平方值的。指定的模型,返回一个向量,其中包含。
2023-06-05 20:58:39
3654
1
原创 最优捕鱼策略
随着k值的增加,n值越来越小,即捕捞强度系数越大,鱼的年产卵量越小。捕捞强度系数k采用17.36,此时最大捕获量为3.887e11。当k为17.58时,f的最大值为1.6056e12。在实现可持续收获的前提下,追求最大产量或最佳效益。
2023-06-01 20:33:59
2617
9
原创 旅行商问题TSP
此 MATLAB 函数 在向量 x 和 y 指定的位置创建一个包含圆形标记的散点图。要在同一组坐标区上绘制多组坐标,请将 x 或 y 中的至少一个指定为矩阵。此 MATLAB 函数 按列顺序将 formatSpec 应用于数组 A1,...An 的所有元素,并将数据写入到一个文本文件。此 MATLAB 函数 返回一个数组,该数组在其行维度和列维度包含 A 的 n 个副本。此 MATLAB 函数 返回行向量,其中包含从 1 到 n 没有重复元素的整数随机排列。网格能量的最小值对应于最佳或者次最佳的路径距离。
2023-05-31 11:50:01
743
原创 投资收益与风险
②在短时期内所购买的各种资产(如股票,证券等)不进行买卖交易,即在买入后就不再卖出。其中r表示平均收益率,q表示风险损失率,p为交易费费率,x为投资金额。④在短时期内所给出的平均收益率、损失率和交易的费率不变。③在投资的过程中,无论盈利与否必须先付交易费。①每种投资是否收益是相互独立的。
2023-05-30 21:48:47
141
原创 水质评价与预测
接收一个N乘Q的矩阵,X和可选的Y的每一行的最小值和最大值,YMIN和YMAX,并返回一个N乘Q的矩阵,Y,和一个允许一致处理值的处理设置,PS。mapminmax通过将每行的最小值和最大值归一化为[YMIN, YMAX]来处理矩阵。将参数作为一个结构: FP.ymin, FP.ymax。是模糊神经网络的关键。
2023-05-29 23:31:48
113
原创 PFC落石模拟
根据Taghavi(2011)的定义,与光滑度的角度测量相对应的距离fdistance为0 < fdistance < 180。销毁,删除落在模型域之外的球、团块和壁面,这些球、团块和壁面的单个外延(即紧紧包括物体的轴对齐的边界框)。可以指定任何数量的分布的半径范围、体积分数和分布类型。当球或团块的中心点落在模型域外时,它们会被传回模型的另一侧。:块状模板中保留的最小和最大卵石的比率,0 < fratio < 1。:这个分布中的团块的体积分数。表示如果存在指定的名为s的前一个文件,那么它将被覆盖。
2023-05-26 18:28:04
2673
2
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人