PSCAD输出数据并用matlab调入的方法

1、PSCAD中右键对应项目“Project Setting”，在Runtime选项卡下

然后在该项目的工作目录下找到.emt文件夹，找到.out和.inf文件

2、将两个文件移到matlab工作目录下，用记事本打开.out文件，删除第一行的文字

Matlab程序中用load(‘XXX.out’)可以载入仿真数据，第一列是时间，第二列开始的数据，由.inf 文档说明

MATLAB数据分析与多项式计算(M)

第7章 MATLAB数据分析与多项式计算 6.1 数据统计处理 6.2 数据插值 6.3 曲线拟合 6.4 离散傅立叶变换 6.5 多项式计算 6.1 数据统计处理 6.1.1 最大值和最小值 MATLAB提供的求数据序列的最大值和最小值的函数分别为max 和min，两个函数的调用格式和操作过程类似。 1．求向量的最大值和最小值 求一个向量X的最大值的函数有两种调用格式，分别是： (1) y=max(X)：返回向量X的最大值存入y，如果X中包含复数元素，则按模取最大值。 (2) [y,I]=max(X)：返回向量X的最大值存入y，最大值的序号存入I，如果X中包含复数元素，则按模取最大值。 求向量X的最小值的函数是min(X)，用法和max(X)完全相同。 例6-1 求向量x的最大值。 命令如下： x=[-43,72,9,16,23,47]; y=max(x) %求向量x中的最大值 [y,l]=max(x) %求向量x中的最大值及其该元素的位置 2．求矩阵的最大值和最小值 求矩阵A的最大值的函数有3种调用格式，分别是： (1) max(A)：返回一个行向量，向量的第i个元素是矩阵A的第i 列上的最大值。 (2) [Y,U]=max(A)：返回行向量Y和U，Y向量记录A的每列的最大值，U向量记录每列最大值的行号。 (3) max(A,[],dim)：dim取1或2。dim取1时，该函数和max(A)完全相同；dim取2时，该函数返回一个列向量，其第i个元素是A矩阵的第i行上的最大值。 求最小值的函数是min，其用法和max完全相同。

例6-2 分别求3×4矩阵x中各列和各行元素中的最大值，并求整个矩阵的最大值和最小值。 3．两个向量或矩阵对应元素的比较 函数max和min还能对两个同型的向量或矩阵进行比较，调用格式为： (1) U=max(A,B)：A,B是两个同型的向量或矩阵，结果U是与A,B 同型的向量或矩阵，U的每个元素等于A,B对应元素的较大者。 (2) U=max(A,n)：n是一个标量，结果U是与A同型的向量或矩阵，U的每个元素等于A对应元素和n中的较大者。 min函数的用法和max完全相同。 例6-3 求两个2×3矩阵x, y所有同一位置上的较大元素构成的新矩阵p。 6.1.2 求和与求积 数据序列求和与求积的函数是sum和prod，其使用方法类似。设X是一个向量，A是一个矩阵，函数的调用格式为： sum(X)：返回向量X各元素的和。 prod(X)：返回向量X各元素的乘积。 sum(A)：返回一个行向量，其第i个元素是A的第i列的元素和。 prod(A)：返回一个行向量，其第i个元素是A的第i列的元素乘积。 sum(A,dim)：当dim为1时，该函数等同于sum(A)；当dim为2时，返回一个列向量，其第i个元素是A的第i行的各元素之和。 prod(A,dim)：当dim为1时，该函数等同于prod(A)；当dim为2时，返回一个列向量，其第i个元素是A的第i行的各元素乘积。 例6-4 求矩阵A的每行元素的乘积和全部元素的乘积。 6.1.3 平均值和中值 求数据序列平均值的函数是mean，求数据序列中值的函数是median。两个函数的调用格式为： mean(X)：返回向量X的算术平均值。 median(X)：返回向量X的中值。

将PSCAD中的数据导入MATLAB

如何将PSCAD/EMTDC中的数据导入MATLAB中呢？ 以接地极线路单线接地故障（将模型命名为WLDanjiedi01）为例进行详细的介绍：1、模型建立完毕，右击选择“Project Settings”出现如下界面 将”Save channels to disk?”选择为“Yes”,并在后面的“Output file”进行输出文件的命名，如例文件名命名为“WLDanjiedi01.out”（最好与模型名称一致），将模型保存至XX位置。 2、模型仿真完毕，在XX位置会生成一个名为“WLDanjiedi01.emt”的文件夹， 文件夹中后缀为“WLDanjiedi01-01.out到WLDanjiedi01-06.out”的文件储存着仿真所得到的数据；名为“WLDanjiedi01.inf”的文件是所有数据的说明，如果需要在MATLAB中进行编程处理数据，则要根据此文件中的说明在MATLAB中进行变量的定义。 3、在MATLAB中的工作窗口如下，

单击“Import data”找到“WLDanjiedi01.emt”目录，界面如下 下拉文件类型（T）选择“All Files(*.*)”出现如下界面 选择“WLDanjiedi01-01.out到WLDanjiedi01-06.out”中所需要的即可，例如导入“WLDanjiedi01-01.out”，选中后点击打开，经过一定时间会出现如下界面

