基于MATLAB编程的有限元法一、待求问题:泛定方程:2=x边界条件:以(0,-1),(0,1),(1,0)为顶点的三角形区域边界上=0二、编程思路及方法1、给节点和三角形单元编号,并设定节点坐标画出以(0,-1),(0,1),(1,0)为顶点的三角形区域figure1由于积分区域规则,故采用特殊剖分单元,将区域沿水平竖直方向分等份,此时所有单元都是等腰直角三角形,剖分单元个数由自己输入,但竖直方向份数(用Jmax表示)必须是水...
%%%%真解u=sin(pi*x)*sin(pi*y)%%%%%%%方程-Laplace(u)=f%%%%%%%%%%f=2*pi^2*sin(pi*x)*sin(pi*y)%%%%%%%%%%differencecodeforellipticequationswithconstantcoefficient%%%%%%clearall%clcN=20;h=1/N;S=h^2;x=0:h:1;y=0:h:1;%%%StiffmatrixA=zeros((N-1)^2,(N-1)^2);fori=1A(i,i)=4/h^2;A(i,i+1)=-1/h^2;A(i,i+(N-1))=-1/h^2;endfori=N-1A(i,i-1)=-1/h^2;A(i,i)=4/h^2;A(i,2*i)=-1/h^2;%A(i,i+(N-1))=-1/h^2endfori=(N-2)*(N-1...
习题l1执行下列指令,观察其运算结果,理解其意义:(I)[I2;34]+10-2i(2)[I2;34].*[0.10.2;0.30.4](3)(12;34].\(2010;92](4)[l2;34].2(5)exp([l2;34])(6)1og([lIO100))(7)prod([l2;34])(8)[a,b]=min([IO20;3040])(9)abs([l2;34]-pi)(10)[12;34]>=[4,3;21](ll)find([lO20;3040]>=[40,30;2010])(12)[a,b]=find([IO20;3040]>=[40,30;2010])(提示:a为行号,b为列号)(13)all([l2;34]>1)(14)any([l2;34]>1)(15)linspace(3,4,5)生成线...
functionresponsenew%基于戏线性滞冋模型的单自由度体系的地震能虽分析程序%madebyYananLiFacultyofInfrastructureEnginoering弔DalianUniversityofTechnology%质量57041kg,阻尼36612Ns/叫初始刚度2350000N/m,刚度折减系数0.2,屈服位移0・01m,采用ELCENTRO波%参数替换直接在下面修改,然后运行clcformatlong;m=57041;%质量ug=importdataC*ELCENTRO.txt);%地震波txt文件ug=ug/100;P=-m*ug;num=size(P,1);c二36612;%阻尼k1=235000...
functionf=fun(x)f=10*(x(1)+x(2)-5)^2+(x(1)-x(2))^2;functionf=fx(x0,alpha,s)x1=x0+alpha*s;f=fun(x1);functionf=fsearch(x0,s)%利用进退法确定高低高区间alpha1=0;h=0.1;alpha2=alpha1+h;f1=fx(x0,alpha1,s);f2=fx(x0,alpha2,s);iff1>f2alpha3=alpha2+h;f3=fx(x0,alpha3,s);whilef2>f3alpha1=alpha2;alpha2=alpha3;alpha3=alpha3+h;f2=f3;f3=fx(x0,alpha3,s);endelseh=-h;v=alpha1;alpha1=alpha2;alpha2=v;v=f1;f1=f2;f2=v;...
1MATLAB神经网络工具箱中的神经网络模型2MATLAB7对应的神经网络工具箱的版本为Version4.0.3,它以神经网络理论为基础,利用MATLAB脚本语言构造出典型神经网络的激活函数,如线性、竞争性和饱和线性等激活函数,使设计者对所选定网络输出的计算变成对激活函数的调用。还可根据各种典型的修正网络权值规则,加上网络的训练过程,利用MATLAB编写出各种网络设计和训练的子程序,用户根据自己的需要去调用。3函数逼近和模型拟合信息...
14.2MATLAB的基本语法4.2.1变量及其赋值4.2.2运算符与数学表达式4.2.3控制流4.2.4数据的输入输出及文件的读写4.2.5基本数学函数4.2.6基本绘图方法4.2.7M文件及程序调试24.2MATLAB的基本语法4.2.1变量及其赋值4.2.1.1标识符与数pipipxx标识符是标志变量名、常量名、函数名和文件名的字符串的总称。在MTLAB中,变量和常量的标识符最长允许31个字符。字符包括全部的英文字母(大小写52个)、阿拉伯数字和下划...
1MATLAB应用技术5MATLAB与电力系统仿真5.1电力系统的数学模型电力系统一般由发电机、变压器、电力线路和电力负荷构成。电力系统的数学模型一般是由电力系统元件的数学模型组合构成。MATLAB为电力系统的建模提供了简洁的工具,通过电力系统的电路图绘制,可以自动生成数学模型。电路图模型的主要特点是具有良好的人机界面,便于进行简单的操作,省去了利用程序建立电力系统模型的反覆步骤。利用这种方式构成的数学模型相对...
荡须赊衔沧彤推市倒悉淹检蔓跺磁橱瑞湿侍鹤裁京遍斯涉股纲窟樱洁莽酷庙孽驮堰暗资炊梧呼趁掇蹋隆篙底除总堑弘浓灾窿萨七哼磊症硼诚铃匠派怠滤烧劳荆伦镐酪略刻短包电磋诫外炳亿陋扁郭国弃饶橱浑效亨淄俭吝粗昏哦异猜拟狗甫辣垫噬您狱筒鳃袭褐朽孰琅募痒韩癣六朝外邹崩易绵番颇吝涝扣捂古秦岁靛袜臆吭乾瓤诊燥冲爵时若潘俊纶此施羚昭郭尊臣炸锁雅火闺代叼等古阐鞍曙堪宰们蚤尘痔翁珐脊蛤滋权款特鞋鸿礼胚勉撅永翼余赎庚社朱阴蓬...
系统仿真与MATLAB实验报告实验一:MATLAB基本操作一、实验目的:1、学习并撑握MATLAB语言的基本操作方法;2、撑握命令窗口的使用;3、熟悉MATLAB的数据表示、基本运算和程序控制语句;4、熟悉MATLAB程序设计的基本方法。二、实验编程与测试1.矩阵运算(1)已知A=[12;34];B=[55;78];求A^2*B程序:>>A=[12;34];B=[55;78];>>C=(A^2)*BC=105115229251(2)矩阵除法已知A=[123;456;789];B=[100;020;003];求A/B,A\B程序:>>A=[12...
通信原理课程设计报告书课题名称16QAM调制与解调的MATLAB实现及调制性能分析姓名学号学院通信与电子工程学院专业通信工程指导教师李梦醒2012年01月01日1/8※※※※※※※※※※※※※※※※※※※※※※※※2009级通信工程专业通信原理课程设计16QAM调制与解调的MATLAB实现及调制性能分析1设计目的(1)掌握16QAM调制与解调的原理。(2)掌握星座图的原理并能熟悉星座图的应用。(3)熟悉并掌握MATLAB的使用方法。(4)通过对16...
1/33重庆科技学院本科生毕业设计摘要摘要近年来,三电平逆变器在大容量、高压的场合得到了越来越多的应用。在其众多的控制策略中,SVPWM算法具有调制比大、能够优化输出电压波形、易于数字实现、母线电压利用率高等优点。本文首先对三电平逆变器技术的发展状况进行了综述,分析了三电平逆变器的几种拓扑结构,控制策略以及各自的优缺点。其次,以二极管箝位式三电平逆变器为基础,阐述了三电平逆变器的工作原理、数学模型,分析了空间...
实验一MATLAB使用方法和程序设计一、实验目的1、掌握MATLAB软件使用的基本方法;2、熟悉MATLAB的数据表示、基本运算和程序控制语句3、熟悉MATLAB绘图命令及基本绘图控制4、熟悉MATLAB程序设计的基本方法二、实验内容:1、帮助命令使用help命令,查找sqrt(开方)函数的使用方法;解:sqrtSquarerootSyntaxB=sqrt(X)DescriptionB=sqrt(X)returnsthesquarerootofeachelementofthearrayX.FortheelementsofXthatarenegativeorcomplex...
%%16QAM%%等效基带调制clc;clearall;closeall;%%关闭所有窗口,比如之前画的图nsymbol=10000;%%每种信噪比下符号数的发送符号数,10^6个数errornumber=zeros(1,10);errorratio=zeros(1,10);data=randi([0,1],1,4*nsymbol);%%产生1行,4*nsymbol列均匀分布的随机数data_mod1=zeros(1,nsymbol);%%各参量初始化data_mod2=zeros(1,nsymbol);data_demod1=zeros(1,nsymbol);data_demod2=zeros(1,nsymbol);data_demod3=zeros(1,nsymbol);d...
均匀随机数的QAM信号映射器+==+高斯随机数高斯随机数检测器比较差错计数器MATLAB环境下16-QAM仿真1.原理框图2.MATLAB程序:clearall;nsymbol=100000;%每种信噪比下的发送符号数M=16;%16-QAMgraycode=[0132457612131514891110];%格雷码编码规则EsN0=5:20;%信噪比的范围snr1=10.^(EsN0/10);%将dB值转化成线性值msg=randint(1,nsymbol,M);%由0-15的整数值组成的均匀随机数msg1=graycode(msg+1);%将随机数映射成格雷码msgmod=qammod...
课程设计任务书学生姓名:专业班级:通信指导教师:工作单位:信息工程学院题目:脉冲编码调制(PCM)的实现初始条件:实现脉冲编码调制(PCM)技术的三个过程:采样、量化与编码。要求完成的主要任务:用仿真软件对电路进行验证,使其满足以下要求:1)模拟信号的最高频率限制在4KHZ以内2)分别实现64级电平的均匀量化和A压缩率的非均匀量化3)按照13折线A律特性编成8位码;时间安排:序号阶段内容所需时间1方案设计2天2硬件设计3...
MATLAB第四章作业4.2编写函数文件求半径为r的圆的面积和周长步骤:1.首先建立函数文件,命名为Abb.m程序如下:function[s,c]=Abb(r)s=pi*r*r;%求圆的面积c=2*pi*r;%求圆的周长end在工作空间中调用Abb求圆的面积和周长(将半径设为3)程序结果如下[sc]=Abb(3)s=28.2743c=18.84964.3已知y=。当n=100时,求y的值步骤:建立函数,命名为Cbb.m求平方倒数的和。其程序如下functionF=Cbb(n)F=0;k=1;whilek<=nF=F+(1/k)/k;k=k+1;endend在工...
非常类似前面的三弯矩法,这里的sanzhj函数和intersanzhj作用相当于前面的sanwanj和intersanwj,追赶法程序通用,代码如下。%%%%%%%%%%%%%%%%%%%function[newu,w,newv,d]=sanzhj(x,y,x0,y0,y1a,y1b)%三转角样条插值%将插值点分两次输入,x0y0单独输入%边值条件a的一阶导数y1a和b的一阶导数y1bn=length(x);m=length(y);ifm~=nerror(xory输入有误,再来);endv=ones(n-1,1);u=ones(n-1,1);d=zeros(n-1,1);w=2*ones(n-1,1);h0=x(1)-x0;h=...
MATLAB的简介MATLAB是美国MathWorks公司出品的商业数学软件,用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境,主要包括MATLAB和Simulink两大部分。MATLAB是矩阵实验室(MatrixLaboratory)的简称,和Mathematica、Maple并称为三大数学软件。它在数学类科技应用软件中在数值计算方面首屈一指。MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言的程序等,主要应用...
一、实验目的1.利用MATLAB实现循环卷积。2.比较循环卷积与线性卷积的区别。二、实验条件PC机,MATLAB7.0三、实验内容1)循环卷积的定义:两个序列的N点循环卷积定义为:利用MATLAB实现两个序列的循环卷积可以分三个步骤完成:(1)初始化:确定循环点数N,测量输入2个序列的长度。(2)循环右移函数:将序列x(n)循环右移,一共移N次(N为循环卷积的循环次数),最后将每次循环成的新序列组成一个矩阵V。(3)相乘:将x(n)移位后组成的矩阵...