摘要：
import pandas as pd def make_df(cols, ind): """一个简单的DataFrame""" # 字典推导式，为每列生成数据 data = {c: [str(c) + str(i) for i in ind] for c in cols} # 创建并返回 Data
阅读全文
posted @ 2024-10-08 16:58
redufa
阅读(48)
推荐(0)
摘要：
目录1. 从列表构造2. 从字典构造3. 从NumPy数组构造4. 从Series对象构造5. 从文件加载6.字典+数组 在Pandas中，可以通过多种方式创建DataFrame对象。以下是五种不同的创建DataFrame的方法： 1. 从列表构造 使用列表的列表（即嵌套列表）来创建DataFram
阅读全文
posted @ 2024-10-08 16:15
redufa
阅读(1088)
推荐(0)
摘要：
目录Pandas 方法总结1. 读取写入2. 数据查看3. 选择过滤4. 数据清洗5. 数据转换6. 数据聚合7. 数据合并8. 时间序列9. 数据处理10 数据重塑11. 数据导出12. 绘图13.其他方法 Pandas 方法总结 Pandas 是一个强大的 Python 数据分析库，它提供了快速
阅读全文
posted @ 2024-10-08 11:33
redufa
阅读(158)
推荐(0)
摘要：
1、伽辽金法 没有分布载荷 \[\frac{d}{d\hat{x}}(AE\frac{d\hat{u}}{d\hat{x}}) \]伽辽金法： \[\int_{0}^{L}\frac{d}{d\hat{x}}(AE\frac{d\hat{u}}{d\hat{x}}){N}_{i}d\hat{x} \
阅读全文
posted @ 2024-10-07 02:00
redufa
阅读(44)
推荐(0)
摘要：
形函数构造 构造单元1的一般近似函数 \(\overline{V(x)}^{(1)}\),由于该单元只有两个节点\(x_1\)和\(x_2\)，我们选择包含两个参数\(\alpha_1\)和\(\alpha_2\)的近似方程 \[\overline{V(x)}^{(1)}=\alpha_1+\alp
阅读全文
posted @ 2024-10-07 01:54
redufa
阅读(29)
推荐(0)
摘要：
[Q] 是两个坐标系的旋转矩阵，为正交矩阵。 \[\{e^{\prime}\}=[Q]\{e\} \\ \]\[[Q]^T=[Q]^{-1} \]对于在坐标系{e}中存在的向量 u 和v 存在关系 \[\{u\}=[a]\{v\} \]其在坐标系\(\{e^{\prime}\}\)存在关系 \[\{
阅读全文
posted @ 2024-10-07 01:54
redufa
阅读(28)
推荐(0)
摘要：
\(l_{ij}\) 是两个坐标系的旋转矩阵，为正交矩阵。 \[e_i^{\prime}=l_{ij}e_j \]对于在坐标系{e}中存在的向量 u 和v 存在关系 \[u_i=a_{ij}\nu_j \]其在坐标系\(\{e^{\prime}\}\)存在关系 \[u'_i=a'_{ij}\nu'_
阅读全文
posted @ 2024-10-07 01:53
redufa
阅读(31)
推荐(0)
摘要：
1、伽辽金法 没有分布载荷 \[\frac{d}{d\hat{x}}(AE\frac{d\hat{u}}{d\hat{x}}) \]伽辽金法： \[\int_{0}^{L}\frac{d}{d\hat{x}}(AE\frac{d\hat{u}}{d\hat{x}}){N}_{i}d\hat{x} \
阅读全文
posted @ 2024-10-07 01:51
redufa
阅读(133)
推荐(0)
摘要：
直接平衡法或刚度法\(^{2}\) 直接平衡法通过基本单元的力平衡条件以及力/变形关系，得出关联节点力和节点位移的刚度矩阵以及单元方程。由于直接平衡法最易应用于线单元或一维单元，可以分别将其应用于弹簧单元、杆单元和梁单元。 直接刚度法\(^{2}\) 节点力与节点位移的关系，此关系是刚度矩阵 \[\
阅读全文
posted @ 2024-10-07 01:51
redufa
阅读(162)
推荐(0)
摘要：
平面应力问题 平面应力问题的平面应力应变关系 \[\begin{Bmatrix}\sigma_{xx}\\\\\sigma_{yy}\\\\\tau_{xy}\end{Bmatrix}=\frac{E}{1-\gamma^2}\begin{bmatrix}1&\gamma&0\\\\\gamma&1
阅读全文
posted @ 2024-10-07 01:50
redufa
阅读(152)
推荐(0)
摘要：
二维或三维的分布积分（格林公式） 分布积分对下式积分 \[\int\int_{\Omega}\Phi\frac{\partial\Psi}{\partial x}\mathrm{d}x\mathrm{d}y \] 首先对变量\(x\)分布积分 \[\int\limits_{X_L}^{X_R}U\m
阅读全文
posted @ 2024-10-07 01:49
redufa
阅读(1216)
推荐(0)
摘要：
在自然坐标系中 , \(\xi_2=1\)和 \(\xi_2=1\),在物理坐标系中为 \(x_1\) 和\(x_2\)，相应的节点位移为\(u_1\) 和\(u_2\) 。 在自然坐标系 下，单元形函数为 \[N_{1}(\xi)=\frac{1}{2}(1-\xi)\\N_{2}(\xi)=\f
阅读全文
posted @ 2024-10-07 01:48
redufa
阅读(197)
推荐(0)
摘要：
变分方法 对连续介质问题，位置函数的\(u\)的泛函为 \[\Pi=\int_{\Omega}\boldsymbol{F}\Bigg(u,\frac{\partial u}{\partial x},\cdots\Bigg)\mathrm{d}\boldsymbol{\Omega}+\int_{\Ga
阅读全文
posted @ 2024-10-07 01:47
redufa
阅读(127)
推荐(0)
摘要：
第一章 1.1 流体的概念 任何固体材料都有一个强度极限，即使合外力和力矩都为零，它的内部也可能会存在着拉力、压力或者剪切力。当这些内应力超过了材料的强度极限时，固体就会被破坏，从而产生运动。微观上体现为断裂处的分子（或原子）之间的化学键被破坏，失去了相互的作用力，不再能保持原有结构形式了。 流体的
阅读全文
posted @ 2024-10-07 01:38
redufa
阅读(547)
推荐(0)
摘要：
目录3.1 几个基本概念3.3 任意斜截面上的应力3.4 主应力及应力（张量）不变量3.5 最大、最小正应力和最大剪应力 3.1 几个基本概念 • 外力 外力指的是我们熟知的机械力、电磁力等，物体因外力作用而变形。作用于物体的外力可分为体积力和表面 力，它们分别简称为体力和面力。 体力是分布在物体体
阅读全文
posted @ 2024-10-07 01:34
redufa
阅读(460)
推荐(0)
摘要：
目录一. 引言二. 塑性应变增量推导三. 弹塑性刚度矩阵推导四. 塑性模量理解五. 小结 一. 引言 弹塑性理论定义材料在荷载作用下的变形是弹性变形和塑性变形之和，其中研究塑性变形需要解决三个方面的问题： ①产生塑性变形的起点； ②产生塑性变形的方向； ③产生塑性变形的大小。 在塑性理论中，描述以上
阅读全文
posted @ 2024-10-07 01:25
redufa
阅读(567)
推荐(0)
摘要：
目录屈服条件2.1 应力偏张量及性质1.应力张量的分解及应力偏量2.应力偏张量的性质2.2 应力空间，π平面和Lode参数1.主应力空间和\(\pi\)平面2. 应力偏量的二维表示2.3 应力偏量和等效应变2.4 初始屈服条件和初始屈服曲面1.屈服条件的一般概念2.屈服条件的简化及屈服面的几何形状性
阅读全文
posted @ 2024-10-07 01:19
redufa
阅读(547)
推荐(0)
摘要：
例1 计算积分 \[I=\int_Cx^2ydx-xy^2dy, \]其中C是上半圆 \(\begin{aligned} & \text{ }x^2+y^2=a^2,y\geqslant0,\text{ }\\ & \end{aligned}\) 逆时间方向 \[\begin{aligned} &
阅读全文
posted @ 2024-10-07 01:02
redufa
阅读(32)
推荐(0)
摘要：
张量（tensor） 这一术语最初是用来描述弹性介质各点应力状态的，后来发展成为力学和物理学的一个有力数学工具，目前力学方面的理论性文献都不同程度地这用了这一工具 由坐标原点和三条不共面的标架直线构成的坐标系称为直线坐标系，如果三标架直线上的单位尺度相同，称为笛卡尔坐标系，否则称为仿射坐标系。 笛卡
阅读全文
posted @ 2024-10-07 01:01
redufa
阅读(269)
推荐(0)
摘要：
二阶非线性自治系统绘制相平面图。假设我们有一个简单的阻尼摆系统，其状态方程可以表示为： \[ dx1/dt = x2 \\ dx2/dt = -cx2 - gsin(x1) \] import numpy as np import matplotlib.pyplot as plt from scip
阅读全文
posted @ 2024-10-04 19:51
redufa
阅读(97)
推荐(0)
摘要：
start_points 和 end_points 数组分别表示向量的起点和终点。 使用 plt.quiver 函数绘制向量（箭头）。 plt.scatter 用于绘制起点和终点（可选）。 plt.legend 添加图例。 plt.title、plt.xlabel 和 plt.ylabel 添加标题
阅读全文
posted @ 2024-10-04 19:37
redufa
阅读(312)
推荐(0)
摘要：
\[\begin{aligned}\frac{\mathrm{d}z\left(t\right)}{\mathrm{d}t}&=Az\left(t\right)+Bu\left(t\right)\\y\left(t\right)&=Cz\left(t\right)+Du\left(t\right)\
阅读全文
posted @ 2024-09-30 17:37
redufa
阅读(237)
推荐(0)
摘要：
单输入单输出(SISO)系统 二阶振动方程 \[m \frac{\mathrm{d}^{2}x\left(t\right)}{\mathrm{d}t^{2}}+b \frac{\mathrm{d}x\left(t\right)}{\mathrm{d}t}+kx\left(t\right)=f(t)
阅读全文
posted @ 2024-09-30 17:37
redufa
阅读(510)
推荐(0)
摘要：
根据状态空间方程的一般表达式，求解矩阵形式的微分方程可以掌握系统状态变量随时间的变化的解为 \[z\left(t\right)=\mathrm{e}^{A\left(t-t_{0}\right)}z\left(t_{0}\right)+\int_{t_{0}}^{t}\mathrm{e}^{A\le
阅读全文
posted @ 2024-09-30 17:36
redufa
阅读(235)
推荐(0)
摘要：
相平面与相轨迹(Phase Portrait)将使用直观的图形来分析微分方程，特别是二阶微分方程。相轨迹描述了系统的状态变量随时间在相平面上的变化轨迹。它的理念也可以拓展到更高维度的系统中。而且不只是线性系统，在非线性系统中也可以利用这一数学工具分析系统的表现。 相平面与相轨迹(Phase Port
阅读全文
posted @ 2024-09-30 17:34
redufa
阅读(481)
推荐(0)
摘要：
matlab 画微分方程 相平面图 在MATLAB中，可以使用quiver函数来绘制微分方程的相平面图。相平面图是用于展示动态系统中状态变量变化的一种图形表示方法，特别适用于二阶微分方程。 以下是一个简单的例子，展示如何在MATLAB中绘制一个线性微分方程的相平面图： % 定义微分方程 dx/dt
阅读全文
posted @ 2024-09-30 14:43
redufa
阅读(870)
推荐(0)
摘要：
机电系统的相似性 物理量符号说明 \(\bullet\) 通过型变量：\(F=\)力，\(T=\)扭矩，\(i=\)电流，\(Q=\)流体体积流速，\(q=\)热流速 \(\bullet\) 跨越型变量：\(v=\)平动速度，\(\omega=\)角速度，\(v=\)电压，\(P=\)压强，\({T
阅读全文
posted @ 2024-09-27 17:18
redufa
阅读(247)
推荐(0)
摘要：
典型环节的幅相特性曲线 1. 比例环节 传递函数 \[G(s)=K \]其频率特性为 \[G(\mathrm{j}\omega)=K+\mathrm{j}0=K\mathrm{e}^{\mathrm{j}0} \]\[A(\omega)=\mid G(\mathrm{j}\omega)\mid=K\
阅读全文
posted @ 2024-09-24 18:02
redufa
阅读(470)
推荐(0)
摘要：
阅读全文
posted @ 2024-09-21 23:55
redufa
阅读(21)
推荐(0)
摘要：
目录 目录目录参考文献1 李荣华 微分方程的数值解法零、 偏微分方程偏微分方程形式抛物线方程：扩散方程双曲线方程：对流方程椭圆方程(1):拉普拉斯方程椭圆方程(2):泊松方程一、常微分方程初值问题常微分方程形式1 常微分方程求解(1)：欧拉法1.1 欧拉法与改进欧拉法matlab代码1）欧拉法2）欧
阅读全文
posted @ 2024-09-21 23:43
redufa
阅读(311)
推荐(0)
摘要：
例1 计算积分 \[I=\int_Cx^2ydx-xy^2dy, \]其中C是上半圆 \(\begin{aligned} & \text{ }x^2+y^2=a^2,y\geqslant0,\text{ }\\ & \end{aligned}\) 逆时间方向 \[\begin{aligned} &
阅读全文
posted @ 2024-09-21 22:38
redufa
阅读(41)
推荐(0)
摘要：
终于成功解决PDF文档批量添加书签的问题 需要的软件 一、常规的pdf阅读器； 二、PDFPatcher（一个好用的第三方软件） 三、EXCEL (用到的函数：RIGHT(), LEFT(), FIND(), 和连接符& 1、利用pdf阅读器，福昕或者Adobe，随便添加一两页书签。（作为后续的书签
阅读全文
posted @ 2024-09-21 16:00
redufa
阅读(659)
推荐(0)
摘要：
从博客园导出md格式，让typora去打开，遇到其中公式标识符的前后换行没有了，用vscode去渲染没事，为了能够正常渲染，用python给前后添加换行符。 import os import re def process_md_file(file_path): with open(file_path
阅读全文
posted @ 2024-09-21 00:02
redufa
阅读(24)
推荐(0)
摘要：
正太分布画图 close all t=linspace(-4,4,1000); pdf=normpdf(t) cdf=normcdf(t) plot(t,pdf,"b",'LineWidth', 2); hold on plot(t,cdf,'k','LineWidth', 3); 生成均匀分布随机
阅读全文
posted @ 2024-09-19 20:49
redufa
阅读(79)
推荐(0)
摘要：
一、获取当前路径 1、使用sys.argv[0] import sys print sys.argv[0] #输出 #本地路径 2、os模块 import os print os.getcwd() #获取当前工作目录路径 print os.path.abspath('.') #获取当前文件目录路径
阅读全文
posted @ 2024-09-19 18:42
redufa
阅读(3711)
推荐(0)
摘要：
' 七个使用公式的例子 ' '1、行标是变量的例子 '2、列标用一维数组输入字符串做列标，行标是变量 '3、用RC做相对引用 '4、RC和变量混合 '5、整列使用RC引用，避免使用循环 '6、整个区域RC引用——相对引用/循环 '7、整个区域RC应用——相对引用/不用循环 Sub 数例() [a1]
阅读全文
posted @ 2024-09-19 18:33
redufa
阅读(101)
推荐(0)
摘要：
频率响应 公式推导 正弦输入的一般形式 \[u(t)=Asin(\omega_i)+Bcos(\omega_i) \]整理 \[u(t)=M_i sin(\omega_i t + \phi_i) \\ 其中 \phi_i=arctan \frac{B}{A} \qquad ,M_i=\sqrt{A^
阅读全文
posted @ 2024-09-19 13:15
redufa
阅读(331)
推荐(0)
摘要：
#include <stdio.h> int main(viod) { int num1,denom1,num2,denom2,result_nun,result_denom; printf("Enter first fraction: "); scanf("%d/%d",&num1,&denom1
阅读全文
posted @ 2023-12-03 16:16
redufa
阅读(9)
推荐(0)
摘要：
【【公开课】Python股票自动交易从零开始~】https://www.bilibili.com/video/BV1SW411A7Ab?p=6&vd_source=056bd9dc74b57a861c5ac342ecab8bbc 1 import requests 2 import pandas 3
阅读全文
posted @ 2023-11-27 21:30
redufa
阅读(145)
推荐(0)
摘要：
重命名有限元文件，提取关键词所在行 import shutil import csv oldname = 'static_linear.bdf' keywords = ['GRID', 'CQUAD4','SPC,','FORCE,'] newname=oldname[0:len(oldname)-
阅读全文
posted @ 2022-03-26 22:31
redufa
阅读(72)
推荐(0)
浙公网安备 33010602011771号