gsls200808的专栏

1、 输入一组整数a，输出其中奇偶数、奇偶数的和、积以及个数。clearclcn=input('输入数字个数');for i=1:n x(i)=input('输入数字：');endj=1;k=1;%y向量存奇数，z向量存偶数for i=1:n if mod(x(i),2) y(j)=x(i); j=j+1; els

Matlab帮助文档 ^ Matrix power. Z = X^y is X to the y power if y is a scalar and X is square. If y is an integer greater than one, the power is computed by repeated squaring. For other valu

最近有位童鞋问笔者下面两段Matlab代码执行结果为啥不一样代码1x=sym('[0 1;-1 0]')y=sym('[1 -1;-2 1]')x*y代码2syms x ysubs(x*y,{x,y},{[0 1;-1 0],[1 -1;-2 1]})代码1执行结果x = [ 0, 1][ -1, 0] y = [ 0, 0][ -2, 1]

对数函数易错点函数名功能易错点示例log2以2为底的对数易写成ln2log210log_210写法为log2(10)log以e为底的对数易写成lnln10ln10写法为log(10)log10以e为底的对数易写成ln10,loglog1035log_{10}35写法为log10(35)指数函数易错点函数名功能易错点示例exp以e为底的指数易写成e^xexe^x写法为exp(x)取整函

均差定义 若已知函数f(x)f(x)在点d0,x1,...xnd_0,x_1,...x_n处的值f(x0),f(x1),...f(xn).f(x_0),f(x_1),...f(x_n).如果i≠j,i≠j,则 一阶均差f[xj,xj+1]=f(xj+1)−f(xj)xj+1−xj(j=0,1,...n−1)f[x_j,x_{j+1}]=\frac{f(x_{j+1})-f(x_j)}{x_{j+1

加工奶制品的生产计划 加工1桶牛奶有两种加工方式，(1)生产12小时得到3公斤A1，获利24元/公斤(2)生产8小时得到4公斤A2，获利16元/公斤 每天： 50桶牛奶 时间480小时 至多加工100公斤A1 制订生产计划，使每天获利最大分析： 设x1桶牛奶生产A1x1桶牛奶生产A2 则生产A1获利 24×3x1 生产A2获利 16×4 x2 每天获利 Max z=72x1+64

符号解用dsolve函数 例 求微分方程y′′+4y′+4y=e−2xy''+4y'+4y=e^{-2x} dsolve('D2y+4*Dy+4*y=exp(-2*x)')结果ans =1/(4*exp(2*x)) + C5/exp(2*t) + (C6*t)/exp(2*t)例 求解微分方程x2y′+xy=y2,y(1)=1x^2y'+xy=y^2,y(1)=1的特解dsolve('x^2*Dy

1．用2种方式创建符号表达式f=sin(x)+e^x方法一clearclcx=sym('x');f=sin(x)+exp(x)结果 f = exp(x) + sin(x) >> 方法二clearclcsyms x;f=sin(x)+exp(x)结果 f = exp(x) + sin(x) >> 2．计算习题 1 中表达式在x=0、x=

1.符号极限limit(F,x,a) 对表达式F，x -> a时的极限limit(F,a)对上面函数的简化，对表达式F，默认自变量 -> a时的极限，这里的默认自变量可以通过findsym(F,1)查找出来，一般是离x字母最近的自变量。例如：clearclcsyms a b x y;f=a+2*b+3*x+4*y;limit(f,2)findsym(f,1)结果ans

数值插值一维插值YI=interp1(X,Y,XI,METHOD) X、Y为采样点和样本值，XI为表示要插值的点，METHOD默认是线性插值'linear'METHOD可以取'nearest'  - 邻近点插值法'linear'   - 线性插值'spline'   - 三次样条插值'pchip'    - 立方插值法'cubic'    - same as '

由于这篇博文公式实在太多，截图又不太好看，所以切换到Markdown编辑器下来写。还好优快云的这个Markdown支持LaTeX公式，方便许多了。 1.计算下列各式： (1)limn→0tanx−sinx1−cos2x\lim_{n \to 0} \frac{\tan x - \sin x}{1 - \cos 2x}%1(1)clearclcsyms x;f=(tan(x)-sin(x

plot3 基本的三维曲线图绘制plot3(x,y,z)，x,y,z均为相同长度的向量,会得到三个向量相同下标构成的的三维坐标(xi,yi,zi)(i=1~n)连的曲线plot3(X,Y,Z)，X,Y,Z均为相同大小的矩阵,均为m*n的矩阵时，plot3命令将绘得m条曲线，每条曲线均为X，Y，Z列向量为（x,y,z）坐标值的曲线例 三维圆柱形螺旋线,x,y,z均为相同长度的向量cl

1.把当前窗口分成四个区域，用不同的颜色和线条分别绘制sin(x),cos(x),e^x,log(x)的函数图形，并加入文字标示和网格。clearclcx=-2*pi:0.01:2*pi;y1=sin(x);y2=cos(x);y3=exp(x);t=0:0.01:20;y4=log(t);subplot(221);plot(x,y1,'b-');title('sin(x)')

matlab中subplot(m,n,p)可以在一个figure中画m×n个子图，p可以指定子图的位置，一般情况下p为单个数字，p为向量时可以合并多个子图为一个子图。看下面这个例子clearclcx=-4*pi+eps:0.01:4*pi;y1=sin(x);y2=cos(x);y3=tan(x);figuresubplot(2,2,1);plot(x,y1);title('

1.选择合适的步长绘制下列函数的图形。(1)ln((1-x)/(1+x)),x∈(-1,1)(2)sqrt(cos(x)),x∈[-π/2,π/2](3)sin(1/t),t∈(-1,0)∪(0,1)(4)sin(x)/x,x∈(-0.5,0)∪(0,0.5)clearclcx1=-1+eps:0.01:1;y1=log((1-x1)/(1+x1));x2=-pi/2:

1.绘制下列图像（1）生成 100 个从 0 到 10 之间的随机整数，绘制图形；（2）生成 10个从 0 到 10 之间的随机整数，绘制点图，要求为绿色圆圈；（3）生成0到1之间36个等差数据作为横坐标，画出6阶魔方阵元素图形（将其按行展成一行），要求为黄色上三角，点划线。%(1)clearclcx=round(rand(1,100)*10);plot(x);%(2)

统计函数sum(A)求和如果A是向量,sum(A)计算每个元素之和如果A是矩阵,sum(A)计算每列元素之和，(即把每列看成一个向量，然会返回一个行向量，结果为每列元素之和)如果A是多维数组,sum(A,dim)把指定的维度对应位置求和，即dim=1时对列求和，dim=2时对列求和B = sum(..., 'double')和B = sum(..., dim,'double')

函数绘图分为二维绘图和三维绘图其中最重要的两个函数为plot和plot3，其它函数基本上是对其功能的扩充1.plot(y)如果y是一个向量，画出的二维图，横轴y的下标，纵轴为y的每一个数，如果y是一个复数，那么相当于plot(real(y),imag(y))，即横轴为实部，纵轴为虚部例1，y为实数clearclcx=[0 0.7 0.3 0.5 0.1 0.6 1.2

SCILIB是一款优秀的数值计算开源软件，与Matlab不同的是，Matlab是商业软件，购买的时候按基本组建和工具箱买，就算是有财力的高校基本也就只购买一两套，毕竟一套几十W呢。说SCILIB是数值计算软件没错，对于符号计算方面，SCILIB也可以通过安装工具箱来实现。SCILIB下载地址：http://www.scilab.org/在人教版高中数学教材中，人教B版在必修一附录和