在无风情况下的喷泉模拟 我的python代码 import numpy as np import random import matplotlib matplotlib.rcParams['font.sans-serif']=[u'simHei'] matplotlib.rcParams['axes.unicode_minus']=False import matplotlib.pyplot as plt tt = np.linspace(0,4,1000) #y = -tt**2+4*tt #…

例12:一只游船上有800(1000)人,一名游客不慎患传染病,12(10)小时后有3人发病,由于船上不能及时隔离,问经过60(30)小时,72小时,患此病的人数.(与人口模型和Logistic模型类似) 先用python和matlab模拟 我的python代码 # -*- coding: utf-8 -*- import numpy as np import random import matplotlib matplotlib.rcParams['font.sans-serif']=[u's…

指派授课问题 现有A.B.C.D四门课程,需由甲.乙.丙.丁四人讲授,并且规定: 每人只讲且必须讲1门课:每门课必须且只需1人讲. 四人分别讲每门课的费用示于表中: 课 费用 人 A B C D 甲 2 10 9 7 乙 15 4 14 8 丙 13 14 16 11 丁 4 15 13 9 带包python代码: from scipy.optimize import linear_sum_assignment import numpy as np #cost =np.array([[4,1,3…

1 (1). 随机游走问题.在-10到10的一维线段上,质点以1/5的概率用左运动1单位,以2/5的概率停止不动,以2/5的概率向右运动2单位,且碰到-10时向右运动3单位,碰到10时向左运动4单位.请画出它的轨迹. 我原来的代码matlab: s=0; now=0; for i=1:400 x=unifrnd(0,1); if x<=0.2 if now-1 <= -10 now=now-1+3; else now=now-1; end elseif x>0.6 if now+1 &g…

证明,在椭圆形的音乐厅内,从一个椭圆的一个焦点发出声音,则另一个焦点听到的声音是最大的. 分析:证明,从椭圆的一个焦点任意发射的直线经过反射后,并经过另一个焦点.            画图,过一个焦点随机画一条直线l1,再做它与椭圆的交点,过该点作椭圆的一条切线l2,在作一条与该切线垂直的线l3,再作一条l1关于l3对称的直线l4,然后计算另一个焦点到l4的距离小于一个很小的数即可 我的matlab代码 % 画出一个椭圆 t=linspace(0,2*pi,1000); A=5; B=4; C…

为期三周的数学建模国赛培训昨天正式结束了,还是有一定的收获的,尤其是在MATLAB的使用上. 1. 一些MATLAB的基础性东西: 元胞数组的使用:http://blog.csdn.net/z1137730824/article/details/39206823 对于任意文件夹的同一格式的图片的批量读取:http://blog.csdn.net/haizimin/article/details/39646595 关于MATLAB在微分/偏微分方程以及其他高等数学问题中的应用. 关于MATLAB在…

1.读入 txt 文件数据. load xxx.txt A=load(‘xxx.txt’) A=dlmread(‘xxx.txt’) A=importdata(‘xxx.txt’) 例:将身高体重的 txt 文件数据读入. 2.数据导出到 txt 文件. dlmwrite(‘xxx.txt’,B); 例:生成一个 100 维的随机向量 B,并将其写入 sj.txt 文件. 3.读入 excel 表数据.  num=xlsread(‘xxx.xls’,sheet,Range) [num,txt,r…

一.向量.矩阵的表示和使用 format long  %小数很多format short %默认4位小数format rat %显示最近的分数format short e %指数格式的数 尾数多少 exp(1) %自然对数的底exp(10) log(x) %x的自然对数 底是e log10(x) %以10为底 asin(pi) atan(pi/3) %反三角函数 向量(vector)一维数值数组.MATLAB 允许你创建列向量和行向量,列向量通过在方 括号内把数值用分号(;)隔开来创建,对元素的…

MATLAB学习笔记——函数式文件 引入函数式文件 说明: 函数式文件主要用于解决计算中的参数传递和函数调用的问题. 函数式的标志是它的第一行为function语句. 函数式文件可以有返回值,也可以没有. 函数式文件执行之后,只保留最后结果,不保留任何中间过程,所定义的变量也仅在函数内部起作用,并随着调用的结束而被清除.即除无特殊声明外,其变量均为局部变量.  注意:存储函数式文件时文件名要与主函数名相同. 语法: 说明: 1.函数定义行定义了函数的输出和输入变量的个数和顺序. 2.H1行和帮助…

抢票是并发执行 多个进程可以访问同一个文件 多个进程共享同一文件,我们可以把文件当数据库,用多个进程模拟多个人执行抢票任务 db.txt {"count": 1} 并发运行,效率高,但竞争写同一文件,数据写入错乱,只有一张票,都卖成功给了10个人 #文件db.txt的内容为:{"count":1} #注意一定要用双引号,不然json无法识别 from multiprocessing import Process import time import json cla…