matlab怎么画出dpdt=(207p^0.17-P)2.7的相图

根据问题的已知条件，对于此微分方程，通过变形可以得到dt/dp形式的方程，然后用dsolve函数解答其解析值，最后用plot函数绘出其相位图。求解过程如下：

>> syms t(p)

>> Dt=diff(t);

>> t=dsolve(Dt==27/(207p^017-p),t(0)==0)；%给出的问题缺初始条件

>> p=0:20;

>> t=eval(t);

>> plot(t,p)

>> xlabel('t'), ylabel('p')

运行结果

给你一个全MATLAB仿真的程序，没用到SIMULINK

close all

clear all

a=newfis('fuzzf');

f1=1;

a=addvar(a,'input','e',[-3f1,3f1]);

a=addmf(a,'input',1,'NB','zmf',[-3f1,-1f1]);

a=addmf(a,'input',1,'NM','trimf',[-3f1,-2f1,0]);

a=addmf(a,'input',1,'NS','trimf',[-3f1,-1f1,1f1]);

a=addmf(a,'input',1,'Z','trimf',[-2f1,0,2f1]);

a=addmf(a,'input',1,'PS','trimf',[-1f1,1f1,3f1]);

a=addmf(a,'input',1,'PM','trimf',[0,2f1,3f1]);

a=addmf(a,'input',1,'PB','smf',[1f1,3f1]);

f2=1;

a=addvar(a,'input','ec',[-3f2,3f2]);

a=addmf(a,'input',2,'NB','zmf',[-3f2,-1f2]);

a=addmf(a,'input',2,'NM','trimf',[-3f2,-2f2,0]);

a=addmf(a,'input',2,'NS','trimf',[-3f2,-1f2,1f2]);

a=addmf(a,'input',2,'Z','trimf',[-2f2,0,2f2]);

a=addmf(a,'input',2,'PS','trimf',[-1f2,1f2,3f2]);

a=addmf(a,'input',2,'PM','trimf',[0,2f2,3f2]);

a=addmf(a,'input',2,'PB','smf',[1f2,3f2]);

f3=15;

a=addvar(a,'output','u',[-3f3,3f3]);

a=addmf(a,'output',1,'NB','zmf',[-3f3,-1f3]);

a=addmf(a,'output',1,'NM','trimf',[-3f3,-2f3,0]);

a=addmf(a,'output',1,'NS','trimf',[-3f3,-1f3,1f3]);

a=addmf(a,'output',1,'Z','trimf',[-2f3,0,2f3]);

a=addmf(a,'output',1,'PS','trimf',[-1f3,1f3,3f3]);

a=addmf(a,'output',1,'PM','trimf',[0,2f3,3f3]);

a=addmf(a,'output',1,'PB','smf',[1f3,3f3]);

rulelist=[1 1 1 1 1;

1 2 1 1 1;

1 3 2 1 1;

1 4 2 1 1;

1 5 3 1 1;

1 6 3 1 1;

1 7 4 1 1;

2 1 1 1 1;

2 2 2 1 1;

2 3 2 1 1;

2 4 3 1 1;

2 5 3 1 1;

2 6 4 1 1;

2 7 5 1 1;

3 1 2 1 1;

3 2 2 1 1;

3 3 3 1 1;

3 4 3 1 1;

3 5 4 1 1;

3 6 5 1 1;

3 7 5 1 1;

4 1 2 1 1;

4 2 3 1 1;

4 3 3 1 1;

4 4 4 1 1;

4 5 5 1 1;

4 6 5 1 1;

4 7 6 1 1;

5 1 3 1 1;

5 2 3 1 1;

5 3 4 1 1;

5 4 5 1 1;

5 5 5 1 1;

5 6 6 1 1;

5 7 6 1 1;

6 1 3 1 1;

6 2 4 1 1;

6 3 5 1 1;

6 4 5 1 1;

6 5 6 1 1;

6 6 6 1 1;

6 7 7 1 1;

7 1 4 1 1;

7 2 5 1 1;

7 3 5 1 1;

7 4 6 1 1;

7 5 6 1 1;

7 6 7 1 1;

7 7 7 1 1];