选择“Next”，接着选择“Finish”即可完成数据的导入，此时MATLAB中的工作窗口如下，出现了“WLDanjiedi01-01”文件夹。 选中“WLDanjiedi01-01”，界面变成如下，单击“Plot（WLDanjiedi01-01）”会生成此文件夹所包含数据的波形图。

数值计算方法实验指导(Matlab版)

《数值计算方法》实验指导 （Matlab 版） 肇庆学院数学与统计学学院 计算方法课程组

1. 实验名称 实验1 算法设计原则验证(之相近数相减、大数吃小数和简化计算步骤) 2. 实验题目 有效数字的损失. 123 )与1000个较小的数(3 10 15)的和，验证 大数吃小数的现象. (3)分别用直接法和秦九韶算法计算多项式 P(x) a 0x n a 1x n 1 在x =1.00037 处的值?验证简化计算步骤能减少运算时间. n 1 对于第(3)题中的多项式P (x )，直接逐项计算需要n (n 1) 2 1 次乘法 和n 次加法，使用秦九韶算法 P(x) (((a °x ajx a 2)x a . 则只需要n 次乘法和n 次加法. 3. 实验目的 验证数值算法需遵循的若干规则. 4. 基础理论 设计数值算法时，应避免两个相近的数相减、防止大数吃小数、简化计算步骤减少运算 次数以减少运算时间并降低舍入误差的积累. 两相近的数相减会损失有效数字的个数， 用一 《数值计算方法》实验 1报告 班级： 20xx 级 XXXXx 班 学号： 20xx2409xxxx 姓名： XXX 成绩: ⑴取 z 1016，计算z 1 Z 和 1/（、z 1 Z)，验证两个相近的数相减会造成 (2)按不同顺序求一个较大的数( a n 1 X a n

个大数依次加小数，小数会被大数吃掉，乘法运算次数太多会增加运算时间. 5.实验环境 操作系统：Win dows xp ；程序设计语言：Matlab 6.实验过程 (1)直接计算并比较； (2)法1 :大数逐个加1000个小数，法2 :先把1000个小数相加再与大数加； (3)将由高次项到低次项的系数保存到数组A[n]中，其中n为多项式次数. 7.结果与分析 (1)计算的~1V Z = _______________________________ ,1/( ~1 < z) ____________________ . 分析： (2)123逐次加1000个3 10 6的和是_________________________ ，先将1000个3 10 6相 加，再用这个和与123相加得_______________________ . 分析： (3)计算__________ 次的多项式： 直接计算的结果是___________________ ，用时___________________ ; 用秦九韶算法计算的结果是____________________ ，用时 ___________________ 分析：

MATLAB导入CAD数据

用AutoCAD绘制平面公式曲线（如渐开线、心形线）、空间公式曲线（如螺旋线）以及公式曲面（如马鞍形曲面）是比较困难的，一般情况下，需要用AutoCAD开发程序编程，但多数程序比较复杂，尤其是公式曲面的绘制程序，需要多层嵌套循环，复杂且运行效率低。 快速且精确地绘制各种公式曲线、曲面恰恰是MATLAB的长项，但是MATLAB绘制的图形却不能直接用于机械零件设计。其中非常关键的一点，就是MATLAB绘制的曲线、曲面分别是由有限个点连接而成的折线和空间网格构成的，而在AutoCAD中绘制的曲线、曲面也是如此。因此，只需要把在MATLAB中绘制的公式曲线、曲面上所有的点坐标数据都提取出来，若能让AutoCAD正确识别，那么我们就可以在AutoCAD中精确地绘制这些曲线、曲面了。 本文介绍了一种快速、精确地绘制各种公式曲线、曲面的方法，即在AutoCAD中通过调用经过Excel处理的MATLAB数据实现。 二、AutoCAD和MATLAB的特点 MATLAB是非常优秀的科学计算、信号处理以及图形显示软件，它有自身的语言，与其他高级语言相比，MATLAB提供了一个人机交互的数学环境，并以矩阵作为基本的数据结构，可大大节省编程时间。另外，MATLAB不仅语法规则简单，容易掌握，调试方便，还可以存储中间结果，这使得MATLAB既可以快捷、精确地绘制各种公式曲线、曲面，又可以很方便地提取中间数据。 在工业设计领域，AutoCAD不仅被广泛应用于平面绘图，也可以用于三维建模，但在曲线、曲面造型方面不是很理想。它是开放型的人机交互系统，有多种语言接口，与外界的数据交换很灵活，这些特点使得它与MATLAB的结合成为可能。 三、结合MATLAB在AutoCAD中绘制曲线、曲面的原理及方法 1.原理 MATLAB中的矩阵数据虽然很容易提取，但由于它不是AutoCAD能识别的格式，因此不能直接被AutoCAD调用，需要先用Excel对从MATLAB中提取的数据进行编辑，转换成AutoCAD可以识别的格式，才能在AutoCAD中绘出曲线、曲面。 2.方法 由于在AutoCAD中绘制平面曲线、空间曲线和曲面的绘制命令不同，且数据结构也不同，因此结合MATLAB的绘制方法也稍有区别。这种绘制方法的关键就是把数据格式转换成AutoCAD的绘制命令所需要的数据格式，只要熟悉AutoCAD的数据结构，就可以举一反三。 (1)利用MATLAB得到公式曲面数据 1)在MATLAB中绘制出曲面 在MATLAB中输入如下命令： [th,r]=meshgrid((0:5:360)*pi/180,0:.05:1); %在极坐标系下设置一个73×21的网格矩阵，即圆周方向分为73份，半径方向分为21份，总共分了1533个点，节点越多，图形越精确% [X,Y]=pol2cart(th,r); %转化为笛卡儿坐标系% Z=X+i.*Y;

数值分析Matlab作业

数值分析编程作业

2012年12月 第二章 14.考虑梯形电阻电路的设计，电路如下： 电路中的各个电流{i1,i2,…,i8}须满足下列线性方程组： 12 123 234 345 456 567 678 78 22/ 2520 2520 2520 2520 2520 2520 250 i i V R i i i i i i i i i i i i i i i i i i i i -= -+-= -+-= -+-= -+-= -+-= -+-= -+= 这是一个三对角方程组。设V=220V，R=27Ω，运用追赶法，求各段电路的电流量。Matlab程序如下： function chase () %追赶法求梯形电路中各段的电流量 a=input('请输入下主对角线向量a='); b=input('请输入主对角线向量b='); c=input('请输入上主对角线向量c='); d=input('请输入右端向量d='); n=input('请输入系数矩阵维数n='); u(1)=b(1); for i=2:n l(i)=a(i)/u(i-1); u(i)=b(i)-c(i-1)*l(i); end y(1)=d(1); for i=2:n y(i)=d(i)-l(i)*y(i-1); end x(n)=y(n)/u(n); i=n-1; while i>0 x(i)=(y(i)-c(i)*x(i+1))/u(i); i=i-1; end x 输入如下:

请输入下主对角线向量a=[0,-2,-2,-2,-2,-2,-2,-2]; 请输入主对角线向量b=[2,5,5,5,5,5,5,5]; 请输入上主对角线向量c=[-2,-2,-2,-2,-2,-2,-2,0]; 请输入方程组右端向量d=[220/27,0,0,0,0,0,0,0]; 请输入系数矩阵阶数n=8 运行结果如下： x = 8.1478 4.0737 2.0365 1.0175 0.5073 0.2506 0.1194 0.0477 第三章 14.试分别用(1)Jacobi 迭代法；(2)Gauss-Seidel 迭代法解线性方程组 1234510123412191232721735143231211743511512x x x x x ?????? ??????---????????????=--?????? --?????? ??????---?????? 迭代初始向量 (0)(0,0,0,0,0)T x =。 （1）雅可比迭代法程序如下： function jacobi() %Jacobi 迭代法 a=input('请输入系数矩阵a='); b=input('请输入右端向量b='); x0=input('请输入初始向量x0='); n=input('请输入系数矩阵阶数n='); er=input('请输入允许误差er='); N=input('请输入最大迭代次数N='); for i=1:n for j=1:n if i==j d(i,j)=a(i,j); else d(i,j)=0; end end end m=eye(5)-d\a; %迭代矩阵 g=d\b; x=m*x0+g; k=1; while k<=N %进行迭代 for i=1:5 if max(abs(x(i)-x0(i))) >er x=m*x+g; k=k+1;

导入包含数据的txt文件到MATLAB中

导入包含数据的txt文件到MATLAB中，并绘制图像 (2013-08-07 17:14:49) 转载▼ 标签： matlab 这回把步骤写得详详细细的，再不会忘记了吧，哇呀呀哎呀 第一步：先把txt文件复制到MATLAB的目录，或者在MATLAB中将路径指向txt文件所在路径。 第二步：右键存有数据的txt文件，选择Import Data... 第三步：Import Data之后就能看到txt里的数据被妥善安放好位置了，然后在Range右边的列表中选择Matrix，再点击绿色的对勾√导入数据：

第四步：导入完数据后，在workplace里能看到名为txt文件名的数组变量，就说明导入成功，这里是a：

第五步:最后就是编写语句了：plot(a(:,2),a(:,3),'o')，回车就会出现以o为点的散点图，如果是：plot(a(:,2),a(:,3),'*')，就得到以*为点的散点图；

绘图说明(本节来自互联网资源)： 1.将数据表的各列数值分别赋予变量x、y、z等，格式如下：x=sheetname(:,1), y=sheetname(:,2), z=sheetname(:,3); 2.用命令plot(x,y,’XXXX’)绘制图形，单引号中的符号表示点线的属性，如线形、颜色、点的形状等，若用双对数坐标画图则命令为loglog(x,y)； 3.在弹出的绘图界面中用菜单View—Property Editor编辑图形属性，如字体大小、数据点形状、横纵坐标名称、绘图区域颜色等； 4.绘图方法2：在数组编辑器上点击Plot Selection按钮，选择图形的类型即可； 5.绘图方法3：菜单File—New—Figure创建新的图形，在图形编辑器中Figure Palette面板点击2D Axes，点击右下角Add Data选择图表类型和坐标轴的数据源，度分布图将坐标轴由线形改为对数即可。 6.hold on/off命令：叠绘命令，切换绘图的保持功能； 7.绘制双纵轴： 7.1 plotyy(x1,y1,x2,y2)：分别用左/右侧y轴表示两条曲线； 7.2 plotyy(x1,y1,x2,y2,FUN)：FUN是字符串格式，用来指定绘图的函数名，可以由多个。

数值分析的matlab实现

第2章牛顿插值法实现 参考文献：[1]岑宝俊. 牛顿插值法在凸轮曲线修正设计中的应用[J]. 机械工程师,2009,10:54-55. 求牛顿插值多项式和差商的MA TLAB 主程序： function[A,C,L,wcgs,Cw]=newpoly(X,Y) n=length(X);A=zeros(n,n);A(:,1) =Y'; s=0.0;p=1.0;q=1.0;c1=1.0; for j=2:n for i=j:n A(i,j)=(A(i,j-1)-A(i-1,j-1))/(X(i)-X(i-j+1)); end b=poly(X(j-1));q1=conv(q,b);c1=c1*j;q=q1; end C=A(n,n);b=poly(X(n));q1=conv(q1,b); for k=(n-1):-1:1 C=conv(C,poly(X(k)));d=length(C);C(d)=C(d)+A(k,k); end L(k,:)=poly2sym(C);Q=poly2sym(q1); syms M wcgs=M*Q/c1;Cw=q1/c1; （1）保存名为newpoly.m 的M 文件 （2）输入MA TLAB 程序 >> X=[242,243,249,250]; >> Y=[13.681,13.526,13.098,13.095]; >> [A,C,L,wcgs,Cw]=newpoly(X,Y) 输出3阶牛顿插值多项式L 及其系数向量C 差商的矩阵A ，插值余项wcgs 及其 ) ()()1(ξ+n n f x R 的系数向量Cw 。 A = 13.6810 0 0 0 13.5260 -0.1550 0 0 13.0980 -0.0713 0.0120 0 13.0950 -0.0030 0.0098 -0.0003 C = 1.0e+003 *

Matlab文件操作及读txt文件(fopen,fseek,fread,fclose)

Matlab文件操作及读txt文件(fopen,fseek,fread,fclose) matlab文件操作 文件操作是一种重要的输入输出方式，即从数据文件读取数据或将结果写入数据文件。MATLAB提供了一系列低层输入输出函数，专门用于文件操作。 1、文件的打开与关闭 1）打开文件 在读写文件之前，必须先用fopen函数打开或创建文件，并指定对该文件进行的操作方式。fopen函数的调用格式为： fid=fopen（文件名，…打开方式?） 说明：其中fid用于存储文件句柄值，如果返回的句柄值大于0，则说明文件打开成功。文件名用字符串形式，表示待打开的数据文件。常见的打开方式如下： λ…r?：只读方式打开文件（默认的方式），该文件必须已存在。 …r+?：读写方式打开文件，打开后先读后写。该文件必须已存在。λλ…w?：打开后写入数据。该文件已存在则更新；不存在则创建。 …w+?：读写方式打开文件。先读后写。该文件已存在则更新；不存在则创建。λ λ…a?：在打开的文件末端添加数据。文件不存在则创建。

…a+?：打开文件后，先读入数据再添加数据。文件不存在则创建。 另外，在这些字符串后添加一个“t”，如…rt?或…wt+?，则将该文件以文本方式打开；如果添加的是“b”，则以二进制格式打开，这也是fopen 函数默认的打开方式。 2）关闭文件 文件在进行完读、写等操作后，应及时关闭，以免数据丢失。关闭文件用fclose函数，调用格式为： sta＝fclose(fid) 说明：该函数关闭fid所表示的文件。sta表示关闭文件操作的返回代码，若关闭成功，返回0，否则返回-1。如果要关闭所有已打开的文件用fclose(…all?)。 2、二进制文件的读写操作 1）写二进制文件 fwrite函数按照指定的数据精度将矩阵中的元素写入到文件中。其调用格式为： COUNT＝fwrite（fid，A，precision） 说明：其中COUNT返回所写的数据元素个数（可缺省），fid为文件句柄，A用来存放写入文件的数据，precision代表数据精度，常用的数据

Matlab的各种数据读取、文件读写等操作汇总

Matlab 的各种数据读取、文件读写等操作汇总 MATLAB 提供了多种方式从磁盘读入文件或将数据输入到工作空间，即读取数据，又叫导入数据；将工作空间的变量存储到磁盘文件中称为存写数据，又叫导出数据。至于选择哪种机制，则根据下面两个因素决定：?用户所执行的 操作是导入数据还是导出数据；?数据的格式为文本格式、 二进制格式还是如HDF 之类的标准格式。将数据导入MATLAB 中最容易的方法就是使用导入数据模板(Import Wizard) ，使用该模板时不需要知道数据的格式，只需指定包含这些数据的文件，然后导入模板会自动处理文件内容。本章重点内容如下：? 文件的打开和关闭? 文本文件的读取?存写ASCII数据?二进制数据的读取? 二进制数据的存写? 使用I/O文件函数进行数据读写?MAT 文件的读写 2.1 文件的打开和关闭2.1.1 文件的打开无论是要读写ASCII 码文件还是二进制文件，都必须先用fopen 函数将其打开，在默认情况下，fopen 以二进制格式打开文件，它的使用语法如下：fopen ('filename', 'mode') 其中filename 表示要读写的文件名称，mode 则表示要对文件进行的处理方式，如下：rt :以只读方式(Reading)打开文件wt:以只写方式(Writing)打开文件at:以追加方式(Appending)打开文件，新内容将从原文件后面续写r+t:以同时读写方式打开文件w+t ：以同时读写创建文件，原文件内容被清除

a+t ：以同时读和追加(Reading and Appdending) 方式，原文件内容被保留，新内容将从原文件的后面开始At ：以读写方式打开或创建文件，适用于对磁带介质文件的操作Wt ：以写入方式打 开或创建文件，原文件内容被清除，适用于磁带介质文件的操作fopen 函数有两个返回值，一个是返回一个文件标志(file Identifier) ，它会作为参数被传入其他对文件进行读写操作的命令，通常是一个非负的整数，可用此标识来对此文件进行各种处理。如果返回的文件标识是-1，则代表fopen无法打开文件，其原因可能是文件不存在，或是用户无法打开此文件权限。另一个返回值就是message ，用于返回无法打开文件的原因。为了安全起见，最好在每次使用fopen 函数时，都测试其返回值是否为有效值。下面以脚本m 文件为例来声明文件的打开。例 2-1 %exam1.m[f,message]=fopen('fileexam1', 'r')if f==-1disp (message); % 显示错误信息end 若文件fileexam1 不存在，则显示如下信息。Cannot open file.existence?permissions?memory?... 例2-2 %exam2.m[f,message]=fopen('fileexam2', 'r');if f==-1disp (message); % 显示错误信息else disp(f);end 若文件fileexam2 存在，则返回f值。 2.1.2 文件的关闭一旦完成文件的读写，最好关闭文件，以便对其进行其他操作。这时就可以使用fclose 函数来关闭文件，其适用语法如下：fclose(f) 。其中 f 为打开文件的标志，若fclose 函数返回值为0 ，则表示成功关闭 f 标志的文件；若返回值为-1，

第2讲 matlab的数值分析

第二讲MATLAB的数值分析 2-1矩阵运算与数组运算 矩阵运算和数组运算是MATLAB数值运算的两大类型，矩阵运算是按矩阵的运算规则进行的，而数组运算则是按数组元素逐一进行的。因此，在进行某些运算（如乘、除）时，矩阵运算和数组运算有着较大的差别。在MATLAB中，可以对矩阵进行数组运算，这时是把矩阵视为数组，运算按数组的运算规则。也可以对数组进行矩阵运算，这时是把数组视为矩阵，运算按矩阵的运算规则进行。 1、矩阵加减与数组加减 矩阵加减与数组加减运算效果一致，运算符也相同，可分为两种情况： （1）若参与运算的两矩阵（数组）的维数相同，则加减运算的结果是将两矩阵的对应元素进行加减，如 A=[1 1 1;2 2 2;3 3 3]; B=A; A+B ans= 2 2 2 4 4 4 6 6 6 （2）若参与运算的两矩阵之一为标量（1*1的矩阵），则加减运算的结果是将矩阵（数组）的每一元素与该标量逐一相加减，如 A=[1 1 1;2 2 2;3 3 3]; A+2 ans= 3 3 3 4 4 4 5 5 5 2、矩阵乘与数组乘 （1）矩阵乘 矩阵乘与数组乘有着较大差别，运算结果也完全不同。矩阵乘的运算符为“*”，运算是按矩阵的乘法规则进行，即参与乘运算的两矩阵的内维必须相同。设A、B为参与乘运算的 =A m×k B k×n。因此，参与运两矩阵，C为A和B的矩阵乘的结果，则它们必须满足关系C m ×n 算的两矩阵的顺序不能任意调换，因为A*B和B*A计算结果很可能是完全不一样的。如：A=[1 1 1;2 2 2;3 3 3]; B=A;

A*B ans= 6 6 6 12 12 12 18 18 18 F=ones(1,3); G=ones(3,1); F*G ans 3 G*F ans= 1 1 1 1 1 1 1 1 1 （2）数组乘 数组乘的运算符为“.*”，运算符中的点号不能遗漏，也不能随意加空格符。参加数组乘运算的两数组的大小必须相等（即同维数组）。数组乘的结果是将两同维数组（矩阵）的对应元素逐一相乘，因此，A.*B和B.*A的计算结果是完全相同的，如： A=[1 1 1 1 1;2 2 2 2 2;3 3 3 3 3]; B=A; A.*B ans= 1 1 1 1 1 4 4 4 4 4 9 9 9 9 9 B.*A ans= 1 1 1 1 1 4 4 4 4 4 9 9 9 9 9 由于矩阵运算和数组运算的差异，能进行数组乘运算的两矩阵，不一定能进行矩阵乘运算。如 A=ones(1,3); B=A; A.*B ans= 1 1 1 A*A ???Error using= =>

实验1 MATLAB使用方法和程序设计

实验1 MATLAB 使用方法和程序设计 一、实验目的 1、掌握MATLAB 软件使用的基本方法。 2、熟悉MATLAB 的数据表示、基本运算和程序控制语句。 3、熟悉MATLAB 绘图命令及基本绘图控制。 4、熟悉MATLAB 程序设计的基本方法。 二、实验内容 1.帮助命令 使用Help 命令，查找sqrt （开方）函数的使用方法。 2、矩阵运算 （1）矩阵乘法 已知A=[1 2;3 4];B=[5 5;7 8]; 求A^2*B 。 (2) 矩阵除法 已知A=[1 2 3;4 5 6;7 8 9];B=[1 0 0;0 2 0;0 0 3]; 求A\B,A/B 。 (3) 矩阵的转置及共轭转置 已知A=[5+i ，2-i ，1；6*i ，4,9-i]；求A.’,A ’ 。 （4）使用冒号选出指定元素 已知A=[1 2 3 ; 4 5 6 ; 7 8 9]; 求A 中第3列前2个元素；A 中所有列第2,3行的元素。 （5）方括号[] 用magic 函数生成一个4阶魔术矩阵，删除该矩阵的第四列。 3、多项式 （1）求多项式p(x)=x 3-2x-4的根。 （2）求 f (x) = (cos x)2 的一次导数。 （3）求微分方程 的通解，并验证。 4、基本绘图命令 （1）绘制余弦曲线cost =y , ]2,0[π∈t 。 （2）在同一坐标系中绘制余弦曲线 0.25)-cos(t =y 和正弦曲线 )5.0sin(-=t y ，]2,0[π∈t （3）用plot3函数绘制三维螺线： 22x dy xy xe dx -+= sin()cos()x t y t z t =??=??=?( 0 < t < 20 )

MATLAB数据输入和输出 C语言转换

数据输入和输出 一、概述 MATLAB可读和写的文件格式： 文件格式文件内容扩展名读、写文件的函数MATLAB 保存的MA TLAB工作区.mat load，save 文本Text 任意textscan，textread 限定文本任意dlmread，dlmwrite 逗号分隔的数值.csv csvread，csvwrite XML XML-格式文本t .xml xmlread，xmlwrite Audio NeXT/SUN sound .au auread，auwrite Microsoft WA VE sound .mav mavread，mavwrite Movie Audio/video .avi aviread Scientific data Data in Common Data Format .cdf cdfread，cdfwrite Flexible Image Transport system data .fits fitsread Data in Hierarchical Data Format .hdf hdfread Spreadsheet Excel worksheet .xls xlsread，xlswrite Lotus 123 worksheet .wk1 wk1read，wk1write Graphics TIFF/PNG/HDF/BMP/JPEG image GIF/PCX/XWD/Cursor/Icon image .tiff .png .hdf .bmp .jpeg .gif .pcx .xwd .cur .ico imread，imwrite 二、使用输入向导（Import Wizard） 从菜单File->Import Data打开Import Wizard；或者命令窗口输入函数uiimport。 从剪贴板开始Import Wizard：Edit->Paste to workspace。 三、保存和加载MAT文件 MAT文件是双精度、二进制、MATLAB格式的文件。 输出到MA T文件：save filename [var1 var2 …] [str*]；可以通过[var1 var2 …]选择性保存变量；也可以使用通配符“*”。 查看MAT文件中的变量：whos –file 文件名。 存放结构数组的某个字段：加上“-struct”选项。 在已经存在的MAT文件上添加数据：-append选项。 禁止压缩和Unicode字符编码文件，在save语句中加入“-v6”或者File->Preferences-> General->MA T-Files->MAT-File save options->Ensure backward compatibility(-v6)。save语句默认为数据压缩。 选择输出格式：“-ascii”、“-tabs”、“-double”“-v4”。 从MA T文件输入数据：load函数。 四、输入文本数据 各种输入函数的数据定界符： csvread 仅逗号dlmread 任何字符（可指定）fscanf 任何字符load 仅空格 textread 任何字符textscan 任何字符详见：doc 函数名。 textscan和textread性能比较：前者有更好的性能，特别是读大文件时；使用前者首先要

数值分析的MATLAB程序

列主元法 function lianzhuyuan(A,b) n=input('请输入n:') %选择阶数A=zeros(n,n); %系数矩阵A b=zeros(n,1); %矩阵b X=zeros(n,1); %解X for i=1:n for j=1:n A(i,j)=(1/(i+j-1)); %生成hilbert矩阵A end b(i,1)=sum(A(i,:)); %生成矩阵b end for i=1:n-1 j=i; top=max(abs(A(i:n,j))); %列主元 k=j; while abs(A(k,j))~=top %列主元所在行 k=k+1; end for z=1:n %交换主元所在行a1=A(i,z); A(i,z)=A(k,z); A(k,z)=a1; end a2=b(i,1); b(i,1)=b(k,1); b(k,1)=a2; for s=i+1:n %消去算法开始m=A(s,j)/A(i,j); %化简为上三角矩阵 A(s,j)=0; for p=i+1:n A(s,p)=A(s,p)-m*A(i,p); end b(s,1)=b(s,1)-m*b(i,1); end end X(n,1)=b(n,1)/A(n,n); %回代开始 for i=n-1:-1:1 s=0; %初始化s for j=i+1:n s=s+A(i,j)*X(j,1);

end X(i,1)=(b(i,1)-s)/A(i,i); end X 欧拉法 clc clear % 欧拉法 p=10; %贝塔的取值 T=10; %t取值的上限 y1=1; %y1的初值 r1=1; %y2的初值 %输入步长h的值 h=input('欧拉法please input number(h=1 0.5 0.25 0.125 0.0625):h=') ; if h>1 or h<0 break end S1=0:T/h; S2=0:T/h; S3=0:T/h; S4=0:T/h; i=1; % 迭代过程 for t=0:h:T Y=(exp(-t)); R=(1/(p-1))*exp(-t)+((p-2)/(p-1))*exp(-p*t); y=y1+h*(-y1); y1=y; r=r1+h*(y1-p*r1); r1=r; S1(i)=Y; S2(i)=R; S3(i)=y; S4(i)=r; i=i+1; end t=[0:h:T]; % 红线为解析解，'x'为数值解 plot(t,S1,'r',t,S3,'x')

MATLAB数据输入和输出

数据输入和输出 一、概述 二、使用输入向导（Import Wizard） 从菜单File->Import Data打开Import Wizard；或者命令窗口输入函数uiimport。 从剪贴板开始Import Wizard：Edit->Paste to workspace。 三、保存和加载MAT文件 MA T文件是双精度、二进制、MA TLAB格式的文件。 输出到MA T文件：save filename [var1 var2 …] [str*]；可以通过[var1 var2 …]选择性保存变量；也可以使用通配符“*”。 查看MA T文件中的变量：whos –file 文件名。 存放结构数组的某个字段：加上“-struct”选项。 在已经存在的MA T文件上添加数据：-append选项。 禁止压缩和Unicode字符编码文件，在save语句中加入“-v6”或者File->Preferences-> General->MA T-Files->MA T-File save options->Ensure backward compatibility(-v6)。save语句默认为数据压缩。 选择输出格式：“-ascii”、“-tabs”、“-double”“-v4”。 从MA T文件输入数据：load函数。 四、输入文本数据 各种输入函数的数据定界符： textscan和textread性能比较：前者有更好的性能，特别是读大文件时；使用前者首先要

打开文件，最后要关闭文件，可从文件任意位置读；前者只输出一个单元数组，不必给每个被读字段指定一个输出参数；前者有更多的数据转换选项和更多的用户设置选项。 五、输出文本数据 六、输入/输出标准图像文件 七、输入/输出音频和视频数据 八、输入/输出电子表数据 九、低级文件输入/输出函数

第6章matlab数据分析与多项式计算_习题答案

第6章 MATLAB数据分析与多项式计算 习题6 一、选择题 1．设A=[1,2,3,4,5;3,4,5,6,7]，则min(max(A))的值是（）。B A．1 B．3 C．5 D．7 2．已知a为3×3矩阵，则运行mean(a)命令是（）。B A．计算a每行的平均值 B．计算a每列的平均值 C．a增加一行平均值 D．a增加一列平均值 3．在MATLAB命令行窗口输入下列命令： >> x=[1,2,3,4]; >> y=polyval(x,1); 则y的值为（）。 D A．5 B．8 C．24 D．10 4．设P是多项式系数向量，A为方阵，则函数polyval(P,A)与函数polyvalm(P,A)的值（）。D A．一个是标量，一个是方阵 B．都是标量 C．值相等 D．值不相等 5．在MATLAB命令行窗口输入下列命令： >> A=[1,0,-2]; >> x=roots(A); 则x(1)的值为（）。 C A．1 B．-2 C． D． 6．关于数据插值与曲线拟合，下列说法不正确的是（）。A A．3次样条方法的插值结果肯定比线性插值方法精度高。 B．插值函数是必须满足原始数据点坐标，而拟合函数则是整体最接近原始数据点，而不一定要必须经过原始数据点。 C．曲线拟合常常采用最小二乘原理，即要求拟合函数与原始数据的均方误差达到极小。 D．插值和拟合都是通过已知数据集来求取未知点的函数值。 二、填空题 1．设A=[1,2,3;10 20 30;4 5 6]，则sum(A)= ，median(A)= 。 [15 27 39]，[4 5 6[ 2．向量[2,0,-1]所代表的多项式是。2x2-1 3．为了求ax2+bx+c=0的根，相应的命令是（假定a、b、c已经赋值）。为了

Matlab数据文件的读写

Matlab数据文件的读写 在编写一个程序时，经常需要从外部读入数据，或者将程序运行的结果保存为文件。MATLAB使用多种格式打开和保存数据。本章将要介绍MATLAB中文件的读写和数据的导入导出。 了解MATLAB的基本数据操作 掌握MATLAB中文本文件的读写方式 掌握MATLAB通过界面导入导出数据 了解MATLAB中的基本输入输出函数 13.1 数据基本操作 本节介绍基本的数据操作，包括工作区的保存、导入和文件打开。 13.1.1 文件的存储 MATLAB支持工作区的保存。用户可以将工作区或工作区中的变量以文件的形式保存，以备在需要时再次导入。保存工作区可以通过菜单进行，也可以通过命令窗口进行。 1. 保存整个工作区 选择File菜单中的Save Workspace As…命令，或者单击工作区浏览器工具栏中的Save，可以将工作区中的变量保存为MAT文件。 2. 保存工作区中的变量 在工作区浏览器中，右击需要保存的变量名，选择Save

As…，将该变量保存为MAT文件。 3. 利用save命令保存 该命令可以保存工作区，或工作区中任何指定文件。该命令的调用格式如下： ● save：将工作区中的所有变量保存在当前工作区中的文件中，文件名为matlab.mat，MAT文件可以通过load函数再次导入工作区，MAT函数可以被不同的机器导入，甚至可以通过其他的程序调用。 ● save('filename')：将工作区中的所有变量保存为文件，文件名由filename指定。如果filename中包含路径，则将文件保存在相应目录下，否则默认路径为当前路径。 ● save('filename', 'var1', 'var2', ...)：保存指定的变量在filename 指定的文件中。 ● save('filename', '-struct', 's')：保存结构体s中全部域作为单独的变量。 ● save('filename', '-struct', 's', 'f1', 'f2', ...)：保存结构体s中的指定变量。 ● save('-regexp', expr1, expr2, ...)：通过正则表达式指定待保存的变量需满足的条件。 ● save('..., 'format')，指定保存文件的格式，格式可以为MAT 文件、ASCII文件等。 13.1.2 数据导入

同济大学数值分析matlab编程题汇编

MATLAB 编程题库 1.下面的数据表近似地满足函数2 1cx b ax y ++=，请适当变换成为线性最小二乘问题，编程求最好的系数c b a ,,，并在同一个图上画出所有数据和函数图像. 625 .0718.0801.0823.0802.0687.0606.0356.0995 .0628.0544.0008.0213.0362.0586.0931.0i i y x ---- 解： x=[-0.931 -0.586 -0.362 -0.213 0.008 0.544 0.628 0.995]'; y=[0.356 0.606 0.687 0.802 0.823 0.801 0.718 0.625]'; A=[x ones(8,1) -x.^2.*y]; z=A\y; a=z(1); b=z(2); c=z(3); xh=-1:0.1:1; yh=(a.*xh+b)./(1+c.*xh.^2); plot(x,y,'r+',xh,yh,'b*')

2.若在Matlab工作目录下已经有如下两个函数文件，写一个割线法程序，求出这两个函数 10 的近似根，并写出调用方式： 精度为10 解： >> edit gexianfa.m function [x iter]=gexianfa(f,x0,x1,tol) iter=0; while(norm(x1-x0)>tol) iter=iter+1; x=x1-feval(f,x1).*(x1-x0)./(feval(f,x1)-feval(f,x0)); x0=x1;x1=x; end >> edit f.m function v=f(x) v=x.*log(x)-1; >> edit g.m function z=g(y) z=y.^5+y-1; >> [x1 iter1]=gexianfa('f',1,3,1e-10) x1 = 1.7632 iter1 = 6 >> [x2 iter2]=gexianfa('g',0,1,1e-10) x2 = 0.7549 iter2 = 8

matlab文件的导入导出

GprMaxV2.0中GprMax2D输入文件的命令(1/2) 原创 修改人 修改时间 版本 微风无尘 微风无尘 2012.3.9 V1.1 实验环境： 操作系统：Windows 7 软件版本：MATLAB 7.1 & GprMaxV2.0 参考文献： GprMaxV2.0软件manual文件夹下的UserGuideV2.pdf。 GprMaxV2.0下载地址：https://www.wendangku.net/doc/264087713.html,/Download.html 说明：翻译得不好，还望大家见谅，因为我也是边看边译的。 3.1 GprMax2D命令的一般注意事项 为了描述GprMax2D/3D命令及其参数，我们作以下约定： f表示浮点数(如1.5或15e-1、1.5e1) i表示整数

c表示字符 str表示字符串 file表示输入文件名 所有空间距离的基本单位为米 所有时间的基本单位为秒 所有频率参数的基本单位是Hz 3.2 GprMax2D 2.0版本共有32条命令：#title: #domain: #dx_dy: #time_step_stability_factor: #time_window: #messages: #number_of_media: #nips_number: #media_file: #geometry_file: #medium: #abc_type:

#abc_order: #abc_stability_factors: #abc_optimization_angles: #abc_mixing_parameters: #pml_layers: #box: #cylinder: #x_segment: #y_segment: #triangle: #analysis: #end_analysis: #tx: #rx: #rx_box: #snapshot: #tx_steps: #rx_steps: #line_source: #excitation_file: