- 博客(59)
- 资源 (2)
- 收藏
- 关注
RSS订阅原创 牛顿迭代法和牛顿-雅可比迭代法求解非线性方程或方程组(各一个简单例子)
【代码】牛顿迭代法和牛顿-雅可比迭代法求解非线性方程或方程组(各一个简单例子)
2025-05-24 17:35:11
206
原创 对colormap的那个N行三列矩阵进行线性插值的代码(三维scatter3中的确可以看到是正确插好值了的)
【代码】对colormap的那个N行三列矩阵进行线性插值的代码(三维scatter3中的确可以看到是正确插好值了的)
2024-12-07 16:08:15
133
原创 给定正多边形外切圆半径和边数，给出正多边形各顶点的xy坐标的matlab函数
给定正多边形外切圆半径和边数，给出正多边形各顶点的xy坐标的matlab函数
2024-06-15 20:45:31
210
原创 群折射率$ n_g$与有效折射率$ n_{eff}$的关系详细推导(假设$ \lambda_0$为真空入射光波长)
群折射率和有效折射率关系公式的详细推导
2023-08-15 15:01:50
4036
原创 二维体光子晶体的平面波展开法代码2
这篇文章接前面《二维体光子晶体的平面波展开法代码》这个是实际网格划分，各标记的含义上篇文章有说明下面是傅里叶系数恢复的相对介电常数分布(虚部)，看起来非常大，那是因为imag(EpsMat)趋于0的原因。下面这个图是实空间相对介电常数分布的伪彩图：下面这个图是TE能带图，可以看到和FEM的结果还是有点偏差，但实际上只要加密网格同时遵守前一篇文章所述的nG参数和precisStruct参数的关系就可以获得误差更小的能带图。
2023-05-10 15:35:15
449
原创 matlab实现DFT(取了一些数据计算结果和fft相差很小，但未详细验证)
DFT的matlab实现，和fft进行了少量对比误差非常小，但没详细验证，慎重。
2023-05-06 21:56:02
197
原创 (PWEM前置知识)解析公式求相对介电常数分布函数的 Fourier 展开系数以及基于此的绘制实空间原胞的相对介电常数分布
(PWEM前置知识)解析公式求相对介电常数分布函数的 Fourier 展开系数以及基于此的绘制实空间原胞的相对介电常数分布
2023-05-03 15:21:32
577
原创 Matlab绘图为什么设置了相同的字体大小结果出的图字体大小还是不一样
matlab相同字号设置但最终图片字号大小不一致(matlab2018b中已验证过)
2023-02-26 10:06:50
1268
1
原创 检查一维数组中的重复元素以及每个重复元素的重复次数(结果以矩阵形式给出)
matlab函数,检测一列数组中的重复元素和重复次数，结果以矩阵形式给出，第一行是重复元素，第二行是对应元素的重复次数
2022-07-27 21:38:42
620
原创 matlab复误差函数的高效计算(有官方背书)
matlab复误差函数的高效计算，可以和mathematica的官方函数对比验证下，有权威官方背书，是参考网站的代码改写而来。
2022-06-27 17:40:46
677
原创 赫尔维兹_勒奇超越函数（matlab自编函数）
function [out] = Lerch_fun(z,s,v)%Lerch_fun(exp(1j*2),3,0.1)% Lerch transcendent%表达式：\Phi(z, s, v)=\sum_{n=0}^{\infty} \frac{z^{n}}{(n+v)^{s}}%syms tf=@(t)((t.^(s-1)).*exp(-v*t))./(1-z.*exp(-t));out=(1/gamma(s))*(integral(f,0,Inf));end和Mathematic
2022-03-01 21:29:14
1006
原创 带选择程序（未测试）
function [band_sto]=choice_band(vec,val,band_index)%根据band_index选择能带%vec,val是最原始的vec,valval=real(sqrt(diag(val)));%w/c,，real可改，决定画能带的实部还是虚部[~,order]=sort(val);%w/cband_sto=zeros(size(vec,1),length(band_index));for count=1:length(band_index)tp=band.
2022-01-05 17:14:49
234
原创 三角格子BZ离散代码(T-K-M-T)（好像高对称点重复了）
function [k0]=BZ_dispersed()global a0delta=0.01;%可调参数T=[0,0];M=[0,2*pi/(sqrt(3)*a0)];K=[2*pi/(3*a0),2*pi/(sqrt(3)*a0)];TK=K-T;tk=TK(1)+TK(2)*1i;list_tk=(T(1)+T(2)*1i)+[0:delta:1]*tk;KM=M-K;km=KM(1)+KM(2)*1i;list_km=(K(1)+K(2)*1i)+[0:delta:
2022-01-01 17:31:14
307
原创 matlab画圆函数
function [] = circle(x,y,r)%画圆函数% circle(0,0,4) ;grid on% x y是中心，r是半径rectangle('Position',[x-r,y-r,2*r,2*r],'Curvature',[1,1])axis equal%为了修饰曲线的颜色，宽度，圈盘填充颜色等，可以设置其他参数等，例如%'edgecolor'，'b'，其中edgecolor表示边框颜色，后面的b是颜色参数值；%facecolor'，'r'，其中facecolor表示内部填
2021-12-31 11:24:04
6390
智能优化算法—免疫算法笔记.md
2020-09-14
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人