a=addrule(a,rulelist);

a1=setfis(a,'DefuzzMethod','mom');%Defuzzy

writefis(a1,'fuzzf');

a2=readfis('fuzzf');

Ulist=zeros(7,7);

for i=1:7

for j=1:7

e(i)=-4+i;

ec(j)=-4+j;

Ulist(i,j)=evalfis([e(i),ec(j)],a2);

end

end

figure(1);

plotfis(a2);

figure(2);

plotmf(a,'input',1);

figure(3);

plotmf(a,'input',2);

figure(4);

plotmf(a,'output',1);

这里简单说明一下：首先是编写2个输入，1个输出的隶属度函数；接下来的是模糊规则，一共49条；然后用解模糊函数得出控制量U，这里输出的U就直接是精确量了，解模糊用到得规则是取隶属度最大的那个数即MOM算法。

显示的三个图形窗口分别是：模糊控制器内部原理图，以及2个输入，1个输出的隶属度函数图。

求解微分方程需要指定初始条件，这里以 f(0) = 01 为例。

解微分方程：

>> f=dsolve('Df=pf^2','f(0)=01','p') 

f = 

-2/(p^2-20)

画f对p的函数图像：

ezplot(f,[0 10])

画df/dp对p的函数图像：

ezplot(diff(f,sym('p')),[0 10])

1、在电脑桌面上打开MATLAB。

2、当你没有把你编辑好的子函数放入正确的文件夹时，不管你输入的程序是否有错，当你执行时，系统都会提示“没有定义”，这个时候就是没有子函数文件放到主函数库文件夹中。

3、首先，输入“edit conv”这个“conv”是我准备编辑的子函数的主函数名，每个人所需的主函数名都不一样，这里可以不要那么在意。然后就可以打开相应的主函数，这是可以看到主函数所在文件夹位置。如图所示，※记住这个位置。

4、然后进行子函数的编辑，首先新建一个新的m文件，然后在编辑框中输入所要的子函数，例如我就是输入了convnem的函数，这是conv的子函数。然后点击保存。

5、接下来，在“conv”所在文件位置，将位置记录下来，将子函数的文件保存到主函数所在的文件夹即可，然后在重新的将你所需的程序（含有子函数的）输入并且执行，这时你会发现原来提示的错误没有了，那么这时调用自己添加的子函数成功了。以后有相应的子函数都可以通过这样的方法进行添加子函数。

用dsolve('Dp=a+bp','t')可以解得：p=-(a - cexp(bt))/b

然后求在最小二乘原则下的最佳的a、b值。

clear;

clc;

p=[1 4 5 6 10 11 ];

t=[86 86 49 130 114 103];

f=@(var) sum((p+(var(1) - var(3)exp(var(2)t))/var(2))^2);%a=var(1),b=var(2),c=var(3)

nvars=3;

[var,fval,exitflag]=ga(f,nvars)

求得：

a= 32038

b= -05185

matlab中调用refpropm如何计算dp/dv，或者是定温度下的dh/dv？

其计算依据，是diff()函数，其使用格式

syms x

y=x^3

diff(y,x) %dy/dx

结果为3x^2

以上就是关于matlab怎么画出dp/dt=(207p^0.17-P)/2.7的相图全部的内容，包括:matlab怎么画出dp/dt=(207p^0.17-P)/2.7的相图、模糊PID控制器的MATLAB仿真程序、MATLAB微分方程图像的问题 df/dp=p*f^2 求f与p函数图像的算法，以及df/dp与p函数图像的算法 急求，救急!等相关内容解答，如果想了解更多相关内容，可以关注我们，你们的支持是我们更新的动力！