- 博客(71)
- 收藏
- 关注
RSS订阅原创 python读写改保存txt文件
我在网上能扒到的代码都是给基本的读入和保存，而我的需求是读入，删掉某一行的内容，在原位置添加相应的内容，找了一段时间没有找到对应的代码，就自己写了一个。代码的注释已经解释得比较清楚了，这个复制后，把对应的文件路径改一下是可以直接用的。
2023-05-13 10:18:18
652
1
原创 Ansys热仿真知识记录
这里举一个简单的例子，由于定子和转子之间是有气隙的，而在进行电磁仿真的时候，空气我们是不会单独进行建模的，所以这也导致了从maxwell里面导出来的模型是没有“空气”这部分3D模块的，我们在spaceclaim里面进行补充，把这部分建出来，如下图所示，黄色高亮部分为新建立的空气的3D模型，包含气隙和槽内气体，这样在进行下一步赋值的时候，热仿真部分才不会报错。是带有绕组端部的3D模型，图片看上去不太清晰，没有关系，我们之后的几何模型是要用spaceclaim去处理的，这里只是为了将模型导出。
2022-09-23 16:55:30
5329
原创 Themal稳态热仿真介绍
简单记录一下吧，themal这个软件虽然只能对稳态热进行仿真，但优势在于它不需要复杂而繁琐的操作，且仿真时长较传统的瞬态有限元仿真软件有大幅度的缩小。首先，在spaceclaim里面建模，我们先建一个2D模型把流程说清楚。模型如下：中间两个圆是铜线，其他几层是绝缘材料（或者换成其他材料），这个很方便设置的。几何模型绘制的时候先画草图，然后电机“PULL”按钮，它就生成面了，接下来把铜和绝缘层的几何构图分开，左侧的菜单树是这个样子的,需要把各组分构成一个component，这样后续仿热的时候才不抱错，comp
2022-07-01 14:32:17
1132
原创 800V电机绝缘与局部放电学习笔记4
目录前言：电机绝缘分类：过电压：上一篇我们讨论了局部放电发生的判准及它的危害，本篇接着讨论电机的绝缘系统与局部放电以及电机端部过压之间的关系。电机绝缘可分为I型散绕结构和II型成型结构。工作电压以700V作为划分。700V以下的，绝缘采用有机材料组成，当发生局部放电现象时，会破坏有机物的分子结构，导致电机绝缘失效。所以在I型结构里面，我们需要控制电机的PDIV大于电机额定运行的电压，设计时通常是留有一定余量的，至于电机PDIV怎么设计和测试，将在之后的笔记中进行介绍。对于额定电压大于700V的，绝缘结构一般
2022-06-24 09:18:41
884
1
原创 800V电机绝缘与局部放电学习笔记3
这部分内容今天下午请教了一下师兄，感觉还有很多东西需要学习，确实是一个很大的方向。前面两篇专栏分别对问题的基本概念和激励电压源进行了讨论，其实我们的目的还是要始终盯着如何把匝间的电应力给预估准了，然后再去判断它是否会发生局部放电现象。本篇内容，我们讨论局部放电的评判标准。帕邢定律，表征均匀电场气体间隙击穿电压、间隙距离和气压间关系的定律。通常写成V=f(pd),p是气压，d是两板间的距离。1889年由F.帕邢根据平行平板电极的间隙击穿试验结果得出。绝缘体的击穿电压（英语：breakdown voltage）
2022-06-21 09:55:28
512
原创 800V电机绝缘与局部放电学习笔记2
第一部分我们介绍了一些基本概念，这篇对于波反射的机理做一个记录，再顺带着说一下高开关频率为什么会让电机的绝缘承受更大的电压应力。也就是对于过压问题中的“源”进行简单的说明和整理。传输线上任意点的电压和电流都是入射波与反射波的叠加，因此通常用反射系数描述反射波与入射波之间的幅度与相位关系。反射系数分为电压反射系数与电流反射系数。我们这里讨论的是过压问题，就关心电压反射系数就行了。简而言之，如果没有波的反射，那么我们输入端给的是多少电压，电机绕组端部电压就是多少，但是这个波到了端部，还要往回走，他们两一叠加，这
2022-06-20 16:33:42
504
原创 800V电机绝缘与局部放电学习笔记1
这部分内容我已经断断续续地学习了一段时间了，也尝试搭建自己的知识架构，但是终究觉得很不扎实，虽然也跟自己才接触这方面没多久有关系，思来想去还是打算写个专题吧，如果记录的笔记存在错误，欢迎大家跟我讨论，如果能帮助到更多的人，那是最好不过了。PWM控制电驱系统中，过冲电压Vos（俗称“电压振荡”）的形成很大一部分来自电回路中无法消除的杂散电感的存在。二极管在关闭的过程，杂散电感因变化的电流而产生感应电压。与传统的Si-IGBT相比，使用SiC-MOSFET会导致更高的dv/dt，从而导致更高的过冲电压。SiC器
2022-06-18 15:23:48
957
原创 matlab fourier变换反变换
求解析表达式的时候，用到的比较多。其实最近在解泊松方程嘛，确实比较头疼这个问题，还好有伙伴一起搞，先简单记录一下这部分代码吧，给的示例，很简单的。%%%test for fouriersyms x t;A=exp(-x^2);test=fourier(A,x,t);b=ifourier(test,t,x);分别做了傅里叶变换和反变换。但是，傅里叶变换需要函数满足狄利克雷条件，并且在整个积分区间满足绝对可积，条件过于苛刻，很多时候就会出现这样的结果：fourier(exp(- k*
2022-04-29 21:22:00
996
1
原创 matlab-heatmap图代码
这个吧，嗯，简单记录一下吧。a = rand(5,5);xname={'x1','x2','x3','x4','x5'};yname={'y1','y2','y3','y4','y5'};h = heatmap(xname,yname,a);h.CellLabelFormat = '%0.3f'; %展示精度% colormap(gca, 'parula') %改颜色colormap(summer);% sortx(h,'y2'); %默认升序排序
2022-04-17 21:41:06
4374
1
原创 matlab矩阵维度转换
其实是基本操作。。。。但我懒得每次都去想一遍，就记录一下。问题是这样的比如我现在有一个矩阵是46*22的，但是我想要把它变成一维的1012*1。就很简单.[x,y]=size(Bw); %确定矩阵维度Bre=zeros(x*y,1);for i=1:y for j=1:x Bre((i-1)*x+j)=Bw(j,i); endend额，代码不解释了。...
2022-04-15 08:30:16
1690
原创 matlab配色
我来当一次搬运工，人家写得太好了（知乎上的）。。。。。我记录一下Bw是我的数据哈，这个大家不要深究了，主要看绘图代码：plot(Bw,'color',[252/255 41/255 30/255],'linewidth',2);代码本身很简单的，后面的矩阵就是三原色比例调控嘛，关键是这个三原色数据咋来的，看这里：Matlab绘图颜色搭配 - 知乎就很香了。然后把图片给大家看看，这个红色与默认的红色有差别的。...
2022-04-11 10:33:00
3108
原创 Matlab基本作图汇总
简单整理一下吧，说明在代码注释里面都有。subplot(221)x=0:0.002:2*pi;y=sin(x);plot(x,y,'linewidth',2) %线xlabel('x');ylabel('y');set(gca,'fontsize',15,'fontname','Times'); %设置坐标轴字体subplot(222)area(x,y) %面xlabel('x');ylabel('y');set(gca,'fontsize',15
2022-04-07 22:10:00
1224
原创 matlab定义变量，对函数积分
很实用的小技巧，每天学习一点点。syms x y; %定义变量x yy=sin(x); %定义函数int(y,x,0,pi) %对函数从0到pi进行积分输出结果：ans=2
2022-04-07 20:08:49
2130
原创 matlab-map图
一直觉得这种图很好看。。。。可能会用到吧[xx,yy]=meshgrid(1:5000,1:100); %x方向为1,2,3，…50，y方向为1,2,3，…50c=xx+yy-0.5*xx;pcolor(xx,yy,c)colorbarshading interp;saveas(gcf,'picex.png')稍微说一下meshgrid的吧，就是咱这种图它除了xy之间有对应关系之外，还有第三个变量c也是有联系的，所以我们这个图，其实是表征的xyc之间的关系。之前写过的用matlab画
2022-04-01 12:02:32
2800
原创 Matlab点画线
这个作图和python还是有点区别的，似乎对命令输入的顺序还有要求。t=[1190.21153.141071.561069.221063.18];w=10:10:50;% scatter(w,t,sz,'r','filled');plot(w,t,'o--','linewidth',2)hold ont1=[1073.021057.811129.71028.181015.6];plot(w,t1,'*--','linewidth',2) %调整线宽，先给线的类型，再
2022-03-14 11:08:20
6125
原创 Matlab散点图绘制及颜色渐变补充
data=readmatrix('72验证点.xlsx');x=data(:,12);y=data(:,11);scatter(x,y,'r','filled');legend('72slots')xlabel('torque per volumn')ylabel('efficiency')saveas(gcf,'opt72.png')
2022-03-09 11:05:50
2465
1
原创 matlab散点图渐变颜色绘制
贴个官方说明：散点图 - MATLAB scatter- MathWorks 中国然后是我的代码：x=data(:,10);y=data(:,9);c=data(:,6);sz = 7; %点大小scatter(x,y,sz,c,'filled')h =colorbar; %渐变色注释栏xlabel('torque per volumn');ylabel('efficiency')set(get(h,'label'),'string','statorDia(m
2022-03-06 19:57:31
6440
原创 Python-另一作图代码风格
这个代码风格跟我之前的不太一样哈，是我一个学弟写的，写得很好，我转载出来分享记录一下，作复杂图像更适合用这种代码风格。import numpy as npimport matplotlib.pyplot as pltimport pandas as pdfrom matplotlib import font_manager#读取实验数据data = pd.ExcelFile(r'488result.xlsx')#data = pd.read_execl(r'C:\Users\Admini
2022-02-24 22:25:00
508
原创 Python-柱状图-横坐标为字符
import numpy as npfrom scipy.fftpack import fftimport pandas as pdimport matplotlibmatplotlib.rcParams['backend'] = 'SVG'# import seaborn# import math as mtimport matplotlib.pyplot as plt# import seaborn as sns# sns.set(style='darkgrid')y=[150,.
2022-01-06 16:21:03
6350
原创 matlab-柱状图汇总操作
fs=500:15:590;yall=[aa;dd];bar(fs,yall);xlabel('Fre(Hz)');ylabel('loss(W)');legend('FEM','Ana','Location','northwest');saveas(gcf,'procompare.png')aa和dd分别是两组柱状图的纵坐标数据。结果如下：
2021-12-23 19:33:43
4500
2
原创 Python-坐标轴、图例字体大小
图片坐标大小很重要。就是那句注释的代码就okplt.plot(time,B,'--',label='Conductor 1',linewidth=2)plt.plot(time,B1,'-',label='Conductor 2',linewidth=2)plt.xlabel('Time(ms)', fontsize = 14)plt.ylabel('B(T)', fontsize = 14)plt.tick_params(labelsize=13) #调整坐标轴数字大小plt.le
2021-12-21 16:44:05
9831
原创 Python饼图、调整字体大小
推荐两个非常好的博客内容：Python实现可视化（二）--label字体大小调整、图例设置https://blog.csdn.net/weixin_39602776/article/details/111314454Python绘制饼图调节字体大小、防止标签重叠解决方法 https://www.cnblogs.com/mengxiaoleng/p/12804117.html然后是我自己的饼图代码，里面用到了字体大小调整这部分的功能：import matplotlibimport matplotl
2021-12-12 12:57:40
9106
转载 latex-特殊符号
转载于我的小伙伴，记录的很全，大家可以去康康Latex数学符号对应表https://www.cnblogs.com/dingdangsunny/p/12312966.html
2021-12-01 22:28:16
299
原创 latex学习1
最近的课程老师建议用latex写课程报告，想一想一直用word排版的确也麻烦，还是学一下latex吧。简单做一个记录。\documentclass{article}%\documentclass{ctexart}\begin{document}tengfeiHU was walking on the street,he saw a beautiful girl.number $n=6$hu's money $m=sum/n$\begin{equation}|E|=2 \pi f
2021-11-28 22:07:16
153
原创 Matlab-矩阵存入txt
数据导出的时候会有这样的需求，简单记录一下：存储的矩阵I长这样是一个2501*4的矩阵。我想原封不动的存入txt里，代码如下：fid=fopen('Iva.txt','wt');[m,n]=size(I);for i = 1:1:m for j=1:1:n if j == n fprintf(fid,'%g\n',I(i,j)); else fprintf(fid,'%g\t',I(i,j)).
2021-11-25 22:35:45
3843
5
原创 热网络法-初步搭建
捣鼓一天多了，还是记录一下吧。前言热学相关的理论我前面也记录过两篇，之后因为有其他事，这方面学习就搁置了，现在重新捡起来，做一些初步的仿真学习，很感谢我的伙伴叮叮当当跟我一起讨论学习这方面的内容。之前其实摘抄过一些关于电机热网络的搭建过程，但那个问题比较复杂，不适合新手入门学习。今天分享一个简单的例子。正文问题是这样的，四个导体，材料是铜。导体通电流会发热，问稳态温度是多少？几何模型每根导体选取1m*1m*1m的尺寸，便于我们后面对导热系数和对流系数做换算。几何模型如下：其实就是
2021-11-23 20:15:47
3878
3
原创 Matlab-多项式拟合
推荐一个拟合及其方便的工具箱cftool，傻瓜式操作。首先，在matlab的工作区里有待拟合的数据。然后命令行输入cftool，就会出现如下界面：在x data和y data区选择拟合数据名称，这里和工作区是对应的。拟合的图像和结果就在下方展示出来了，十分方便。...
2021-11-11 20:45:29
1126
原创 Matlab-FFT
搬运工，完整的多版本FFT代码在这里，分享给大家。叮叮当当sunny-FFThttps://www.cnblogs.com/dingdangsunny/p/12573744.html#_label3function [Fre,Amp,Ph] = FFT(data,Fs,ampDB,isDetrend) % 快速傅里叶变换 % data:波形数据 % Fs:采样率 % ampDB:逻辑值，是否进行对数变换，默认为false % isDetrend:逻辑值，是否
2021-11-11 14:25:22
779
原创 Python-双坐标作图（加图例）
其实就是很简单的一个双坐标绘制，由于画了两条曲线，对应图列还是需要加的。data1=pd.read_csv('br.csv')data2=pd.read_csv('bt.csv')By=list(data1['Br []'])Bx=list(data2['Bt []'])fig = plt.figure(figsize=(10,6))ax1=fig.add_subplot(111)ax1.plot(By,label='Br')ax1.set_ylabel('Br(T)')ax2 = a
2021-11-07 21:44:55
4438
2
原创 Python-点加标注
现学现卖，当一次搬运工，我觉得写得挺好的。plt.annotate()函数解析（最清晰的解释）https://blog.csdn.net/qq_36387683/article/details/101377416下面是我的代码：drawloss=[11.432,0.495,0.0260,0.0020,0.0004,0.0003,]names = ['conductor1', 'conductor2', 'conductor3', 'conductor4','conductor5','conduc
2021-11-07 19:48:08
1703
原创 Matlab-标准遗传算法
写了两个月多目标了，突然上课课程作业叫写单目标。之前用Python写过一次，这次用matlab写吧，也暴露了遗传算法的一些固有问题，可以记录一下：先上代码，讨论放到最后：%%get initial parameterNP=40; L=10; Pc=0.7; Pm=0.1; G=50; %迭代次数 O2=[0 3.14];wid=[2 10]; %rotor maglen=[27 35];%% generate Primary populationfo
2021-11-06 19:58:21
965
原创 Matlab-作图垂线、文字标注
figure(2)plot(h2,kr,'LineWidth',2)hold onplot(thickness,kr1,'r*')line([1.65,1.65],[1,1.45],'Color','g','LineWidth',2) %垂直线text(1.65,1.1,' \leftarrow skin depth') %文字标注xlabel('thickness(mm)')ylabel('ratio AC to DC ')legend('Analyze','FEM')% titl.
2021-11-05 19:45:00
2170
原创 Matlab-读取CSV数据、作图
一些基本操作，记录一下：matlab比较依赖它的矩阵机制（这一点和Python很不同），特别是在处理数据的时候。比如在读取csv文件的时候，matlab是有csvread语句的，但是官方不推荐，我们就不使用了。上代码了，中间处理数据的代码我给省略掉了哈。data=readmatrix('skin.csv'); %读取数据%作图代码f1=100:100:1000;plot(f,kr,'-.','LineWidth',2);hold onplot(f1,ues,'r*')legend
2021-11-04 20:12:22
4750
原创 Python-代码分模块
这段时间matlab用多了，居然把python的分模块操作给忘了。。。卑微Python:#%%Matlab：%%绿色的这部分就是我们划分的模块。
2021-11-03 19:13:03
984
原创 Matlab-输出函数表达
利用matlab处理完数据之后，得到了对应的表达式，现在想把表达式输出出来怎么办呢，一条简单的语句：sprintf('y1=%.3f*sin(2*pi*%d*time)',is,f1)跟C很像哈。
2021-10-29 16:34:57
5272
原创 Matlab求定积分
操作很简单，Matlab yyds>> syms x %申明变量>> f=sin(x)^3*cos(x)^2 f = cos(x)^2*sin(x)^3 %定义积分函数 >> int(f,x,0,pi/2) %求解 ans = 2/15
2021-10-26 16:29:38
2670
原创 Python—图片大小调节及横坐标为文字的折线图
很多时候图片为了展示，是对作图的大小和坐标有要求的。这里咱就改一下。说明都在注释里面了。才发现今天是1024，跟着过个节，凑凑热闹。import matplotlib.pyplot as pltimport seaborn as snssns.set(style='darkgrid')y=[1.14,1.28,1.5,2]names = ['design1', 'design2', 'design3', 'design4'] #这是横坐标展示的文字x = range(len(names)
2021-10-24 10:24:27
4024
1
原创 电机振动噪声（NVH）气隙磁场推导
其实我之前振动噪声系列中，涉及气隙磁密的分析的，都是给一个推导的逻辑链，说什么知道磁动势了，知道气隙磁导了，就可以推出气隙磁密了。这话吧，不能说是错的，但是就跟大海一样，全是水。言归正传，之前的笔记旨在给大家搭建一个整体的理论框架，建立起这么一个印象就足够了。因为你拿着之前的笔记去细究，仔细一看其实哪都没说明白，但是振动噪声关心的问题好像又说了很多了。接下来，我想从源头上，给大家梳理梳理，电机振动噪声那么关心的气隙磁场究竟是怎么推出来的
2021-10-23 21:18:46
2820
4
原创 电机振动噪声（NVH）——整数槽分数槽谐波分析补充
在之前的博客里，对径向力做了分类。但是没有定量地对力的阶数做讨论。这篇博客就把这个坑给填上。其实数学的推导算不上困难，阶数推导的难点在于整数槽的电机和分数槽的电机它们的谐波次数不一样的，需要分开讨论。这一点也是困扰了我很长时间。
2021-10-20 16:16:30
5227
2
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人