微环谐振器 matlab 代码

合集下载

matlab命令大全 全面

matlab命令大全 全面

分类:VC MatlabA aabs 绝对值、模、字符的ASCII码值acos 反余弦acosh 反双曲余弦acot 反余切acoth 反双曲余切acsc 反余割acsch 反双曲余割align 启动图形对象几何位置排列工具all 所有元素非零为真angle 相角ans 表达式计算结果的缺省变量名any 所有元素非全零为真area 面域图argnames 函数M文件宗量名asec 反正割asech 反双曲正割asin 反正弦asinh 反双曲正弦assignin 向变量赋值atan 反正切atan2 四象限反正切atanh 反双曲正切autumn 红黄调秋色图阵axes 创建轴对象的低层指令axis 控制轴刻度和风格的高层指令B bbar 二维直方图bar3 三维直方图bar3h 三维水平直方图barh 二维水平直方图base2dec X进制转换为十进制bin2dec 二进制转换为十进制blanks 创建空格串bone 蓝色调黑白色图阵box 框状坐标轴break while 或for 环中断指令brighten 亮度控制C ccapture (3版以前)捕获当前图形cart2pol 直角坐标变为极或柱坐标cart2sph 直角坐标变为球坐标cat 串接成高维数组caxis 色标尺刻度cd 指定当前目录cdedit 启动用户菜单、控件回调函数设计工具cdf2rdf 复数特征值对角阵转为实数块对角阵ceil 向正无穷取整cell 创建元胞数组cell2struct 元胞数组转换为构架数组celldisp 显示元胞数组内容cellplot 元胞数组内部结构图示char 把数值、符号、内联类转换为字符对象chi2cdf 分布累计概率函数chi2inv 分布逆累计概率函数chi2pdf 分布概率密度函数chi2rnd 分布随机数发生器chol Cholesky分解clabel 等位线标识cla 清除当前轴class 获知对象类别或创建对象clc 清除指令窗clear 清除内存变量和函数clf 清除图对象clock 时钟colorcube 三浓淡多彩交叉色图矩阵colordef 设置色彩缺省值colormap 色图colspace 列空间的基close 关闭指定窗口colperm 列排序置换向量comet 彗星状轨迹图comet3 三维彗星轨迹图compass 射线图compose 求复合函数cond (逆)条件数condeig 计算特征值、特征向量同时给出条件数condest 范-1条件数估计conj 复数共轭contour 等位线填色等位线contourfcontour3 三维等位线contourslice 四维切片等位线图conv 多项式乘、卷积cool 青紫调冷色图copper 古铜调色图cos 余弦cosh 双曲余弦cot 余切coth 双曲余切cplxpair 复数共轭成对排列csc 余割csch 双曲余割cumsum 元素累计和cumtrapz 累计梯形积分cylinder 创建圆柱D ddblquad 二重数值积分deal 分配宗量deblank 删去串尾部的空格符dec2base 十进制转换为X进制dec2bin 十进制转换为二进制dec2hex 十进制转换为十六进制deconv 多项式除、解卷delaunay Delaunay 三角剖分del2 离散Laplacian差分demo Matlab演示det 行列式diag 矩阵对角元素提取、创建对角阵diary Matlab指令窗文本内容记录diff 数值差分、符号微分digits 符号计算中设置符号数值的精度dir 目录列表disp 显示数组display 显示对象内容的重载函数dlinmod 离散系统的线性化模型dmperm 矩阵Dulmage-Mendelsohn 分解dos 执行DOS 指令并返回结果double 把其他类型对象转换为双精度数值drawnow 更新事件队列强迫Matlab刷新屏幕dsolve 符号计算解微分方程E eecho M文件被执行指令的显示edit 启动M文件编辑器eig 求特征值和特征向量eigs 求指定的几个特征值end 控制流FOR等结构体的结尾元素下标eps 浮点相对精度error 显示出错信息并中断执行errortrap 错误发生后程序是否继续执行的控制erf 误差函数erfc 误差补函数erfcx 刻度误差补函数erfinv 逆误差函数errorbar 带误差限的曲线图etreeplot 画消去树eval 串演算指令evalin 跨空间串演算指令exist 检查变量或函数是否已定义exit 退出Matlab环境exp 指数函数expand 符号计算中的展开操作expint 指数积分函数expm 常用矩阵指数函数expm1 Pade法求矩阵指数expm2 Taylor法求矩阵指数expm3 特征值分解法求矩阵指数eye 单位阵ezcontour 画等位线的简捷指令ezcontourf 画填色等位线的简捷指令ezgraph3 画表面图的通用简捷指令ezmesh 画网线图的简捷指令ezmeshc 画带等位线的网线图的简捷指令ezplot 画二维曲线的简捷指令ezplot3 画三维曲线的简捷指令ezpolar 画极坐标图的简捷指令ezsurf 画表面图的简捷指令ezsurfc 画带等位线的表面图的简捷指令F ffactor 符号计算的因式分解feather 羽毛图feedback 反馈连接feval 执行由串指定的函数fft 离散Fourier变换fft2 二维离散Fourier变换fftn 高维离散Fourier变换fftshift 直流分量对中的谱fieldnames 构架域名figure 创建图形窗fill3 三维多边形填色图find 寻找非零元素下标findobj 寻找具有指定属性的对象图柄findstr 寻找短串的起始字符下标findsym 机器确定内存中的符号变量finverse 符号计算中求反函数fix 向零取整flag 红白蓝黑交错色图阵fliplr 矩阵的左右翻转flipud 矩阵的上下翻转flipdim 矩阵沿指定维翻转floor 向负无穷取整flops 浮点运算次数flow Matlab提供的演示数据fmin 求单变量非线性函数极小值点(旧版)fminbnd 求单变量非线性函数极小值点fmins 单纯形法求多变量函数极小值点(旧版)fminunc 拟牛顿法求多变量函数极小值点fminsearch 单纯形法求多变量函数极小值点fnder 对样条函数求导fnint 利用样条函数求积分fnval 计算样条函数区间内任意一点的值fnplt 绘制样条函数图形fopen 打开外部文件for 构成for环用format 设置输出格式fourier Fourier 变换fplot 返函绘图指令fprintf 设置显示格式fread 从文件读二进制数据fsolve 求多元函数的零点full 把稀疏矩阵转换为非稀疏阵funm 计算一般矩阵函数函数计算器图形用户界面funtool fzero 求单变量非线性函数的零点G ggamma 函数gammainc 不完全函数gammaln 函数的对数gca 获得当前轴句柄gcbo 获得正执行回调的对象句柄gcf 获得当前图对象句柄gco 获得当前对象句柄geomean 几何平均值get 获知对象属性getfield 获知构架数组的域getframe 获取影片的帧画面ginput 从图形窗获取数据global 定义全局变量gplot 依图论法则画图gradient 近似梯度gray 黑白灰度grid 画分格线griddata 规则化数据和曲面拟合gtext 由鼠标放置注释文字guide 启动图形用户界面交互设计工具H hharmmean 调和平均值help 在线帮助helpwin 交互式在线帮助helpdesk 打开超文本形式用户指南hex2dec 十六进制转换为十进制hex2num 十六进制转换为浮点数hidden 透视和消隐开关hilb Hilbert矩阵hist 频数计算或频数直方图histc 端点定位频数直方图histfit 带正态拟合的频数直方图hold 当前图上重画的切换开关horner 分解成嵌套形式hot 黑红黄白色图饱和色图hsvI iif-else-elseif 条件分支结构ifft 离散Fourier反变换ifft2 二维离散Fourier反变换ifftn 高维离散Fourier反变换ifftshift 直流分量对中的谱的反操作ifourier Fourier反变换i, j 缺省的虚单元变量ilaplace Laplace反变换imag 复数虚部image 显示图象imagesc 显示亮度图象imfinfo 获取图形文件信息imread 从文件读取图象imwrite 把imwrite 把图象写成文件ind2sub 单下标转变为多下标inf 无穷大info MathWorks公司网点地址inline 构造内联函数对象inmem 列出内存中的函数名input 提示用户输入inputname 输入宗量名int 符号积分int2str 把整数数组转换为串数组interp1 一维插值interp2 二维插值interp3 三维插值interpn N维插值interpft 利用FFT插值intro Matlab自带的入门引导inv 求矩阵逆invhilb Hilbert矩阵的准确逆ipermute 广义反转置isa 检测是否给定类的对象ischar 若是字符串则为真isequal 若两数组相同则为真isempty 若是空阵则为真isfinite 若全部元素都有限则为真isfield 若是构架域则为真若是全局变量则为真isglobal ishandle 若是图形句柄则为真ishold 若当前图形处于保留状态则为真isieee 若计算机执行IEEE规则则为真isinf 若是无穷数据则为真isletter 若是英文字母则为真islogical 若是逻辑数组则为真ismember 检查是否属于指定集isnan 若是非数则为真isnumeric 若是数值数组则为真isobject 若是对象则为真isprime 若是质数则为真isreal 若是实数则为真isspace 若是空格则为真issparse 若是稀疏矩阵则为真isstruct 若是构架则为真isstudent 若是Matlab学生版则为真iztrans 符号计算Z反变换J j , K kjacobian 符号计算中求Jacobian 矩阵jet 蓝头红尾饱和色jordan 符号计算中获得Jordan标准型keyboard 键盘获得控制权kron Kronecker乘法规则产生的数组L llaplace Laplace变换lasterr 显示最新出错信息lastwarn 显示最新警告信息leastsq 解非线性最小二乘问题(旧版)legend 图形图例lighting 照明模式line 创建线对象lines 采用plot 画线色linmod 获连续系统的线性化模型linmod2 获连续系统的线性化精良模型linspace 线性等分向量ln 矩阵自然对数load 从MAT文件读取变量自然对数loglog10 常用对数log2 底为2的对数loglog 双对数刻度图形logm 矩阵对数logspace 对数分度向量lookfor 按关键字搜索M文件lower 转换为小写字母lsqnonlin 解非线性最小二乘问题lu LU分解M mmad 平均绝对值偏差magic 魔方阵maple &nb, sp; 运作Maple格式指令mat2str 把数值数组转换成输入形态串数组material 材料反射模式max 找向量中最大元素mbuild 产生EXE文件编译环境的预设置指令mcc 创建MEX或EXE文件的编译指令mean 求向量元素的平均值median 求中位数menuedit 启动设计用户菜单的交互式编辑工具mesh 网线图meshz 垂帘网线图meshgrid 产生格点矩阵methods 获知对指定类定义的所有方法函数mex 产生MEX文件编译环境的预设置指令mfunlis 能被mfun计算的MAPLE经典函数列表mhelp 引出Maple的在线帮助min 找向量中最小元素mkdir 创建目录mkpp 逐段多项式数据的明晰化mod 模运算more 指令窗中内容的分页显示movie 放映影片动画moviein 影片帧画面的内存预置mtaylor 符号计算多变量Taylor级数展开N n求数组维数ndimsNaN 非数(预定义)变量nargchk 输入宗量数验证nargin 函数输入宗量数nargout 函数输出宗量数ndgrid 产生高维格点矩阵newplot 准备新的缺省图、轴nextpow2 取最接近的较大2次幂nnz 矩阵的非零元素总数nonzeros 矩阵的非零元素norm 矩阵或向量范数normcdf 正态分布累计概率密度函数normest 估计矩阵2范数norminv 正态分布逆累计概率密度函数normpdf 正态分布概率密度函数normrnd 正态随机数发生器notebook 启动Matlab和Word的集成环境null 零空间num2str 把非整数数组转换为串numden 获取最小公分母和相应的分子表达式nzmax 指定存放非零元素所需内存O oode1 非Stiff 微分方程变步长解算器ode15s Stiff 微分方程变步长解算器ode23t 适度Stiff 微分方程解算器ode23tb Stiff 微分方程解算器ode45 非Stiff 微分方程变步长解算器odefile ODE 文件模板odeget 获知ODE 选项设置参数odephas2 ODE 输出函数的二维相平面图odephas3 ODE 输出函数的三维相空间图odeplot ODE 输出函数的时间轨迹图odeprint 在Matlab指令窗显示结果odeset 创建或改写ODE选项构架参数值ones 全1数组optimset 创建或改写优化泛函指令的选项参数值orient 设定图形的排放方式orth 值空间正交化P ppack 收集Matlab内存碎块扩大内存pagedlg 调出图形排版对话框patch 创建块对象path 设置Matlab搜索路径的指令pathtool 搜索路径管理器pause 暂停pcode 创建预解译P码文件pcolor 伪彩图peaks Matlab提供的典型三维曲面permute 广义转置pi (预定义变量)圆周率pie 二维饼图pie3 三维饼图pink 粉红色图矩阵pinv 伪逆plot 平面线图plot3 三维线图plotmatrix 矩阵的散点图plotyy 双纵坐标图poissinv 泊松分布逆累计概率分布函数poissrnd 泊松分布随机数发生器pol2cart 极或柱坐标变为直角坐标polar 极坐标图poly 矩阵的特征多项式、根集对应的多项式poly2str 以习惯方式显示多项式poly2sym 双精度多项式系数转变为向量符号多项式polyder 多项式导数polyfit 数据的多项式拟合polyval 计算多项式的值polyvalm 计算矩阵多项式pow2 2的幂ppval 计算分段多项式pretty 以习惯方式显示符号表达式print 打印图形或SIMULINK模型printsys 以习惯方式显示有理分式prism 光谱色图矩阵procread 向MAPLE输送计算程序profile 函数文件性能评估器propedit 图形对象属性编辑器pwd 显示当前工作目录Q qquad 低阶法计算数值积分quad8 高阶法计算数值积分(QUADL)quit 推出Matlab 环境quiver 二维方向箭头图quiver3 三维方向箭头图R rrand 产生均匀分布随机数randn 产生正态分布随机数randperm 随机置换向量range 样本极差rank 矩阵的秩rats 有理输出rcond 矩阵倒条件数估计real 复数的实部reallog 在实数域内计算自然对数realpow 在实数域内计算乘方realsqrt 在实数域内计算平方根realmax 最大正浮点数realmin 最小正浮点数rectangle 画长方框rem 求余数repmat 铺放模块数组reshape 改变数组维数、大小residue 部分分式展开return 返回ribbon 把二维曲线画成三维彩带图rmfield 删去构架的域roots 求多项式的根rose 数扇形图rot90 矩阵旋转90度rotate 指定的原点和方向旋转rotate3d 启动三维图形视角的交互设置功能round 向最近整数圆整rref 简化矩阵为梯形形式rsf2csf 实数块对角阵转为复数特征值对角阵rsums Riemann和S ssave 把内存变量保存为文件scatter3 三维散点图sec 正割sech 双曲正割semilogx X轴对数刻度坐标图semilogy Y轴对数刻度坐标图series 串联连接set 设置图形对象属性setfield 设置构架数组的域setstr 将ASCII码转换为字符的旧版指令sign 根据符号取值函数signum 符号计算中的符号取值函数sim 运行SIMULINK模型simget 获取SIMULINK模型设置的仿真参数simple 寻找最短形式的符号解simplify 符号计算中进行简化操作simset 对SIMULINK模型的仿真参数进行设置simulink 启动SIMULINK模块库浏览器sin 正弦sinh 双曲正弦size 矩阵的大小slice 立体切片图solve 求代数方程的符号解spalloc 为非零元素配置内存sparse 创建稀疏矩阵spconvert 把外部数据转换为稀疏矩阵spdiags 稀疏对角阵spfun 求非零元素的函数值sph2cart 球坐标变为直角坐标sphere 产生球面spinmap 色图彩色的周期变化spline 样条插值spones 用1置换非零元素sprandsym 稀疏随机对称阵sprank 结构秩spring 紫黄调春色图sprintf 把格式数据写成串spy 画稀疏结构图sqrt 平方根sqrtm 方根矩阵squeeze 删去大小为1的孤维sscanf 按指定格式读串stairs 阶梯图std 标准差step 阶跃响应指令str2double 串转换为双精度值str2mat 创建多行串数组str2num 串转换为数strcat 接成长串strcmp 串比较strjust 串对齐strmatch 搜索指定串strncmp 串中前若干字符比较strrep 串替换strtok 寻找第一间隔符前的内容struct 创建构架数组struct2cell 把构架转换为元胞数组strvcat 创建多行串数组sub2ind 多下标转换为单下标subexpr 通过子表达式重写符号对象subplot 创建子图subs 符号计算中的符号变量置换subspace 两子空间夹角sum 元素和summer 绿黄调夏色图superiorto 设定优先级surf 三维着色表面图surface 创建面对象surfc 带等位线的表面图surfl 带光照的三维表面图surfnorm 空间表面的法线svd 奇异值分解svds 求指定的若干奇异值switch-case-otherwise 多分支结构sym2poly 符号多项式转变为双精度多项式系数向量symmmd 对称最小度排序symrcm 反向Cuthill-McKee排序syms 创建多个符号对象T ttan 正切tanh 双曲正切taylortool 进行Taylor逼近分析的交互界面text 文字注释tf 创建传递函数对象启动计时器tictitle 图名toc 关闭计时器trapz 梯形法数值积分treelayout 展开树、林treeplot 画树图tril 下三角阵trim 求系统平衡点trimesh 不规则格点网线图trisurf 不规则格点表面图triu 上三角阵try-catch 控制流中的Try-catch结构type 显示M文件U uuicontextmenu 创建现场菜单uicontrol 创建用户控件uimenu 创建用户菜单unmkpp 逐段多项式数据的反明晰化unwrap 自然态相角upper 转换为大写字母V vvar 方差varargin 变长度输入宗量varargout 变长度输出宗量vectorize 使串表达式或内联函数适于数组运算ver 版本信息的获取view 三维图形的视角控制voronoi V oronoi多边形vpa 任意精度(符号类)数值W wwarning 显示警告信息what 列出当前目录上的文件whatsnew 显示Matlab中Readme文件的内容which 确定函数、文件的位置while 控制流中的While环结构white 全白色图矩阵whitebg 指定轴的背景色who 列出内存中的变量名whos 列出内存中变量的详细信息winter 蓝绿调冬色图启动内存浏览器workspaceX x , Y y , Z zxlabel X轴名xor 或非逻辑yesinput 智能输入指令ylabel Y轴名zeros 全零数组zlabel Z轴名zoom 图形的变焦放大和缩小ztrans 符号计算Z变换MATLAB命令大全管理命令和函在线帮助文hel装入超文本说do文件的目录列MEwhaM文列typ条目搜索关键lookfo通hel定位函数和文whic运行演示程Dem的搜索路MATLA控Path管理变量和工作空间列出当前变量Who列出当前变量(长表)Whos从磁盘文件中恢复变量Load保存工作空间变量Save从内存中清除变量和函数Clear整理工作空间内存Pack矩阵的尺寸Size向量的长度Length显示矩阵或disp与文件和操作系统有关的命改变当前工作目c目录列Di删除文Delet获取环境变量Geten操作系统命执DO操作系统命令并返回结Uni执UNI任保DiarMATLA控制命令窗设置命令行编Cedi 清命令窗Cl光标置左上Hom设置输出格Forma底稿文件内使用的回显命Ech在命令窗口中控制分页输morTLAM启动和退MATLA退Qui文时所执行Startu引MATLA文Matlabr主启一般信公司信MathworkMATLA系统信息Inf的订购用成MATLASubscribhostiMA TLA主服务程序的识别代Whatsne在说明书中未包含的新信Ve版本信操作符和特殊字矩阵乘*数组乘.矩阵数组.左除或反斜右除或斜数组.张量KroneckeKro冒圆括( 方括[ 小数父目.继逗分注感叹转置或引赋相= 关系操作< 逻辑逻辑逻辑逻辑异xo逻辑函检查变量或函数是否存Exis向量的任一元为真,则其值为An向量的所有元为真,则其值为All找出非零元素的索引Fin三角函正Si双曲正Sin反正Asi反双曲正AsinCo余Cos双曲余Aco反余Acos反双曲余Ta正Tan双曲正Ata反正Atan四象限反正Atan反双曲正Se正Sec双曲正Asech 反双曲正Cs余Csc双曲余Acs反余Acsc反双曲余Co余Cot双曲余Aco反余Acot反双曲余指数函指Ex自然对Log常用对Log1平方Sqr复数函绝对Ab相Argl复共Con复数虚Imag复数实Rea数值函朝零方向取Fi朝负无穷大方向取Floo朝正无穷大方向取Cei朝最近的整数取Roun除后取Re符号函Sig 基本矩零矩Zero”矩全One单位矩EyRan均匀分布的随机数矩Rand正态分布的随机数矩Logspac 对数间隔的向Meshgri三维图形数规则间隔的向特殊变量和常当前的答An相对浮点精Ep最大浮点Realma最小浮点Realmin圆周P虚数单I,无穷In非数Na浮点运算次Flop函数输入变量Nargi函数输出变量Nargou计算机类Compute算术标准时,其值为IsieeIEE当计算机采简明的答Wh版本VersioMATLA时间和日挂Cloc日Dat计时函Etim秒表开始计Ti计时函To时间(以秒为单位CputimCP矩阵操建立和提取对角Dia矩阵作左右翻Flipl矩阵作上下翻Flipu改变矩阵大Reshap矩阵旋9Rot9提取矩阵的下三角部Tri提取矩阵的上三角部Tri矩阵的索引号,重新排列矩友矩Compa矩HadamarHadamard矩HankeHanke矩HilHilber矩InvhilHilberKroKronecke张量Magi魔方矩ToeplitToeplit矩VandeVandermond矩矩阵分计算矩阵条件Con计算矩阵或向量范Nor逆条件值估Rcond Linpac 计算矩阵Ran计算矩阵行列式De计算矩阵的Trac零矩Nul正交Ort线性方线性方程求分CholeskCho高斯消元法求系数L矩阵求In分解正交三角矩阵分解QQ矩阵伪Pin特征值和奇异求特征值和特征向Ei求特征多项Pol形HessberHes广义特征Q变复对角矩阵为实分块对角形Cdf2rdf分SchuSchu矩阵均衡处理以提高特征值精Balanc奇异值分Svd矩阵函矩阵指Exp文实Expmexp 通过泰勒级数求矩阵指Expm通过特征值和特征向量求矩阵指Expm矩阵对Log矩阵开平方Sqrt 一般矩阵的计Fun泛函——非线性数值方低阶法求解常微分方Ode2低阶法求解常微分方程并绘出结果图Ode23高阶法求解常微分方Ode4低阶法计算数值积QuaQuad高阶法计算数值积Fmi 单变量函数的极小变Fmin多变量函数的极小Fzer找出单变量函数的零Fplo函数绘多项式函求多项式Root构造具有指定根的多项Pol带矩阵变量的多项式计Polyval部分分式展开(留数计算Residu数据的多项式拟Polyfi微分多项Polyde多项式乘Conv多项式除Decon建立和控制图形窗建立图Figur获取当前图形的句Gc清除当前图Cl关闭图Clos 建立和控制坐标在标定位置上建立坐标Subplo在任意位置上建立坐标Axe获取当前坐标系的句Gc清除当前坐标Cl控制坐标系的刻度和形Axi控制伪彩色坐标刻Caxi保持当前图Hol句柄图形对建立图形窗Figur建立坐标Axe建立曲Lin建立文本Tex建立图形填充Patc建立曲Surfac 建立图Imag建立用户界面控Uicontro建立用户界面菜Uime句柄图形操设置对Se获取对象特Ge重置对象特Rese删除对Delete文性质Newplo预nextplo获取当前对象的句Gc填充未完成绘图事Drawno寻找指定特征值的对Findob打印和存打印图形或保存图Prin配置本地打印机缺省Printop设置纸张取Orien屏幕抓取当前图Captur图基线性图Plo对数坐标图Loglo轴为对数坐标Semilog半对数坐标图形轴为对数坐标半对数坐标图形Semilog绘制二维多边形填充Fil图特极坐标Pola条形Ba离散序列图或杆Ste阶梯Stair误差条Errorba直方His角度直方Ros区域Compas箭头Feathe绘图函Fplo星点Come图形注图形标Title轴标Xlabe轴标Ylabe文本注Tex用鼠标放置文Gtex网格Gri编程语MATLA增加新的函Functio 表达式构成的字TLA执行EvaM执行由字串指定的函Feva定义全局变Globa程序控制条件执行语I命令配合使Elsi命令配合使Elseii语句的结EniFor,whil重复执行指定次数(循环Fo重复执行不定次数(循环Whil终止循环的执Brea返回引用的函Retur显示信息并终止函数的执Erro交互输提示用户输Inpu像底稿文件一样使用键盘输Keyboar产生由用户输入选择的菜Men等待用户响Paus建立用户界面菜Uimen建立用户界面控Uicontro一般字符串函中有关字符串函数的说MATLAStrings变字符串为数Ab变数值为字符Setst当变量为字符串时其值为Isst空Blank删除尾部的空Deblan 从各个字符串中形成文本矩Str2ma表达式组成的TLAEva执行M字符串比比较字符Strcm在一字符串中查找另一个子Findst变字符串为大Uppe变字符串为小Lowe当变量为字母时,其值为Islette当变量为空白字符时,其值为Isspac字符串与数值之间变变数值为字符Num2st变整数为字符Int2st变字符串为数Str2nu变数值为格式控制下的字符Sprint变字符串为格式控制下的数Sscan十进制与十六进制数之间变标准下的浮点变十六进制Hex2nuIEE变十六制数为十进制Hex2de变十进制数为十六进制Dec2he建追加系统动态特Appen变量状态作为输Augstat从方框图中构造状态空间系Blkbuild系统的闭Cloo方框图建Connec两个多项式的卷Con从增益矩阵中形成离散状态估计Desti从增益矩阵中形成离散控制器和估计Dre产生随机离散模Drmode从增益矩阵中形成连续状态估计Esti反馈系统连FeedbacOrd产生二阶系统近时延PadPad并行系统连Paralle从增益矩阵中形成连续控制器和估计Re产生随机连续模Rmode串行系统连Serie从模型中删除输入、输出或状Ssdelet 从大系统中选择子系ssselec模型变变连续系统为离散系C2利用指定方法变连续为离散系C2d 带一延时变连续为离散系C2d变离散为连续系D2利用指定方法变离散为连续系D2c变根值表示为多项式表Pol部分分式展Residu变状态空间表示为传递函数表Ss2t变状态空间表示为零极点表Ss2z变传递函数表示为状态空间表Tf2s变传递函数表示为零极点表Tf2z变零极点表示为传递函数表Zp2tf变零极点表示为状态空间表Zp2s模型简平衡实Balrea离散平衡实Dbalrea离散模型降Dmodre 最小实现和零极点对Minrea模型降Modre模型实正则形Cano可控阶梯Ctrb可观阶梯Obsv采用相似变Ss2s模型特相对于白噪声的连续协方差响Cova可控性矩Ctr阻尼系数和固有频Dam 连续稳态(直流)增Dcgai相对于白噪声的离散协方差响Dcova离散阻尼系数和固有频Ddam离散系统增Ddcgai离散可控性和可观Dgra按幅值排序离散特征Dsor特征值和特征向Ei按实部排列连续特征Esor可控性和可观Gra可观性矩Obs按格式显示系Printsy多项式之Root传递零Tzero利用随机扰动法传递零Tzero时域响离散时间单位冲激响Dimpuls离散时间零输入响Dinitia任意输入下的离散时间仿Dlsi离散时间阶跃响Dste变换仿Filte单输入单输冲激响Impuls连续时间零输入响Initia任意输入下的连续时间仿Lsi低级时间响应函Ltit阶跃响Ste阶跃函Stepfu频域响图(频域响应BodBod离DbodBod离DnicholNicholNyquisDnyquis离离散奇异值频域DsigmBodFbod连续系统的快拉普拉斯变换频率响Freq变换频率响FreqLtif低级频率响应函Prod 计算数组元素的连乘积。

Matlab振动程序-代码作业

Matlab振动程序-代码作业

一、课题任务要求随着机械工业不断向自动化、高精度、智能化等方向的发展,在机械设备运行及生产过程中进行参量测试、分析与诊断等处理过程已成为必要环节,许多信号处理方法如时域统计分析、相关分析、相干分析、频谱分析等已经被广泛被应用与机械工程测试领域。

本文为机械测试信号的时域和频域分析,其中时域分析包括对信号最大值、最小值、中值、方差的分析,频域分析包括FFT分析、频谱分析、倒谱分析。

在进行上述分析之前先要对振动信号进行拟合。

机械振动分为确定性振动和随机振动,确定性振动又分为周期振动和非周期振动,周期振动又进一步分为简谐振动和复杂的周期振动。

所以可以根据上述的分类来拟合振动信号。

在设计信号的处理程序时,用MATLAB中的相关函数来对所拟合的振动信号进行时域分析和频域分析,并对绘出的频谱图进行说明。

二、技术路线对机械振动信号的时域和频域采集,根据所拟合的振动信号,选取所需要的时域性能指标和频域分析的性能指标对振动信号进行分析。

其中时域分析包括对信号最大值、最小值、中值、方差的分析,频域分析包括FFT分析、频谱分析、倒谱分析。

现构造一个振动信号(在该程序中以两个衰减振动分量和一个随机数rand 之和来拟合振动信号),再利用MATLAB 中的函数mean ()、min ()、max ()、std ()对离散序列中的平均值、最大值、最小值、标准差等时域性能进行分析,通过调用函数fft (y );psd (y );rcep (y )对该振动信号进行频域内的性能分析。

在设计过程中的理论知识有离散傅立叶变换(DFT )、功率谱的概念和意义以及倒谱的概念和意义。

①DFT 的定义和意义:DFT 的定义式为:DFT 的意义:DFT 的意义在于它表示信号中的各个频率的分量的振动幅值的大小,亦即该分量对于振动信号影响的大小。

通过DFT 的快速算法FFT 可以很方便的将振动信号的各个分量的幅值比重计算出来。

即可以把信号的主频分量提取出来。

微环谐振器

微环谐振器
A1 A2
S
B 1 t j B 2 j A1 t A2
B1 B1
B2 A2 B1 A1
散射矩阵法
B 1 t B j 2 j A1 t A2
B2 A2 B1 A1
0 exp[ j( R ) 2R ] A2 q 0 exp[ j( R )2R ] B 2
单环双波导
A4 A1 p 2 q p1 B 4 B1
A1
B1 A2 B2 A3 B3 B1 A4
并联双环滤波器
A41 A11 p 2 q p1 B 41 B 11 A42 A12 p 2 q p1 B 42 B 12
A12 B11 exp[ j( R )L]
A41 B 42 exp[ j( R )L]
1
ni (ni i ) i (ni j ) j ni i sin i ni j cos i

E r E r 0e i[k ( •sin r x cos r y )t ]
1
Et Et 0e i[k ( •sin t x cos t y )t ]
B11 A11 A21 A12 B12 A22 B22 A32 B32 A31 A42 B42
L
B21
B31
B41
A41
并联双环滤波器
B 11 MA11 NA41 B 41 NA11 MA41
2
B 12 MA12 NA42 B 42 NA12 MA42
B11 A11
i ( H ) E

Matlab源程序代码

Matlab源程序代码

Matlab源程序代码正弦波的源程序:(一),用到的函数1,f2t函数function x=f2t(X)global dt df t f T N%x=f2t(X)%x为时域的取样值矢量%X为x的傅氏变换%X与x长度相同并为2的整幂%本函数需要一个全局变量dt(时域取样间隔) X=[X(N/2+1:N),X(1:N/2)];x=ifft(X)/dt;end2,t2f函数。

function X=t2f(x)global dt df N t f T%X=t2f(x)%x为时域的取样值矢量%X为x的傅氏变换%X与x长度相同,并为2的整幂。

%本函数需要一个全局变量dt(时域取样间隔) H=fft(x);X=[H(N/2+1:N),H(1:N/2)]*dt;end(二),主程序。

1,%(1)绘出正弦信号波形及频谱global dt df t f Nclose allk=input('取样点数=2^k, k取10左右');if isempty(k), k=10; endf0=input('f0=取1(kz)左右');if isempty(f0), f0=1; endN=2^k;dt=0.01; %msdf=1/(N*dt); %KHzT=N*dt; %截短时间Bs=N*df/2; %系统带宽f=[-Bs+df/2:df:Bs]; %频域横坐标t=[-T/2+dt/2:dt:T/2]; %时域横坐标s=sin(2*pi*f0*t); %输入的正弦信号S=t2f(s); %S是s的傅氏变换a=f2t(S); %a是S的傅氏反变换a=real(a);as=abs(S);subplot(2,1,1) %输出的频谱plot(f,as,'b');gridaxis([-2*f0,+2*f0,min(as),max(as)]) xlabel('f (KHz)') ylabel('|S(f)| (V/KHz)') %figure(2) subplot(2,1,2) plot(t,a,'black') %输出信号波形画图gridaxis([-2/f0,+2/f0,-1.5,1.5])xlabel('t(ms)')ylabel('a(t)(V)')gtext('频谱图')最佳基带系统的源程序:(一),用到的函数f2t函数和t2f函数。

matlab有限元法计算谐振腔模式

matlab有限元法计算谐振腔模式

matlab有限元法计算谐振腔模式摘要:一、引言二、MATLAB有限元法基本原理1.有限元法简介2.MATLAB在有限元法中的应用三、谐振腔基本概念1.谐振腔定义2.谐振腔模式及其特点四、MATLAB计算谐振腔模式的步骤1.准备数据2.编写MATLAB代码3.运行代码及结果分析五、实例演示1.实例描述2.编写与运行MATLAB代码3.结果分析与讨论六、结论与展望正文:一、引言随着科技的飞速发展,微波技术在各个领域得到了广泛应用。

微波谐振腔作为微波技术中的基本元件,其内部的模式特性研究具有重要意义。

MATLAB 作为一种强大的科学计算软件,可以方便地应用于有限元法的计算。

本文将介绍利用MATLAB有限元法计算谐振腔模式的基本原理和方法,并通过实例进行详细说明。

二、MATLAB有限元法基本原理1.有限元法简介有限元法(Finite Element Method,简称FEM)是一种求解偏微分方程(组)的数值方法。

它将求解区域划分为若干个小的子区域,在每个子区域里用一个简单的函数来表示未知函数,然后通过求解离散化的代数方程组得到近似解。

2.MATLAB在有限元法中的应用MATLAB提供了丰富的有限元工具箱(PDE Toolbox),可以方便地应用于求解偏微分方程(组)。

利用MATLAB有限元法计算谐振腔模式,主要包括以下几个步骤:准备数据、编写MATLAB代码、运行代码及结果分析。

三、谐振腔基本概念1.谐振腔定义谐振腔是指一种封闭的空间,在其中电磁波能够维持稳定的振荡。

谐振腔内的电磁波满足特定的边界条件,这些条件决定了谐振腔内的模式。

2.谐振腔模式及其特点谐振腔模式是指在谐振腔内满足边界条件的电磁波振动形式。

谐振腔模式的特性包括:纵向对称性、横向分布、模态参量等。

不同模式的电磁波在谐振腔内的分布特点不同,因此研究谐振腔模式对于设计微波器件具有重要意义。

四、MATLAB计算谐振腔模式的步骤1.准备数据根据谐振腔的物理参数(如尺寸、边界条件等),准备输入数据。

MATLAB 光学谐振腔的设计

MATLAB 光学谐振腔的设计

A1:光学谐振腔的设计设计任务:有一平面镜和曲率半径为R的凹面镜,画出光束发散角与腔长L的关系曲线。

设计要求:(1)具有友好输入输出界面;(2)可输入不同凹面镜曲率半径值,查看结果变化情况。

function varargout = A1(varargin)% A1 M-file for A1.fig% A1, by itself, creates a new A1 or raises the existing% singleton*.%% H = A1 returns the handle to a new A1 or the handle to% the existing singleton*.%% A1('CALLBACK',hObject,eventData,handles,...) calls the local% function named CALLBACK in A1.M with the given input arguments.%% A1('Property','Value',...) creates a new A1 or raises the% existing singleton*. Starting from the left, property value pairs are% applied to the GUI before A1_OpeningFunction gets called. An% unrecognized property name or invalid value makes property application% stop. All inputs are passed to A1_OpeningFcn via varargin.%% *See GUI Options on GUIDE's Tools menu. Choose "GUI allows only one% instance to run (singleton)".%% See also: GUIDE, GUIDA TA, GUIHANDLES% Copyright 2002-2003 The MathWorks, Inc.% Edit the above text to modify the response to help A1% Last Modified by GUIDE v2.5 21-Oct-2010 18:17:18% Begin initialization code - DO NOT EDITgui_Singleton = 1;gui_State = struct('gui_Name', mfilename, ...'gui_Singleton', gui_Singleton, ...'gui_OpeningFcn', @A1_OpeningFcn, ...'gui_OutputFcn', @A1_OutputFcn, ...'gui_LayoutFcn', [] , ...'gui_Callback', []);if nargin && ischar(varargin{1})gui_State.gui_Callback = str2func(varargin{1});endif nargout[varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:}); elsegui_mainfcn(gui_State, varargin{:});end% End initialization code - DO NOT EDIT% --- Executes just before A1 is made visible.function A1_OpeningFcn(hObject, eventdata, handles, varargin)% This function has no output args, see OutputFcn.% hObject handle to figure% eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDA TA) % varargin command line arguments to A1 (see V ARARGIN)% Choose default command line output for A1handles.output = hObject;% Update handles structureguidata(hObject, handles);% UIWAIT makes A1 wait for user response (see UIRESUME)% uiwait(handles.figure1);% --- Outputs from this function are returned to the command line. function varargout = A1_OutputFcn(hObject, eventdata, handles)% varargout cell array for returning output args (see V ARARGOUT); % hObject handle to figure% eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDA TA)% Get default command line output from handles structure varargout{1} = handles.output;% --- Executes on button press in pushbutton1.function pushbutton1_Callback(hObject, eventdata, handles)% hObject handle to pushbutton1 (see GCBO)% handles structure with handles and user data (see GUIDA TA)handles.R=str2double(get(handles.editR,'string'));handles.B=str2double(get(handles.editB,'string'));guidata(hObject,handles);L=0.1*10^1:0.1:10^2;f=sqrt(L*handles.R*10^2-L.^2);w=sqrt(handles.B*10^(-4)*f/pi);jiao=2*handles.B*10^(-4)/pi./w;plot(L,jiao);grid on;xlabel('腔长l(cm)');ylabel('散射角(rad)');title('光束发散角与腔长的关系曲线');function editR_Callback(hObject, eventdata, handles)% hObject handle to editR (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDA TA)% Hints: get(hObject,'String') returns contents of editR as text% str2double(get(hObject,'String')) returns contents of editR as a double% --- Executes during object creation, after setting all properties.function editR_CreateFcn(hObject, eventdata, handles)% hObject handle to editR (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles empty - handles not created until after all CreateFcns called% Hint: edit controls usually have a white background on Windows.% See ISPC and COMPUTER.if ispcset(hObject,'BackgroundColor','white');elseset(hObject,'BackgroundColor',get(0,'defaultUicontrolBackgroundColor')); endfunction editB_Callback(hObject, eventdata, handles)% hObject handle to editB (see GCBO)% handles structure with handles and user data (see GUIDA TA)% Hints: get(hObject,'String') returns contents of editB as text% str2double(get(hObject,'String')) returns contents of editB as a double% --- Executes during object creation, after setting all properties.function editB_CreateFcn(hObject, eventdata, handles)% hObject handle to editB (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles empty - handles not created until after all CreateFcns called% Hint: edit controls usually have a white background on Windows.% See ISPC and COMPUTER.if ispcset(hObject,'BackgroundColor','white');elseset(hObject,'BackgroundColor',get(0,'defaultUicontrolBackgroundColor')); end% --- Executes on button press in pushbutton2.function pushbutton2_Callback(hObject, eventdata, handles)close(gcf)% hObject handle to pushbutton2 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDA TA)% --- Executes on key press over pushbutton2 with no controls selected. function pushbutton2_KeyPressFcn(hObject, eventdata, handles)% hObject handle to pushbutton2 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDA TA)% --- Executes on key press over pushbutton1 with no controls selected. function pushbutton1_KeyPressFcn(hObject, eventdata, handles)% hObject handle to pushbutton1 (see GCBO)% handles structure with handles and user data (see GUIDA TA)。

基于微环谐振腔产生克尔光频梳的关键理论研究

基于微环谐振腔产生克尔光频梳的关键理论研究

摘要摘要光学频率梳技术的出现,为实现光学频率的精准测量提供了可靠的解决方案。

近年来,一种基于微谐振腔的克尔光频梳,以其卓越的光频梳性能和简单的实验配置成为研究的热点。

自从在微环芯谐振腔中观察到克尔光频梳后,相关实验成果不断得到报道,相应的理论研究也逐渐得以完善。

但是,目前克尔光频梳的产生方案仍存在很多不足,从而限制了进一步对克尔光频梳的应用研究。

本文对基于微环谐振腔的克尔光频梳的产生原理与机制,性质以及相关应用展开了深入的理论研究,主要工作和创新如下:1.在微腔内由克尔非线性效应引起光学参量振荡的基础上,深入分析克尔光频梳在频域上的演变过程。

通过微环内光场传输方程建立描述微环内克尔光频梳产生过程的LLE模型,包括微环损耗、色散和非线性效应等。

利用MATLAB实现对克尔光频梳产生过程的数值模拟,根据其演变过程将其分为四个阶段,分析不同阶段中光频梳的表现形式和产生机制,并与近期报道的相关实验结果对比。

2.目前实现相位锁定光梳主要通过扫描泵浦光频率在微腔内构造耗散腔孤子。

在MATLAB中通过扫描泵浦失谐值,实现对耗散腔孤子产生过程的数值模拟。

根据仿真结果,将微环内耗散腔孤子的演变过程划分为四个阶段,分析不同阶段下相应光频梳的噪声特性和相干特性,揭示耗散腔孤子以及相应相位锁定光梳产生的根本机制。

此外,在对耗散腔孤子产生机制的深刻理解上,进一步通过降低泵浦光功率以减少微腔内的耗散腔孤子数量,从而改善克尔光频梳的平坦度。

3.由于克尔光频梳具有频谱覆盖范围广、重复频率高等优点,本文提出了一种基于克尔光频梳的新型光子信道化接收方案。

首先,在OptiSystem中建立基于克尔光频梳的信道化接收系统,通过实现对宽带RF信号的频率检测验证了该信道化接收方案的可行性;其次,通过实现对加载在宽带RF信号上的模拟基带信号的接收,进一步验证了本文提出的基于克尔光频梳的光子信道化接收系统具有良好的接收性能。

关键词:克尔光频梳,克尔非线性效应,LLE模型,耗散腔孤子,信道化ABSTRACTThe advent of optical frequency comb technology provides a reliable solution for accurate measurement of optical frequency and time. In recent years, Kerr frequency comb with excellent optical comb performance and simple experimental configuration based on micro-resonators has aroused general interests. Since the Kerr optical frequency combs have been observed in the micro-resonator, the theoretical basis has been improved a lot and the experiment has made great progress. However, there are still many shortcomings in the current Kerr optical frequency combs generation, which severe limits its further application. We are about carrying out the key theoretical research on the principle and mechanism of Kerr combs generation, and the characteristic and related applications of Kerr combs based on the micro-ring resonator. The main work and innovation are as follows:1.Based on the research of optical parametric oscillation caused by Kerr nonlinear effect in the micro-resonator, we analyze the evolution of Kerr optical comb in frequency domain. The LLE model based on the optical pulses transmission equation in the micro-resonator including the loss, dispersion and nonlinear effects, describe the evolution of Kerr frequency combs. The simulation of Kerr frequency combs is realized by MATLAB. It is divided into four stages according to the simulation results, and the manifestations and mechanisms of Kerr combs at different stages are analyzed and compared with the recent experimental results.2.By changing the detuning value of the pump C.W., the simulation of the dissipative cavity soliton’s evolution in the micro-ring resonator is realized. According to the simulation results, the evolution process of the micro-ring dissipative cavity soliton is divided into four stages, too. The noise characteristics and coherence characteristics of the corresponding optical combs at different stages are analyzed, and the fundamental mechanism of the dissipative cavity soliton and the corresponding phase-locked optical frequency combs are revealed. In addition, the flatness of the Kerr combs is improved by properly tuning the pumped optical power.3.In this paper, we propose a new photonic channelization receiver scheme based on Kerr's optical comb, which has the advantage of a wide spectrum coverage and high repetition frequency. Firstly, the channelization receiver system based on Kerr opticalfrequency comb is established in OptiSystem. The feasibility of the channelized reception scheme is verified by realizing the frequency detection of the specific broadband RF signal. Then, the good receiving performance has been proved by the result of detecting the microwave signal loaded on the broadband RF signal, which validates the photonic channelization receiver system based on Kerr combs further.Keywords:Kerr frequency combs, Kerr nonlinearity, LLE model, Dissipative cavity soliton, Channelization目录第一章绪论 (1)1.1引言 (1)1.2光频梳的研究背景 (1)1.3 克尔光频梳的研究现状 (4)1.4 本论文的主要内容和结构 (7)第二章微环谐振腔中产生克尔光频梳的理论基础 (9)2.1 微谐振腔中的非线性光学效应原理 (9)2.1.1 克尔非线性效应 (10)2.1.2 微谐振腔的热动态 (12)2.1.3 微谐振腔中的参量振荡与调制不稳定性 (13)2.2 克尔光频梳的形成过程分析 (15)2.3 微谐振腔中克尔光频梳的模型建立 (17)2.4 本章小结 (19)第三章克尔光频梳的仿真实现及其性质研究 (20)3.1 微环谐振腔中克尔光频梳的仿真实现 (20)3.1.1 LLE模型的仿真实现 (20)3.1.2 不同泵浦光条件下克尔光频梳的性质研究 (24)3.2 耗散腔孤子的形成及其性质研究 (27)3.2.1 耗散腔孤子形成过程的仿真实现 (27)3.2.2 耗散腔孤子的性质研究 (30)3.3 关于克尔光频梳平坦度的关键理论研究 (33)3.4 本章小结 (39)第四章克尔光频梳在光子信道化接收机中的应用研究 (40)4.1 基于光频梳的光子信道化接收方案 (40)4.2 基于克尔光频梳的新型光子信道化接收方案 (41)4.3本章小结 (51)第五章总结与展望 (52)5.1本文总结 (52)5.2工作展望 (52)致谢 (54)参考文献 (55)第一章 绪论第一章 绪论1.1引言随着信息化时代的快速发展,现代电磁环境环境变得越来越复杂,其所含信号种类繁多,带宽覆盖范围日益增大,并且同一时间段内存在多个信号[1]。

(完整word版)matlab函数大全(非常实用),推荐文档

(完整word版)matlab函数大全(非常实用),推荐文档

信源函数randerr 产生比特误差样本randint 产生均匀分布的随机整数矩阵randsrc 根据给定的数字表产生随机矩阵wgn 产生高斯白噪声信号分析函数biterr 计算比特误差数和比特误差率eyediagram 绘制眼图scatterplot 绘制分布图symerr 计算符号误差数和符号误差率信源编码compand mu律/A律压缩/扩张dpcmdeco DPCM(差分脉冲编码调制)解码dpcmenco DPCM编码dpcmopt 优化DPCM参数lloyds Lloyd法则优化量化器参数quantiz 给出量化后的级和输出值误差控制编码bchpoly 给出二进制BCH码的性能参数和产生多项式convenc 产生卷积码cyclgen 产生循环码的奇偶校验阵和生成矩阵cyclpoly 产生循环码的生成多项式decode 分组码解码器encode 分组码编码器gen2par 将奇偶校验阵和生成矩阵互相转换gfweight 计算线性分组码的最小距离hammgen 产生汉明码的奇偶校验阵和生成矩阵rsdecof 对Reed-Solomon编码的ASCII文件解码rsencof 用Reed-Solomon码对ASCII文件编码rspoly 给出Reed-Solomon码的生成多项式syndtable 产生伴随解码表vitdec 用Viterbi法则解卷积码(误差控制编码的低级函数)bchdeco BCH解码器bchenco BCH编码器rsdeco Reed-Solomon解码器rsdecode 用指数形式进行Reed-Solomon解码rsenco Reed-Solomon编码器rsencode 用指数形式进行Reed-Solomon编码调制与解调ademod 模拟通带解调器ademodce 模拟基带解调器amod 模拟通带调制器amodce 模拟基带调制器apkconst 绘制圆形的复合ASK-PSK星座图ddemod 数字通带解调器ddemodce 数字基带解调器demodmap 解调后的模拟信号星座图反映射到数字信号dmod 数字通带调制器dmodce 数字基带调制器modmap 把数字信号映射到模拟信号星座图(以供调制)qaskdeco 从方形的QASK星座图反映射到数字信号qaskenco 把数字信号映射到方形的QASK星座图专用滤波器hank2sys 把一个Hankel矩阵转换成一个线性系统模型hilbiir 设计一个希尔伯特变换IIR滤波器rcosflt 升余弦滤波器rcosine 设计一个升余弦滤波器(专用滤波器的低级函数)rcosfir 设计一个升余弦FIR滤波器rcosiir 设计一个升余弦IIR滤波器信道函数awgn 添加高斯白噪声伽罗域计算gfadd 伽罗域上的多项式加法gfconv 伽罗域上的多项式乘法gfcosets 生成伽罗域的分圆陪集gfdeconv 伽罗域上的多项式除法gfdiv 伽罗域上的元素除法gffilter 在质伽罗域上用多项式过滤数据gflineq 在至伽罗域上求Ax=b的一个特解gfminpol 求伽罗域上元素的最小多项式gfmul 伽罗域上的元素乘法gfplus GF(2^m)上的元素加法gfpretty 以通常方式显示多项式gfprimck 检测多项式是否是基本多项式gfprimdf 给出伽罗域的MATLAB默认的基本多项式gfprimfd 给出伽罗域的基本多项式gfrank 伽罗域上矩阵求秩gfrepcov GF(2)上多项式的表达方式转换gfroots 质伽罗域上的多项式求根gfsub 伽罗域上的多项式减法gftrunc 使多项式的表达最简化gftuple 简化或转换伽罗域上元素的形式工具函数bi2de 把二进制向量转换成十进制数de2bi 把十进制数转换成二进制向量erf 误差函数erfc 余误差函数istrellis 检测输入是否MA TLAB的trellis结构(structure)marcumq 通用Marcum Q 函数oct2dec 八进制数转十进制数poly2trellis 把卷积码多项式转换成MATLAB的trellis描述vec2mat 把向量转换成矩阵——————————————————————————————————————————————————A aabs 绝对值、模、字符的ASCII码值acos 反余弦acosh 反双曲余弦acot 反余切acoth 反双曲余切acsc 反余割acsch 反双曲余割align 启动图形对象几何位置排列工具all 所有元素非零为真angle 相角ans 表达式计算结果的缺省变量名any 所有元素非全零为真area 面域图argnames 函数M文件宗量名asec 反正割asech 反双曲正割asin 反正弦asinh 反双曲正弦assignin 向变量赋值atan 反正切atan2 四象限反正切atanh 反双曲正切autumn 红黄调秋色图阵axes 创建轴对象的低层指令axis 控制轴刻度和风格的高层指令B bbar 二维直方图bar3 三维直方图bar3h 三维水平直方图barh 二维水平直方图base2dec X进制转换为十进制bin2dec 二进制转换为十进制blanks 创建空格串bone 蓝色调黑白色图阵box 框状坐标轴break while 或for 环中断指令brighten 亮度控制C ccapture (3版以前)捕获当前图形cart2pol 直角坐标变为极或柱坐标cart2sph 直角坐标变为球坐标cat 串接成高维数组caxis 色标尺刻度cd 指定当前目录cdedit 启动用户菜单、控件回调函数设计工具cdf2rdf 复数特征值对角阵转为实数块对角阵ceil 向正无穷取整cell 创建元胞数组cell2struct 元胞数组转换为构架数组celldisp 显示元胞数组内容cellplot 元胞数组内部结构图示char 把数值、符号、内联类转换为字符对象chi2cdf 分布累计概率函数chi2inv 分布逆累计概率函数chi2pdf 分布概率密度函数chi2rnd 分布随机数发生器chol Cholesky分解clabel 等位线标识cla 清除当前轴class 获知对象类别或创建对象clc 清除指令窗clear 清除内存变量和函数clf 清除图对象clock 时钟colorcube 三浓淡多彩交叉色图矩阵colordef 设置色彩缺省值colormap 色图colspace 列空间的基close 关闭指定窗口colperm 列排序置换向量comet 彗星状轨迹图comet3 三维彗星轨迹图compass 射线图compose 求复合函数cond (逆)条件数condeig 计算特征值、特征向量同时给出条件数condest 范-1条件数估计conj 复数共轭contour 等位线contourf 填色等位线contour3 三维等位线contourslice 四维切片等位线图conv 多项式乘、卷积cool 青紫调冷色图copper 古铜调色图cos 余弦cosh 双曲余弦cot 余切coth 双曲余切cplxpair 复数共轭成对排列csc 余割csch 双曲余割cumsum 元素累计和cumtrapz 累计梯形积分cylinder 创建圆柱D ddblquad 二重数值积分deal 分配宗量deblank 删去串尾部的空格符dec2base 十进制转换为X进制dec2bin 十进制转换为二进制dec2hex 十进制转换为十六进制deconv 多项式除、解卷delaunay Delaunay 三角剖分del2 离散Laplacian差分demo Matlab演示det 行列式diag 矩阵对角元素提取、创建对角阵diary Matlab指令窗文本内容记录diff 数值差分、符号微分digits 符号计算中设置符号数值的精度dir 目录列表disp 显示数组display 显示对象内容的重载函数dlinmod 离散系统的线性化模型dmperm 矩阵Dulmage-Mendelsohn 分解dos 执行DOS 指令并返回结果double 把其他类型对象转换为双精度数值drawnow 更新事件队列强迫Matlab刷新屏幕dsolve 符号计算解微分方程E eecho M文件被执行指令的显示edit 启动M文件编辑器eig 求特征值和特征向量eigs 求指定的几个特征值end 控制流FOR等结构体的结尾元素下标eps 浮点相对精度error 显示出错信息并中断执行errortrap 错误发生后程序是否继续执行的控制erf 误差函数erfc 误差补函数erfcx 刻度误差补函数erfinv 逆误差函数errorbar 带误差限的曲线图etreeplot 画消去树eval 串演算指令evalin 跨空间串演算指令exist 检查变量或函数是否已定义exit 退出Matlab环境exp 指数函数expand 符号计算中的展开操作expint 指数积分函数expm 常用矩阵指数函数expm1 Pade法求矩阵指数expm2 Taylor法求矩阵指数expm3 特征值分解法求矩阵指数eye 单位阵ezcontour 画等位线的简捷指令ezcontourf 画填色等位线的简捷指令ezgraph3 画表面图的通用简捷指令ezmesh 画网线图的简捷指令ezmeshc 画带等位线的网线图的简捷指令ezplot 画二维曲线的简捷指令ezplot3 画三维曲线的简捷指令ezpolar 画极坐标图的简捷指令ezsurf 画表面图的简捷指令ezsurfc 画带等位线的表面图的简捷指令F ffactor 符号计算的因式分解feather 羽毛图feedback 反馈连接feval 执行由串指定的函数fft 离散Fourier变换fft2 二维离散Fourier变换fftn 高维离散Fourier变换fftshift 直流分量对中的谱fieldnames 构架域名figure 创建图形窗fill3 三维多边形填色图find 寻找非零元素下标findobj 寻找具有指定属性的对象图柄findstr 寻找短串的起始字符下标findsym 机器确定内存中的符号变量finverse 符号计算中求反函数fix 向零取整flag 红白蓝黑交错色图阵fliplr 矩阵的左右翻转flipud 矩阵的上下翻转flipdim 矩阵沿指定维翻转floor 向负无穷取整flops 浮点运算次数flow Matlab提供的演示数据fmin 求单变量非线性函数极小值点(旧版)fminbnd 求单变量非线性函数极小值点fmins 单纯形法求多变量函数极小值点(旧版)fminunc 拟牛顿法求多变量函数极小值点fminsearch 单纯形法求多变量函数极小值点fnder 对样条函数求导fnint 利用样条函数求积分fnval 计算样条函数区间内任意一点的值fnplt 绘制样条函数图形fopen 打开外部文件for 构成for环用format 设置输出格式fourier Fourier 变换fplot 返函绘图指令fprintf 设置显示格式fread 从文件读二进制数据fsolve 求多元函数的零点full 把稀疏矩阵转换为非稀疏阵funm 计算一般矩阵函数funtool 函数计算器图形用户界面fzero 求单变量非线性函数的零点G ggamma 函数gammainc 不完全函数gammaln 函数的对数gca 获得当前轴句柄gcbo 获得正执行"回调"的对象句柄gcf 获得当前图对象句柄gco 获得当前对象句柄geomean 几何平均值get 获知对象属性getfield 获知构架数组的域getframe 获取影片的帧画面ginput 从图形窗获取数据global 定义全局变量gplot 依图论法则画图gradient 近似梯度gray 黑白灰度grid 画分格线griddata 规则化数据和曲面拟合gtext 由鼠标放置注释文字guide 启动图形用户界面交互设计工具H hharmmean 调和平均值help 在线帮助helpwin 交互式在线帮助helpdesk 打开超文本形式用户指南hex2dec 十六进制转换为十进制hex2num 十六进制转换为浮点数hidden 透视和消隐开关hilb Hilbert矩阵hist 频数计算或频数直方图histc 端点定位频数直方图histfit 带正态拟合的频数直方图hold 当前图上重画的切换开关horner 分解成嵌套形式hot 黑红黄白色图hsv 饱和色图I iif-else-elseif 条件分支结构ifft 离散Fourier反变换ifft2 二维离散Fourier反变换ifftn 高维离散Fourier反变换ifftshift 直流分量对中的谱的反操作ifourier Fourier反变换i, j 缺省的"虚单元"变量ilaplace Laplace反变换imag 复数虚部image 显示图象imagesc 显示亮度图象imfinfo 获取图形文件信息imread 从文件读取图象imwrite 把imwrite 把图象写成文件ind2sub 单下标转变为多下标inf 无穷大info MathWorks公司网点地址inline 构造内联函数对象inmem 列出内存中的函数名input 提示用户输入inputname 输入宗量名int 符号积分int2str 把整数数组转换为串数组interp1 一维插值interp2 二维插值interp3 三维插值interpn N维插值interpft 利用FFT插值intro Matlab自带的入门引导inv 求矩阵逆invhilb Hilbert矩阵的准确逆ipermute 广义反转置isa 检测是否给定类的对象ischar 若是字符串则为真isequal 若两数组相同则为真isempty 若是空阵则为真isfinite 若全部元素都有限则为真isfield 若是构架域则为真isglobal 若是全局变量则为真ishandle 若是图形句柄则为真ishold 若当前图形处于保留状态则为真isieee 若计算机执行IEEE规则则为真isinf 若是无穷数据则为真isletter 若是英文字母则为真islogical 若是逻辑数组则为真ismember 检查是否属于指定集isnan 若是非数则为真isnumeric 若是数值数组则为真isobject 若是对象则为真isprime 若是质数则为真isreal 若是实数则为真isspace 若是空格则为真issparse 若是稀疏矩阵则为真isstruct 若是构架则为真isstudent 若是Matlab学生版则为真iztrans 符号计算Z反变换J j , K kjacobian 符号计算中求Jacobian 矩阵jet 蓝头红尾饱和色jordan 符号计算中获得Jordan标准型keyboard 键盘获得控制权kron Kronecker乘法规则产生的数组L llaplace Laplace变换lasterr 显示最新出错信息lastwarn 显示最新警告信息leastsq 解非线性最小二乘问题(旧版)legend 图形图例lighting 照明模式line 创建线对象lines 采用plot 画线色linmod 获连续系统的线性化模型linmod2 获连续系统的线性化精良模型linspace 线性等分向量ln 矩阵自然对数load 从MAT文件读取变量log 自然对数log10 常用对数log2 底为2的对数loglog 双对数刻度图形logm 矩阵对数logspace 对数分度向量lookfor 按关键字搜索M文件lower 转换为小写字母lsqnonlin 解非线性最小二乘问题lu LU分解M mmad 平均绝对值偏差magic 魔方阵maple &nb, sp; 运作Maple格式指令mat2str 把数值数组转换成输入形态串数组material 材料反射模式max 找向量中最大元素mbuild 产生EXE文件编译环境的预设置指令mcc 创建MEX或EXE文件的编译指令mean 求向量元素的平均值median 求中位数menuedit 启动设计用户菜单的交互式编辑工具mesh 网线图meshz 垂帘网线图meshgrid 产生"格点"矩阵methods 获知对指定类定义的所有方法函数mex 产生MEX文件编译环境的预设置指令mfunlis 能被mfun计算的MAPLE经典函数列表mhelp 引出Maple的在线帮助min 找向量中最小元素mkdir 创建目录mkpp 逐段多项式数据的明晰化mod 模运算more 指令窗中内容的分页显示movie 放映影片动画moviein 影片帧画面的内存预置mtaylor 符号计算多变量Taylor级数展开N nndims 求数组维数NaN 非数(预定义)变量nargchk 输入宗量数验证nargin 函数输入宗量数nargout 函数输出宗量数ndgrid 产生高维格点矩阵newplot 准备新的缺省图、轴nextpow2 取最接近的较大2次幂nnz 矩阵的非零元素总数nonzeros 矩阵的非零元素norm 矩阵或向量范数normcdf 正态分布累计概率密度函数normest 估计矩阵2范数norminv 正态分布逆累计概率密度函数normpdf 正态分布概率密度函数normrnd 正态随机数发生器notebook 启动Matlab和Word的集成环境null 零空间num2str 把非整数数组转换为串numden 获取最小公分母和相应的分子表达式nzmax 指定存放非零元素所需内存O oode1 非Stiff 微分方程变步长解算器ode15s Stiff 微分方程变步长解算器ode23t 适度Stiff 微分方程解算器ode23tb Stiff 微分方程解算器ode45 非Stiff 微分方程变步长解算器odefile ODE 文件模板odeget 获知ODE 选项设置参数odephas2 ODE 输出函数的二维相平面图odephas3 ODE 输出函数的三维相空间图odeplot ODE 输出函数的时间轨迹图odeprint 在Matlab指令窗显示结果odeset 创建或改写ODE选项构架参数值ones 全1数组optimset 创建或改写优化泛函指令的选项参数值orient 设定图形的排放方式orth 值空间正交化P ppack 收集Matlab内存碎块扩大内存pagedlg 调出图形排版对话框patch 创建块对象path 设置Matlab搜索路径的指令pathtool 搜索路径管理器pause 暂停pcode 创建预解译P码文件pcolor 伪彩图peaks Matlab提供的典型三维曲面permute 广义转置pi (预定义变量)圆周率pie 二维饼图pie3 三维饼图pink 粉红色图矩阵pinv 伪逆plot 平面线图plot3 三维线图plotmatrix 矩阵的散点图plotyy 双纵坐标图poissinv 泊松分布逆累计概率分布函数poissrnd 泊松分布随机数发生器pol2cart 极或柱坐标变为直角坐标polar 极坐标图poly 矩阵的特征多项式、根集对应的多项式poly2str 以习惯方式显示多项式poly2sym 双精度多项式系数转变为向量符号多项式polyder 多项式导数polyfit 数据的多项式拟合polyval 计算多项式的值polyvalm 计算矩阵多项式pow2 2的幂ppval 计算分段多项式pretty 以习惯方式显示符号表达式print 打印图形或SIMULINK模型printsys 以习惯方式显示有理分式prism 光谱色图矩阵procread 向MAPLE输送计算程序profile 函数文件性能评估器propedit 图形对象属性编辑器pwd 显示当前工作目录Q qquad 低阶法计算数值积分quad8 高阶法计算数值积分(QUADL) quit 推出Matlab 环境quiver 二维方向箭头图quiver3 三维方向箭头图R rrand 产生均匀分布随机数randn 产生正态分布随机数randperm 随机置换向量range 样本极差rank 矩阵的秩rats 有理输出rcond 矩阵倒条件数估计real 复数的实部reallog 在实数域内计算自然对数realpow 在实数域内计算乘方realsqrt 在实数域内计算平方根realmax 最大正浮点数realmin 最小正浮点数rectangle 画"长方框"rem 求余数repmat 铺放模块数组reshape 改变数组维数、大小residue 部分分式展开return 返回ribbon 把二维曲线画成三维彩带图rmfield 删去构架的域roots 求多项式的根rose 数扇形图rot90 矩阵旋转90度rotate 指定的原点和方向旋转rotate3d 启动三维图形视角的交互设置功能round 向最近整数圆整rref 简化矩阵为梯形形式rsf2csf 实数块对角阵转为复数特征值对角阵rsums Riemann和S ssave 把内存变量保存为文件scatter 散点图scatter3 三维散点图sec 正割sech 双曲正割semilogx X轴对数刻度坐标图semilogy Y轴对数刻度坐标图series 串联连接set 设置图形对象属性setfield 设置构架数组的域setstr 将ASCII码转换为字符的旧版指令sign 根据符号取值函数signum 符号计算中的符号取值函数sim 运行SIMULINK模型simget 获取SIMULINK模型设置的仿真参数simple 寻找最短形式的符号解simplify 符号计算中进行简化操作simset 对SIMULINK模型的仿真参数进行设置simulink 启动SIMULINK模块库浏览器sin 正弦sinh 双曲正弦size 矩阵的大小slice 立体切片图solve 求代数方程的符号解spalloc 为非零元素配置内存sparse 创建稀疏矩阵spconvert 把外部数据转换为稀疏矩阵spdiags 稀疏对角阵spfun 求非零元素的函数值sph2cart 球坐标变为直角坐标sphere 产生球面spinmap 色图彩色的周期变化spline 样条插值spones 用1置换非零元素sprandsym 稀疏随机对称阵sprank 结构秩spring 紫黄调春色图sprintf 把格式数据写成串spy 画稀疏结构图sqrt 平方根sqrtm 方根矩阵squeeze 删去大小为1的"孤维" sscanf 按指定格式读串stairs 阶梯图std 标准差stem 二维杆图step 阶跃响应指令str2double 串转换为双精度值str2mat 创建多行串数组str2num 串转换为数strcat 接成长串strcmp 串比较strjust 串对齐strmatch 搜索指定串strncmp 串中前若干字符比较strrep 串替换strtok 寻找第一间隔符前的内容struct 创建构架数组struct2cell 把构架转换为元胞数组strvcat 创建多行串数组sub2ind 多下标转换为单下标subexpr 通过子表达式重写符号对象subplot 创建子图subs 符号计算中的符号变量置换subspace 两子空间夹角sum 元素和summer 绿黄调夏色图superiorto 设定优先级surf 三维着色表面图surface 创建面对象surfc 带等位线的表面图surfl 带光照的三维表面图surfnorm 空间表面的法线svd 奇异值分解svds 求指定的若干奇异值switch-case-otherwise 多分支结构sym2poly 符号多项式转变为双精度多项式系数向量symmmd 对称最小度排序symrcm 反向Cuthill-McKee排序syms 创建多个符号对象T ttan 正切tanh 双曲正切taylortool 进行Taylor逼近分析的交互界面text 文字注释tf 创建传递函数对象tic 启动计时器title 图名toc 关闭计时器trapz 梯形法数值积分treelayout 展开树、林treeplot 画树图tril 下三角阵trim 求系统平衡点trimesh 不规则格点网线图trisurf 不规则格点表面图triu 上三角阵try-catch 控制流中的Try-catch结构type 显示M 文件U uuicontextmenu 创建现场菜单uicontrol 创建用户控件uimenu 创建用户菜单unmkpp 逐段多项式数据的反明晰化unwrap 自然态相角upper 转换为大写字母V vvar 方差varargin 变长度输入宗量varargout 变长度输出宗量vectorize 使串表达式或内联函数适于数组运算ver 版本信息的获取view 三维图形的视角控制voronoi Voronoi多边形vpa 任意精度(符号类)数值W wwarning 显示警告信息what 列出当前目录上的文件whatsnew 显示Matlab中Readme文件的内容which 确定函数、文件的位置while 控制流中的While环结构white 全白色图矩阵whitebg 指定轴的背景色who 列出内存中的变量名whos 列出内存中变量的详细信息winter 蓝绿调冬色图workspace 启动内存浏览器X x , Y y , Z zxlabel X轴名xor 或非逻辑yesinput 智能输入指令ylabel Y轴名zeros 全零数组zlabel Z轴名zoom 图形的变焦放大和缩小ztrans 符号计算Z变换v。

基于串行双环微环谐振器的滤波特性

基于串行双环微环谐振器的滤波特性

基于串行双环微环谐振器的滤波特性伴随科技的发展,人们对光的认知深度和广度不断增加。

光学技术的不断发展,对社会以及科学的发展有巨大的促进作用,因此光信号的传输和处理成为了光学运用的关键。

我们必须找到一种能够控制光产生,运输以及处理的工具,其中微环谐振器就是一个很好的控制光传输和处理的工具。

微环谐振器具有结构紧凑,高度集成,功能丰富的特点,是搭建大规模及超大规模集成光子回路的最具有前途的代表之一。

目前,它已被用于延时线,通信滤波器,微激光器和光学传感器等,成为光通信和集成光电子领域的研究热点之一。

硅电子器件制备工艺由于技术成熟,批量生产成本低等优势,硅基微环谐振器刚好充分利用这一技术。

本文中,利用建立的微环谐振器的数学模型即散射矩阵模型,着重研究基于硅微环谐振器的光子回路的滤波特性,包括自由光谱范围,3dB带宽等参数。

同时利用了FDTD算法和MATLAB进行适当的仿真,对理论分析的结果加以验证。

关键字:微环,谐振器,滤波,滤波特性,散射矩阵绪言很多年前,Marcatili 在自己的文章中清晰地介绍了微环波导谐振腔的工作原理。

一系列频率的光束由波导1的输入端口输入,在波导中传播,若输入光波的频率符合微环谐振器的谐振频率,该频率的光就可以在微环谐振腔内与新耦合的光波发生干涉。

该频率的光波在波导2的下载端口输出,其他频率的信号则在波导1的输出端口输出。

虽然提出微环谐振器的概念已经很长时间,但在一段时间内,由于平面工艺水平没有大的突破,所以微环谐振器的研究和发展一直停滞不前。

后来随着平面工艺技术的提升,微环谐振器才重新得到科研技术工作者的重视和研究,并得以迅速发展。

作为搭建光子回路的一个重要组成部分,微环谐振器已被运用于许多方面,比如半导体激光器,因为微环谐振器可以承担起很多功能,其已成为光集成的一个极为重要的组成部分。

随着集成光学研究领域的不断发展,特别是近几年微纳光集成工艺取得的突破性进展,环形谐振器逐渐引起全球各科研团体的兴趣和重视。

基于matlab程序的16qam基带信号处理编程.

基于matlab程序的16qam基带信号处理编程.

摘要:正交振幅调制技术(QAM) 由于其高的频带利用率和相对低的误码率而被定为很多数字通信系统的数字传输标准。

该论文讨论了16QAM 调制原理,用软件MATLAB 产生16QAM 信号的算法。

关键词:正交振幅调制;噪声; 多径衰落;内插滤波一、 引言正交振幅调制技术(QAM)是幅度、相位联合调制的技术,它同时利用了载波的幅度和相位来传递信息比特,因此在相邻星座点最小距离相同的条件下,QAM 星座图中可以容纳更多的星座点,即可实现更高的频带利用率。

多进制正交振幅调制技术由于其高的频带利用率和相对低的误码率而被定为很多数字通信系统的数字传输标准。

软件无线电是以可编程的DSP 或CPU 为中心,将模块化、标准化的硬件单元以总线方式连接件起来,构成通用的基本硬件平台,并通过软件加载来实现各种无线通信功能的开放式的体系结构。

本文讨论了16QAM 调制原理并给出了基于Matlab 的16QAM 数字调制的软件算法及其编程实现。

二、 正交振幅调制(QAM )原理正交振幅调制的一般表达式为:2()Re[()()]()cos(2)()sin(2)c j f t MQAM mc ms mc c ms c S t A jA g t e A g t f t A g t f t πππ=+=-上式中,mc A 、ms A 为离散的振幅值,m=1,2…,M ,M 为mc A 和ms A 的个数。

对于 16QAM ,M=4。

()g t 是发送信号脉冲波形。

本报告中16QAM 的调制原理方框图如下图所示:输入数据经过串并变换分成两路,再分别经过映射变换成mc A 和ms A ,mc A 和ms A 决定了已调QAM 信号在信号星座中的星座点的位置。

为减少码间干扰和已调数字信号的带宽,mc A 和ms A 分别通过脉冲成形低通滤波器,脉冲成形低通滤波器是一个平方根升余弦滤波器。

滤波后的输出分别经过半带滤波器、CIC 滤波器实现内插,改变速率,最后经过NCO(数控频率振荡器)产生正交上变频用的正/余弦信号。

Matlab系统仿真代码技术文档

Matlab系统仿真代码技术文档

流程示意图(a)发射端框图CSI流程示意图 (1)1 主程序ofdm_wav_trans_system.m (2)2 D/A转换模块:Get_wav_data.m (9)3 长短帧设定:construct_flag.m (10)4 子载波参数设定:bitload.m (11)5 数据读取Fetch_data_old.m (11)6 数据读取:Fetch_data.m (12)7 加扰模块Scramble_data_old.m (13)8 加扰模块Scramble_data (14)9 卷积码encd216.m (14)10 交织模块1 index_change.m (16)11 交织模块2 index_change2.m (18)12 交织模块3 Data_interleave.m (20)13 导频模块:add_prefix.m (21)14 QAM模块:modulate.m (22)15 成型滤波器模块:shape_filtering.m (24)16 信道模拟模块:channel.m (25)17 匹配滤波器:match_filtering.m (26)18 欠抽样恢复模块:sample_decision.m (28)19 同步模块:fresym_combine_est.m (30)20 解帧模块:decompose_frame.m (32)21 频偏估计:get_sample_clock_offset.m (33)22 频偏补偿:recov_sample_clock.m (34)23 QAM解映射:demodulate.m (35)24 解交织模块1:index_normal2.m (37)25 解交织模块2:index_normal.m (38)26 解交织模块3:Data_Deinterleave.m (39)27 V iterbi译码模块:decd216.m (40)1 主程序ofdm_wav_trans_system.m功能:描述数据在数字域的变化,将音频文件转为数字,组帧,加扰,解扰clear all;tic;%帧参数global GUARD_num; GUARD_num=16; %循环前缀所占global FFT_num; FFT_num=64; %每符号子载波数global OFDM_symbols_per_frame; OFDM_symbols_per_frame =25; %每帧数据符号数global OFDM_prefix_per_frame; OFDM_prefix_per_frame=5; %每帧前导符号数global OFDM_wavbits_per_frame1;OFDM_wavbits_per_frame=1408; %每帧比特数global OFDM_wavbits_per_frame2; OFDM_wavbits_per_frame2=1440;global OFDM_scrambits_per_frame;OFDM_scrambits_per_frame=1440; %每帧扰码后比特数global Total_bits_per_symbol; Total_bits_per_symbol=144; %每个符号的比特数36×4scram_polynomial=[1 0 0 0 1 0 0 1]; %扰码多项式SNR=10; %单位:DBrolloff_factor=0.25;up_sampling_factor=4;%帧结构% 长训练符号LTS = [0 0 0 0 0 0 1 1 -1 -1 1 1 -1 1 -1 1 ...1 1 1 1 1 -1 -1 1 1 -1 1 -1 1 1 1 1 ...0 1 -1 -1 1 1 -1 1 -1 1 -1 -1 -1 -1 -1 1 ...1 -1 -1 1 -1 1 -1 1 1 1 1 0 0 0 0 0];% 短训练符号STS = sqrt(13/6)* ...[0 0 0 0 0 0 0 0 1+j 0 0 0 -1-j 0 0 0 ...1+j 0 0 0 -1-j 0 0 0 -1-j 0 0 0 1+j 0 0 0 ...0 0 0 0 -1-j 0 0 0 -1-j 0 0 0 1+j 0 0 0 ...1+j 0 0 0 1+j 0 0 0 1+j 0 0 0 0 0 0 0]; ...global SIGNAL1;% 信令符号SIGNAL1 = zeros(1,144);global SIGNAL2;SIGNAL2 = ones(1,144);% 导频PILOT=[-1 1 1 1];Pattern_Datas_Symbols_Index = [8:15 17:26 40:49 51:58];Pattern_Pilot_Symbols_Index = [16 27 39 50];% 功率放大器参数HPA_para_AmAm_alpha = 2.1587;HPA_para_AmAm_beta = 1.1517;HPA_para_AmPm_alpha = 4.0033;HPA_para_AmPm_beta = 9.1040;HPA_para_InputScaling = 1; % Not dbHPA_para_OutputScaling = 1;% Not dbHPA_para=[HPA_para_AmAm_alpha,HPA_para_AmAm_beta;HPA_para_AmPm_alpha,HPA_para_AmPm_beta; HPA_para_InputScaling,HPA_para_OutputScaling];chan_multipath_imp_res_upsampling=[1 0.89 0.81 0.74 0.68 0.62 0.57 0.52 0.48 0.46 0.44 0.43 0.42 ...0.41 0.405 0.402 0.40 0.392 0.387 0.383 0.38 0.36 0.34 0.33 0.32 ...0.27 0.23 0.19 0.17 0.12 0.09 0.06]; %多径效应trans_carrier_num=length(Pattern_Datas_Symbols_Index); %数字载波的子载波数energy_per_carrier=1; %符号能量load ENC2.matload ENC4.mat %qpsk对照表load ENC16.mat %16QAM对照表load ENC64.mat %64QAM对照表load ENC256.mat %256QAM对照表% I/Q ImbalenceIQ_Imb_am = 0.025;IQ_Imb_ph = 1*pi*2/360;[Wav_data,sample_rate,per_sample_bits] =Get_wav_data('tada.wav'); %将声音文件转化为数据trans_frame_wav_data=[];data_index=1;%flag_in=[0 0 0 0 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1];flag_in=[0 0 0 0 1];flag=flag_in;%flag=construct_flag(flag_in); % 产生flag序列% trans_iter_num=ceil(length(Trans_wav_data)/OFDM_wavbits_per_frame); %所需传输帧数rec_wav_data=[];[trans_bits_alloctrans_energy_alloc]=bitload(Pattern_Datas_Symbols_Index,Total_bits_per_symbol,tran s_carrier_num,energy_per_carrier);%chan_multipath_fre=fft(chan_multipath_imp_res,FFT_num);%频域信道响应for Trans_frame_index =1:1:length(flag)flag_value=flag(Trans_frame_index);[Frame_wav_data,get_signal_data,new_begin_index] =Fetch_data(Wav_data,data_index,flag_value); % 基于缓冲池的状态决定帧长% 并根据flag决定信令data_index=new_begin_index;trans_frame_wav_data=[trans_frame_wav_data Frame_wav_data];T_data=[get_signal_data Frame_wav_data];T_data=Scramble_data(T_data);% 加扰,信令144比特+1440比特数据Scram_signal_data=T_data(1:144);% 前面144比特为加扰后的信令Scram_Frame_wav_data=T_data(145:1584);% 后面是1440比特加扰后的数据coded_Frame_wav_data=encd216(Scram_Frame_wav_data);Frame_scram_data=[Scram_signal_data coded_Frame_wav_data]; %两次交织总效果%temp_data=index_change( Frame_scram_data);% 第一次交织%trans_wav_data=index_change2(temp_data);% 第二次交织% Construct OFDM frameOFDM_trans_frame=zeros(25,80);OFDM_trans_frame=add_prefix(OFDM_trans_frame,STS,LTS); %生成帧头% 包括信令在内的21个符号,每个符号144比特,经过16QAM调制到相应的信道,并插入导频,这里信令也插入了导频data_begin_index=1;for k=1:OFDM_symbols_per_frame-OFDM_prefix_per_frame+1trans_bits_per_symbol=trans_wav_data(data_begin_index:data_begin_index+Total_bits_per_symbol-1);trans_16QAM_per_symbol = zeros(1,FFT_num);trans_16QAM_per_symbol=modulate(trans_bits_per_symbol,trans_bits_alloc,trans_energy_alloc,s2,s4,s16,s64,s256);trans_16QAM_per_symbol(Pattern_Pilot_Symbols_Index)=PILOT;%加入导频符号trans_tmdata_per_symbol=ifft(trans_16QAM_per_symbol,FFT_num);%进行ifft变换,相当于系统模块中OFDM调制OFDM_trans_frame(k+4,:)=[trans_tmdata_per_symbol(49:64)trans_tmdata_per_symbol]; % 数据帧加循环前缀data_begin_index=data_begin_index+Total_bits_per_symbol;endOFDM_trans_frame=transpose(OFDM_trans_frame); %矩阵转置,变成时间域表示形式% 以下步骤形成一个数据帧tx_frame_signal =reshape(OFDM_trans_frame,1,OFDM_symbols_per_frame*(GUARD_num+FFT_num));% 形成一个2000个采样点的数据帧,参考批注4。

硅基微环谐振器中耦合损耗的计算方法[发明专利]

硅基微环谐振器中耦合损耗的计算方法[发明专利]

(10)申请公布号 CN 102692784 A (43)申请公布日 2012.09.26C N 102692784 A*CN102692784A*(21)申请号 201210204101.1(22)申请日 2012.06.20G02F 3/00(2006.01)G02F 1/365(2006.01)(71)申请人哈尔滨工业大学地址150000 黑龙江省哈尔滨市南岗区西大直街92号(72)发明人左宝君 胡海力 范志刚 陈守谦刘建军(54)发明名称硅基微环谐振器中耦合损耗的计算方法(57)摘要硅基微环谐振器中耦合损耗的计算方法,涉及一种基于逐步近似法的TE 模式的硅基微环谐振器耦合损耗的计算方法。

鉴于目前没有针对硅基微环谐振器的耦合损耗计算方法的问题,本发明的硅基微环谐振器中耦合损耗的计算方法,首先利用逐步近似方法,将微环谐振器的环形波导近似为多个相互平行的直线型波导组;然后,利用波导方程和TE 模式的边界条件,对波导组分别计算光波在其中的模式转换系数和透射系数,计算在每个波导组由于模式转换效应造成的能量损耗;最终,整个耦合区域的耦合损耗就是每个波导组能量损耗的总和。

本发明的耦合损耗计算方法,填补了计算耦合损耗的理论方法的空白。

本发明方法用MATLAB 即可实现,计算时间快。

(51)Int.Cl.权利要求书2页 说明书6页 附图3页(19)中华人民共和国国家知识产权局(12)发明专利申请权利要求书 2 页 说明书 6 页 附图 3 页1.硅基微环谐振器中耦合损耗的计算方法,其特征在于所述计算方法如下:首先,利用逐步近似方法,将微环谐振器的环形波导近似为多个相互平行的直线型波导组;然后,利用波导方程和TE模式的边界条件,对波导组分别计算光波在其中的模式转换系数和透射系数,计算在每个波导组由于模式转换效应造成的能量损耗;最终,整个耦合区域的耦合损耗就是每个波导组能量损耗的总和。

2.根据权利要求1所述的硅基微环谐振器中耦合损耗的计算方法,其特征在于所述微环谐振器为硅基超紧凑型微环谐振器,光波导的横截面尺寸最小为220nm×445nm,光波长为1.55~1.70μm。

matlab 简谐振动 程序

matlab 简谐振动 程序

matlab 简谐振动程序
简谐振动是物理学中重要的概念,可以用来描述许多物理系统的运动,比如弹簧振子和摆动等。

在MATLAB中,我们可以编写程序来模拟简谐振动的运动。

以下是一个简单的MATLAB程序,用来模拟简谐振动的运动:
matlab.
% 定义参数。

m = 1; % 质量。

k = 1; % 弹簧常数。

x0 = 1; % 初始位移。

v0 = 0; % 初始速度。

omega = sqrt(k/m); % 角频率。

% 定义时间变量。

t = 0:0.1:10; % 时间范围。

% 计算位移。

x = x0 cos(omega t) + (v0/omega) sin(omega t);
% 绘制位移-时间图。

plot(t, x);
xlabel('时间');
ylabel('位移');
title('简谐振动的位移-时间图');
在这个程序中,我们首先定义了简谐振动的参数,包括质量m、弹簧常数k、初始位移x0和初始速度v0。

然后我们计算了在给定时
间范围内的位移x,并绘制了位移-时间图。

当你运行这个程序时,你会得到一个简谐振动的位移-时间图,它展示了振动的周期性运动。

你可以根据自己的需要修改参数和时间范围,来观察不同条件下简谐振动的运动情况。

除了这个简单的例子,你还可以编写更复杂的程序来模拟不同类型的简谐振动,比如受阻尼力和外力作用下的简谐振动。

希望这个简单的程序可以帮助你开始理解如何在MATLAB中模拟简谐振动。

matlab命令大全 (全面)

matlab命令大全 (全面)

分类:VC MatlabA aabs 绝对值、模、字符的ASCII码值acos 反余弦acosh 反双曲余弦acot 反余切acoth 反双曲余切acsc 反余割acsch 反双曲余割align 启动图形对象几何位置排列工具all 所有元素非零为真angle 相角ans 表达式计算结果的缺省变量名any 所有元素非全零为真area 面域图argnames 函数M文件宗量名asec 反正割asech 反双曲正割asin 反正弦asinh 反双曲正弦assignin 向变量赋值atan 反正切atan2 四象限反正切atanh 反双曲正切autumn 红黄调秋色图阵axes 创建轴对象的低层指令axis 控制轴刻度和风格的高层指令B bbar 二维直方图bar3 三维直方图bar3h 三维水平直方图barh 二维水平直方图base2dec X进制转换为十进制bin2dec 二进制转换为十进制blanks 创建空格串bone 蓝色调黑白色图阵box 框状坐标轴break while 或for 环中断指令brighten 亮度控制C ccapture (3版以前)捕获当前图形cart2pol 直角坐标变为极或柱坐标cart2sph 直角坐标变为球坐标cat 串接成高维数组caxis 色标尺刻度cd 指定当前目录cdedit 启动用户菜单、控件回调函数设计工具cdf2rdf 复数特征值对角阵转为实数块对角阵ceil 向正无穷取整cell 创建元胞数组cell2struct 元胞数组转换为构架数组celldisp 显示元胞数组内容cellplot 元胞数组内部结构图示char 把数值、符号、内联类转换为字符对象chi2cdf 分布累计概率函数chi2inv 分布逆累计概率函数chi2pdf 分布概率密度函数chi2rnd 分布随机数发生器chol Cholesky分解clabel 等位线标识cla 清除当前轴class 获知对象类别或创建对象clc 清除指令窗clear 清除内存变量和函数clf 清除图对象clock 时钟colorcube 三浓淡多彩交叉色图矩阵colordef 设置色彩缺省值colormap 色图colspace 列空间的基close 关闭指定窗口colperm 列排序置换向量comet 彗星状轨迹图comet3 三维彗星轨迹图compass 射线图compose 求复合函数cond (逆)条件数condeig 计算特征值、特征向量同时给出条件数condest 范-1条件数估计conj 复数共轭contour 等位线contourf 填色等位线contour3 三维等位线contourslice 四维切片等位线图conv 多项式乘、卷积cool 青紫调冷色图copper 古铜调色图cos 余弦cosh 双曲余弦cot 余切coth 双曲余切cplxpair 复数共轭成对排列csc 余割csch 双曲余割cumsum 元素累计和cumtrapz 累计梯形积分cylinder 创建圆柱D ddblquad 二重数值积分deal 分配宗量deblank 删去串尾部的空格符dec2base 十进制转换为X进制dec2bin 十进制转换为二进制dec2hex 十进制转换为十六进制deconv 多项式除、解卷delaunay Delaunay 三角剖分del2 离散Laplacian差分demo Matlab演示det 行列式diag 矩阵对角元素提取、创建对角阵diary Matlab指令窗文本内容记录diff 数值差分、符号微分digits 符号计算中设置符号数值的精度dir 目录列表disp 显示数组display 显示对象内容的重载函数dlinmod 离散系统的线性化模型dmperm 矩阵Dulmage-Mendelsohn 分解dos 执行DOS 指令并返回结果double 把其他类型对象转换为双精度数值drawnow 更新事件队列强迫Matlab刷新屏幕dsolve 符号计算解微分方程E eecho M文件被执行指令的显示edit 启动M文件编辑器eig 求特征值和特征向量eigs 求指定的几个特征值end 控制流FOR等结构体的结尾元素下标eps 浮点相对精度error 显示出错信息并中断执行errortrap 错误发生后程序是否继续执行的控制erf 误差函数erfc 误差补函数erfcx 刻度误差补函数erfinv 逆误差函数errorbar 带误差限的曲线图etreeplot 画消去树eval 串演算指令evalin 跨空间串演算指令exist 检查变量或函数是否已定义exit 退出Matlab环境exp 指数函数expand 符号计算中的展开操作expint 指数积分函数expm 常用矩阵指数函数expm1 Pade法求矩阵指数expm2 Taylor法求矩阵指数expm3 特征值分解法求矩阵指数eye 单位阵ezcontour 画等位线的简捷指令ezcontourf 画填色等位线的简捷指令ezgraph3 画表面图的通用简捷指令ezmesh 画网线图的简捷指令ezmeshc 画带等位线的网线图的简捷指令ezplot 画二维曲线的简捷指令ezplot3 画三维曲线的简捷指令ezpolar 画极坐标图的简捷指令ezsurf 画表面图的简捷指令ezsurfc 画带等位线的表面图的简捷指令F ffactor 符号计算的因式分解feather 羽毛图feedback 反馈连接feval 执行由串指定的函数fft 离散Fourier变换fft2 二维离散Fourier变换fftn 高维离散Fourier变换fftshift 直流分量对中的谱fieldnames 构架域名figure 创建图形窗fill3 三维多边形填色图find 寻找非零元素下标findobj 寻找具有指定属性的对象图柄findstr 寻找短串的起始字符下标findsym 机器确定内存中的符号变量finverse 符号计算中求反函数fix 向零取整flag 红白蓝黑交错色图阵fliplr 矩阵的左右翻转flipud 矩阵的上下翻转flipdim 矩阵沿指定维翻转floor 向负无穷取整flops 浮点运算次数flow Matlab提供的演示数据fmin 求单变量非线性函数极小值点(旧版)fminbnd 求单变量非线性函数极小值点fmins 单纯形法求多变量函数极小值点(旧版)fminunc 拟牛顿法求多变量函数极小值点fminsearch 单纯形法求多变量函数极小值点fnder 对样条函数求导fnint 利用样条函数求积分fnval 计算样条函数区间内任意一点的值fnplt 绘制样条函数图形fopen 打开外部文件for 构成for环用format 设置输出格式fourier Fourier 变换fplot 返函绘图指令fprintf 设置显示格式fread 从文件读二进制数据fsolve 求多元函数的零点full 把稀疏矩阵转换为非稀疏阵funm 计算一般矩阵函数funtool 函数计算器图形用户界面fzero 求单变量非线性函数的零点G ggamma 函数gammainc 不完全函数gammaln 函数的对数gca 获得当前轴句柄gcbo 获得正执行"回调"的对象句柄gcf 获得当前图对象句柄gco 获得当前对象句柄geomean 几何平均值get 获知对象属性getfield 获知构架数组的域getframe 获取影片的帧画面ginput 从图形窗获取数据global 定义全局变量gplot 依图论法则画图gradient 近似梯度gray 黑白灰度grid 画分格线griddata 规则化数据和曲面拟合gtext 由鼠标放置注释文字guide 启动图形用户界面交互设计工具H hharmmean 调和平均值help 在线帮助helpwin 交互式在线帮助helpdesk 打开超文本形式用户指南hex2dec 十六进制转换为十进制hex2num 十六进制转换为浮点数hidden 透视和消隐开关hilb Hilbert矩阵hist 频数计算或频数直方图histc 端点定位频数直方图histfit 带正态拟合的频数直方图hold 当前图上重画的切换开关horner 分解成嵌套形式hot 黑红黄白色图hsv 饱和色图I iif-else-elseif 条件分支结构ifft 离散Fourier反变换ifft2 二维离散Fourier反变换ifftn 高维离散Fourier反变换ifftshift 直流分量对中的谱的反操作ifourier Fourier反变换i, j 缺省的"虚单元"变量ilaplace Laplace反变换imag 复数虚部image 显示图象imagesc 显示亮度图象imfinfo 获取图形文件信息imread 从文件读取图象imwrite 把imwrite 把图象写成文件ind2sub 单下标转变为多下标inf 无穷大info MathWorks公司网点地址inline 构造内联函数对象inmem 列出内存中的函数名input 提示用户输入inputname 输入宗量名int 符号积分int2str 把整数数组转换为串数组interp1 一维插值interp2 二维插值interp3 三维插值interpn N维插值interpft 利用FFT插值intro Matlab自带的入门引导inv 求矩阵逆invhilb Hilbert矩阵的准确逆ipermute 广义反转置isa 检测是否给定类的对象ischar 若是字符串则为真isequal 若两数组相同则为真isempty 若是空阵则为真isfinite 若全部元素都有限则为真isfield 若是构架域则为真isglobal 若是全局变量则为真ishandle 若是图形句柄则为真ishold 若当前图形处于保留状态则为真isieee 若计算机执行IEEE规则则为真isinf 若是无穷数据则为真isletter 若是英文字母则为真islogical 若是逻辑数组则为真ismember 检查是否属于指定集isnan 若是非数则为真isnumeric 若是数值数组则为真isobject 若是对象则为真isprime 若是质数则为真isreal 若是实数则为真isspace 若是空格则为真issparse 若是稀疏矩阵则为真isstruct 若是构架则为真isstudent 若是Matlab学生版则为真iztrans 符号计算Z反变换J j , K kjacobian 符号计算中求Jacobian 矩阵jet 蓝头红尾饱和色jordan 符号计算中获得Jordan标准型keyboard 键盘获得控制权kron Kronecker乘法规则产生的数组L llaplace Laplace变换lasterr 显示最新出错信息lastwarn 显示最新警告信息leastsq 解非线性最小二乘问题(旧版)legend 图形图例lighting 照明模式line 创建线对象lines 采用plot 画线色linmod 获连续系统的线性化模型linmod2 获连续系统的线性化精良模型linspace 线性等分向量ln 矩阵自然对数load 从MAT文件读取变量log 自然对数log10 常用对数log2 底为2的对数loglog 双对数刻度图形logm 矩阵对数logspace 对数分度向量lookfor 按关键字搜索M文件lower 转换为小写字母lsqnonlin 解非线性最小二乘问题lu LU分解M mmad 平均绝对值偏差magic 魔方阵maple &nb, sp; 运作Maple格式指令mat2str 把数值数组转换成输入形态串数组material 材料反射模式max 找向量中最大元素mbuild 产生EXE文件编译环境的预设置指令mcc 创建MEX或EXE文件的编译指令mean 求向量元素的平均值median 求中位数menuedit 启动设计用户菜单的交互式编辑工具mesh 网线图meshz 垂帘网线图meshgrid 产生"格点"矩阵methods 获知对指定类定义的所有方法函数mex 产生MEX文件编译环境的预设置指令mfunlis 能被mfun计算的MAPLE经典函数列表mhelp 引出Maple的在线帮助min 找向量中最小元素mkdir 创建目录mkpp 逐段多项式数据的明晰化mod 模运算more 指令窗中内容的分页显示movie 放映影片动画moviein 影片帧画面的内存预置mtaylor 符号计算多变量Taylor级数展开N nndims 求数组维数NaN 非数(预定义)变量nargchk 输入宗量数验证nargin 函数输入宗量数nargout 函数输出宗量数ndgrid 产生高维格点矩阵newplot 准备新的缺省图、轴nextpow2 取最接近的较大2次幂nnz 矩阵的非零元素总数nonzeros 矩阵的非零元素norm 矩阵或向量范数normcdf 正态分布累计概率密度函数normest 估计矩阵2范数norminv 正态分布逆累计概率密度函数normpdf 正态分布概率密度函数normrnd 正态随机数发生器notebook 启动Matlab和Word的集成环境null 零空间num2str 把非整数数组转换为串numden 获取最小公分母和相应的分子表达式nzmax 指定存放非零元素所需内存O oode1 非Stiff 微分方程变步长解算器ode15s Stiff 微分方程变步长解算器ode23t 适度Stiff 微分方程解算器ode23tb Stiff 微分方程解算器ode45 非Stiff 微分方程变步长解算器odefile ODE 文件模板odeget 获知ODE 选项设置参数odephas2 ODE 输出函数的二维相平面图odephas3 ODE 输出函数的三维相空间图odeplot ODE 输出函数的时间轨迹图odeprint 在Matlab指令窗显示结果odeset 创建或改写ODE选项构架参数值ones 全1数组optimset 创建或改写优化泛函指令的选项参数值orient 设定图形的排放方式orth 值空间正交化P ppack 收集Matlab内存碎块扩大内存pagedlg 调出图形排版对话框patch 创建块对象path 设置Matlab搜索路径的指令pathtool 搜索路径管理器pause 暂停pcode 创建预解译P码文件pcolor 伪彩图peaks Matlab提供的典型三维曲面permute 广义转置pi (预定义变量)圆周率pie 二维饼图pie3 三维饼图pink 粉红色图矩阵pinv 伪逆plot 平面线图plot3 三维线图plotmatrix 矩阵的散点图plotyy 双纵坐标图poissinv 泊松分布逆累计概率分布函数poissrnd 泊松分布随机数发生器pol2cart 极或柱坐标变为直角坐标polar 极坐标图poly 矩阵的特征多项式、根集对应的多项式poly2str 以习惯方式显示多项式poly2sym 双精度多项式系数转变为向量符号多项式polyder 多项式导数polyfit 数据的多项式拟合polyval 计算多项式的值polyvalm 计算矩阵多项式pow2 2的幂ppval 计算分段多项式pretty 以习惯方式显示符号表达式print 打印图形或SIMULINK模型printsys 以习惯方式显示有理分式prism 光谱色图矩阵procread 向MAPLE输送计算程序profile 函数文件性能评估器propedit 图形对象属性编辑器pwd 显示当前工作目录Q qquad 低阶法计算数值积分quad8 高阶法计算数值积分(QUADL)quit 推出Matlab 环境quiver 二维方向箭头图quiver3 三维方向箭头图R rrand 产生均匀分布随机数randn 产生正态分布随机数randperm 随机置换向量range 样本极差rank 矩阵的秩rats 有理输出rcond 矩阵倒条件数估计real 复数的实部reallog 在实数域内计算自然对数realpow 在实数域内计算乘方realsqrt 在实数域内计算平方根realmax 最大正浮点数realmin 最小正浮点数rectangle 画"长方框"rem 求余数repmat 铺放模块数组reshape 改变数组维数、大小residue 部分分式展开return 返回ribbon 把二维曲线画成三维彩带图rmfield 删去构架的域roots 求多项式的根rose 数扇形图rot90 矩阵旋转90度rotate 指定的原点和方向旋转rotate3d 启动三维图形视角的交互设置功能round 向最近整数圆整rref 简化矩阵为梯形形式rsf2csf 实数块对角阵转为复数特征值对角阵rsums Riemann和S ssave 把内存变量保存为文件scatter 散点图scatter3 三维散点图sec 正割sech 双曲正割semilogx X轴对数刻度坐标图semilogy Y轴对数刻度坐标图series 串联连接set 设置图形对象属性setfield 设置构架数组的域setstr 将ASCII码转换为字符的旧版指令sign 根据符号取值函数signum 符号计算中的符号取值函数sim 运行SIMULINK模型simget 获取SIMULINK模型设置的仿真参数simple 寻找最短形式的符号解simplify 符号计算中进行简化操作simset 对SIMULINK模型的仿真参数进行设置simulink 启动SIMULINK模块库浏览器sin 正弦sinh 双曲正弦size 矩阵的大小slice 立体切片图solve 求代数方程的符号解spalloc 为非零元素配置内存sparse 创建稀疏矩阵spconvert 把外部数据转换为稀疏矩阵spdiags 稀疏对角阵spfun 求非零元素的函数值sph2cart 球坐标变为直角坐标sphere 产生球面spinmap 色图彩色的周期变化spline 样条插值spones 用1置换非零元素sprandsym 稀疏随机对称阵sprank 结构秩spring 紫黄调春色图sprintf 把格式数据写成串spy 画稀疏结构图sqrt 平方根sqrtm 方根矩阵squeeze 删去大小为1的"孤维"sscanf 按指定格式读串stairs 阶梯图std 标准差stem 二维杆图step 阶跃响应指令str2double 串转换为双精度值str2mat 创建多行串数组str2num 串转换为数strcat 接成长串strcmp 串比较strjust 串对齐strmatch 搜索指定串strncmp 串中前若干字符比较strrep 串替换strtok 寻找第一间隔符前的内容struct 创建构架数组struct2cell 把构架转换为元胞数组strvcat 创建多行串数组sub2ind 多下标转换为单下标subexpr 通过子表达式重写符号对象subplot 创建子图subs 符号计算中的符号变量置换subspace 两子空间夹角sum 元素和summer 绿黄调夏色图superiorto 设定优先级surf 三维着色表面图surface 创建面对象surfc 带等位线的表面图surfl 带光照的三维表面图surfnorm 空间表面的法线svd 奇异值分解svds 求指定的若干奇异值switch-case-otherwise 多分支结构sym2poly 符号多项式转变为双精度多项式系数向量symmmd 对称最小度排序symrcm 反向Cuthill-McKee排序syms 创建多个符号对象T ttan 正切tanh 双曲正切taylortool 进行Taylor逼近分析的交互界面text 文字注释tf 创建传递函数对象tic 启动计时器title 图名toc 关闭计时器trapz 梯形法数值积分treelayout 展开树、林treeplot 画树图tril 下三角阵trim 求系统平衡点trimesh 不规则格点网线图trisurf 不规则格点表面图triu 上三角阵try-catch 控制流中的Try-catch结构type 显示M 文件U uuicontextmenu 创建现场菜单uicontrol 创建用户控件uimenu 创建用户菜单unmkpp 逐段多项式数据的反明晰化unwrap 自然态相角upper 转换为大写字母V vvar 方差varargin 变长度输入宗量varargout 变长度输出宗量vectorize 使串表达式或内联函数适于数组运算ver 版本信息的获取view 三维图形的视角控制voronoi Voronoi多边形vpa 任意精度(符号类)数值W wwarning 显示警告信息what 列出当前目录上的文件whatsnew 显示Matlab中Readme文件的内容which 确定函数、文件的位置while 控制流中的While环结构white 全白色图矩阵whitebg 指定轴的背景色who 列出内存中的变量名whos 列出内存中变量的详细信息winter 蓝绿调冬色图workspace 启动内存浏览器X x , Y y , Z zxlabel X轴名xor 或非逻辑yesinput 智能输入指令ylabel Y轴名zeros 全零数组zlabel Z轴名zoom 图形的变焦放大和缩小ztrans 符号计算Z变换MATLAB命令大全Prod 计算数组元素的连乘积。

Matlab Simulink建模与仿真例题源代码

Matlab Simulink建模与仿真例题源代码

〔實例3.1〕試設計一個模擬低通濾波器,f p = 2400Hz,f s = 5000Hz,R p = 3 dB,R s =25dB。

分別用巴特沃斯和橢圓濾波器原型,求出其3dB 截止頻率和濾波器階數,傳遞函數,並作出幅頻、相頻特性曲線。

巴特沃斯濾波器設計的程序代碼如下:% ch3example1A.mclear;f_p=2400; f_s=5000; R_p=3; R_s=25; % 設計要求指標[n, fn]=buttord(f_p,f_s,R_p,R_s, 's'); % 計算階數和截止頻率Wn=2*pi*fn; % 轉換為角頻率[b,a]=butter(n, Wn, 's'); % 計算H(s)f=0:100:10000; % 計算頻率點和頻率範圍s=j*2*pi*f; % s=jw=j*2*pi*fH_s=polyval(b,s)./polyval(a,s); % 計算相應頻率點處H(s)的值figure(1);subplot(2,1,1); plot(f, 20*log10(abs(H_s))); % 幅頻特性axis([0 10000 -40 1]);xlabel('頻率Hz');ylabel('幅度dB');subplot(2,1,2); plot(f, angle(H_s)); % 相頻特性xlabel('頻率Hz');ylabel('相角rad');figure(2); freqs(b,a); % 也可用指令freqs直接畫出H(s)的頻率響應曲線。

橢圓濾波器設計的程序代碼如下:〔程序代碼〕ch3example1B.m% ch3example1B.mclear;f_p=2400; f_s=5000; R_p=3; R_s=25; % 設計要求指標[n, fn]=ellipord(f_p,f_s,R_p,R_s,'s'); % 計算階數和截止頻率Wn=2*pi*fn; % 轉換為角頻率[b,a]=ellip(n,R_p,R_s,Wn,'s'); % 計算H(s)f=0:100:10000; % 計算頻率點和頻率範圍s=j*2*pi*f; % s=jw=j*2*pi*fH_s=polyval(b,s)./polyval(a,s); % 計算相應頻率點處H(s)的值figure(1);subplot(2,1,1); plot(f, 20*log10(abs(H_s))); % 幅頻特性axis([0 10000 -40 1]);xlabel('頻率Hz');ylabel('幅度dB');subplot(2,1,2); plot(f, angle(H_s)); % 相頻特性xlabel('頻率Hz');ylabel('相角rad');figure(2); freqs(b,a); % 也可用指令freqs直接畫出H(s)的頻率響應曲線。

matlab常用命令大全

matlab常用命令大全

matlab命令一、常用对象操作:除了一般windows窗口的常用功能键外。

1、!dir 可以查看当前工作目录的文件。

!dir& 可以在dos状态下查看。

2、who 可以查看当前工作空间变量名,whos 可以查看变量名细节。

3、功能键:功能键快捷键说明方向上键 Ctrl+P 返回前一行输入方向下键 Ctrl+N 返回下一行输入方向左键 Ctrl+B 光标向后移一个字符方向右键 Ctrl+F 光标向前移一个字符Ctrl+方向右键 Ctrl+R 光标向右移一个字符Ctrl+方向左键 Ctrl+L 光标向左移一个字符home Ctrl+A 光标移到行首End Ctrl+E 光标移到行尾Esc Ctrl+U 清除一行Del Ctrl+D 清除光标所在的字符Backspace Ctrl+H 删除光标前一个字符Ctrl+K 删除到行尾Ctrl+C 中断正在执行的命令4、clc可以命令窗口显示的内容,但并不清除工作空间。

二、函数及运算1、运算符:+:加,-:减, *:乘, /:除,\:左除 ^:幂,‘:复数的共轭转置,():制定运算顺序。

2、常用函数表:sin( ) 正弦(变量为弧度)Cot( ) 余切(变量为弧度)sind( ) 正弦(变量为度数)Cotd( ) 余切(变量为度数)asin( ) 反正弦(返回弧度)acot( ) 反余切(返回弧度)Asind( ) 反正弦(返回度数)acotd( ) 反余切(返回度数)cos( ) 余弦(变量为弧度)exp( ) 指数cosd( ) 余弦(变量为度数)log( ) 对数acos( ) 余正弦(返回弧度)log10( ) 以10为底对数acosd( ) 余正弦(返回度数)sqrt( ) 开方tan( ) 正切(变量为弧度)realsqrt( ) 返回非负根tand( ) 正切(变量为度数)abs( ) 取绝对值atan( ) 反正切(返回弧度)angle( ) 返回复数的相位角atand( ) 反正切(返回度数)mod(x,y) 返回x/y的余数sum( ) 向量元素求和3、其余函数可以用help elfun和help specfun命令获得。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

微环谐振器 matlab 代码以下是微环谐振器的 matlab 代码:
%%微环谐振器
%%1.定义信号的频率范围及步长
f_start = 1e7;
f_end = 1.2e7;
f_step = 1e4;
%%2.定义微环谐振器的参数
diameter = 10e-6; %微环直径
n = 1; %环数
g = 200e-9; %微环与波导之间的间隙
neff = 2.3; %波导的等效折射率
alpha = 0.1; %微环的损耗
L = 100e-6; %微环的长度
%%3.计算微环谐振器的参数
lambda = 1550e-9; %波长
radius = diameter/2; %微环半径
circumference = 2*pi*radius; %微环周长
beta = 2*pi*n/circumference; %微环的模式 delta = beta*neff*2*pi/L; %微环的色散
Q = pi*radius*neff/alpha; %微环的品质因数 %%4.计算微环谐振器的响应
f = f_start:f_step:f_end;
w = 2*pi*f;
S = Q./(1+1i*Q.*(w/(beta*neff)-1/delta)); %%5.绘制微环谐振器的响应图像
figure;
plot(f/1e6,20*log10(abs(S)));
xlabel('频率(MHz)');
ylabel('响应(dB)');
title('微环谐振器的响应');
grid on;
可以根据需要进行进一步的修改和优化。

相关文档
最新文档