matlab
daijiangtao
这个作者很懒,什么都没留下…
展开
-
MATLAB中floor、round、ceil、fix区别
Matlab取整函数有: fix, floor, ceil, round.具体应用方法如下:fix朝零方向取整,如fix(-1.3)=-1; fix(1.3)=1;floor,顾名思义,就是地板,所以是取比它小的整数,即朝负无穷方向取整,如floor(-1.3)=-2; floor(1.3)=1;floor(-1.8)=-2,floor(1.8)=1ceil,与floor相反,它的意思转载 2015-07-03 10:47:44 · 620 阅读 · 0 评论 -
matlab中find函数的使用
语法:1. ind = find(X)2. ind = find(X, k)3. ind = find(X, k, 'first')4. ind = find(X, k, 'last')5. [row,col] = find(X, ...)6. [row,col,v] = find(X, ...)说明:1. ind = find(X) 找出矩阵转载 2015-07-26 14:35:30 · 961 阅读 · 0 评论 -
ode求解器的事件(Event)属性
检测事件matlab微分方程如何设定变量的范围 如dy1=y2dy2=y1+1其中y1的范围为0【解】M文件:function [value,isterminal,direction] = events1(t,y)value = y(1)-4;isterminal= 1;direction = 0;命令窗口:dy = @(t,转载 2015-07-07 09:58:18 · 13119 阅读 · 3 评论 -
Matlab在曲线中如何标示极值点?
>> x =1:0.1:10;>> y=sin(x)./x;>> plot(x,y)就是把它的最大值最小值用虚线标示出来也可以。解答1: x =1:0.1:10;y=sin(x)./x;plot(x,y)[a,b]=findpeaks(y);hold onplot(x(b),a,'ro')[c,d]=findpeaks(-y);hold on转载 2015-07-07 10:15:30 · 30134 阅读 · 2 评论 -
时滞微分方程的Matlab解法之一dde23
ddex1histz = @(t) 2*ones(2,1);ddex1dez = @(t,y,Z) [y(1)*(1 + 0.1*sin(t)-0.1*Z(1,1) - y(2)/(1+y(1)) ); y(2)*( (2+sin(t))*10^(-5) + 9*Z(1,2)/(1+Z(1,2)) - Z(2,1) )]; sol = dde23(ddex1dez转载 2015-07-03 16:12:43 · 38409 阅读 · 5 评论 -
时滞微分方程的matlab解法之二dde23---另外一些例子
1. function testclear;clc;lags=[1]; tspan=[0,100];sol=dde23(@ddefun,lags,@history,tspan); plot(sol.x,sol.y); title('ddefun');xlabel('t');ylabel('y'); legend('y_1','y_2','y_3','y_4'原创 2016-08-19 15:20:05 · 17127 阅读 · 1 评论 -
时滞微分方程求解之三ddesd--变时滞
考虑下面的变时滞微分方程试绘制上述含有时变时滞微分方程的图像分析:该方程应该在某个t0时间之后成立,初始值必须是定义在t0之前的一个关于t的单值向量函数phi(t)。我假设t0=0吧,phi(t)=[1;-1];matlab程序:function ddeext0 = 0;tfinal = 5;tspan = [t0, tfinal];sol = dd转载 2016-08-20 12:27:51 · 15165 阅读 · 6 评论