[MATLAB/编程]报童的诀窍/报童问题-图解法和二分法
MATLAB编程题目 : 报童的诀窍
关于每天报纸购进量的优化模型:
已知b为每份报纸的购进价,a为零售价,c为退回价(a>b>c),每天报纸的需求量为r份的概率是f( r )(r=0,1,2,…)。求每天购进量n份,使日平均收入,即
达到最大。
视r为连续变量,f( r )转化为概率密度函数p( r ),则所求n*满足
实验要求:
已知b=0.75,a=1,c=0.6,r服从均值=500(份),均方差=50(份)的正态分布。报童每天应购进多少份报纸才能使平均收入最高,这个最高收入是多少?
1.的图解法(编程)
参考实验12-1的实验要求2,编写程序,求出n*。
初始根区间为[450,550],增量取0.01,随着根区间变小可改为0.001。
clc;
clear;syms y1 y2 a b c z
a = 1;
b = 0.75;
c = 0.6;
z = 450:0.01:550;y1 = (a-b)/(a-c);
y2 = normcdf(z,500,50)-normcdf(0,500,50); %分布函数plot([450,550],[y1,y1])
hold on;
plot(z,y2);grid on;% 输出交点
fsolve(@(z) (a-b)/(a-c) - normcdf(z,500,50)-normcdf(0,500,50),[450 550])
下图可见,交点在520左侧
补充:
关于匿名函数
@运算符创建句柄,@(x)为函数的入口,指定函数当中的变量x为函数的输入,在接受函数和进行求解时,对象为x。
如下例子,
sqr = @(z) 2.*z
sqr(3)% 输入的变量默认输入到z当中
-> 6
关于fsolve
fsolve(@(x) f(x))
解决的问题为f(x)=0
如下例子
fsolve(@(z) (a-b)/(a-c) - normcdf(z,500,50)-normcdf(0,500,50),[450 550])
解释 : 在z = [450,550]之间,寻找(a-b)/(a-c) - normcdf(z,500,50)-normcdf(0,500,50) = 0
的解
2.的数值解法(二分法/编程)
参考实验12-1的实验要求3,编写程序,求出n和G(n)。
初始根区间为[450,550]。
erfen.m
function y=erfen(a,b) %二分法
if fun(a)*fun(b)<0 c=(a+b)/2;while abs(fun(c))>1.0e-4 % 二分分界点if fun(a)*fun(c)<0b=c; c=(a+b)/2;elseif fun(c)*fun(b)<0a=c; c=(a+b)/2;endendy=c;
elseif fun(a)==0y=a;
elseif fun(b)==0y=b;
elsedisp('区间中不存在根');
end
return;
main.m
function main()
clc;
n = erfen(450,550) a=1;b=0.75;c=0.6;
r=n+1;
while (a-b)*n*normpdf(r,500,50)>exp(-6) r=r+1;
end
r=n+1:r;
G = sum((a-b)*n*normpdf(r,500,50)); % 离散值求和
r=0:n;
G = G+sum(((a-b)*r-(b-c)*(n-r)).*normpdf(r,500,50)) % 离散值求和
return;
fun.m
function y=fun(z) %方程
a = 1;
b = 0.75;
c = 0.6;
y = (a-b)/(a-c) - normcdf(z,500,50) - normcdf(0,500,50); % 构造函数差,方便二分查找
return;
1.计算正态变量的概率密度函数的调用形式为:Y=normpdf(X,mu,sigma)
正态变量的概率密度函数为
其中:X是x的一组值,Y对应一组函数值。
mu为μ,sigma为σ。
当μ=0,σ=1时,为标准正态变量的概率密度函数。
2.计算正态变量的分布函数的调用形式为:P=normcdf(X,mu,sigma)
正态变量的分布函数为且
标准正态变量的概率密度函数对应标准正态变量的分布函数。
[MATLAB/编程]报童的诀窍/报童问题-图解法和二分法相关推荐
- 概率模型——轧钢中的浪费、报童的诀窍
实验13-1 轧钢中的浪费 1.(不知道要不要画图) l=2;sigma=0.2; a=-3; b=2.5; deta=0.5; z=a:deta:b; fz=(1-normcdf(z,0,1))./ ...
- Matlab编程与数据类型 -- 单元数组
Matlab编程与数据类型 – 单元数组
- matlab 天线设计 泰勒加权_微带天线设计尺寸MATLAB编程及其仿真验证
龙源期刊网 http://www.qikan.com.cn 微带天线设计尺寸 MATLAB 编程及其仿真 验证 作者:杨小敏 母玉泽 严月 郭小康 马波 张栋 莫骄弟 来源:<中国科技博览> ...
- MATLAB编程与应用系列-第2章 数组及矩阵的创建及操作(4)
本系列教程来源于出版设计<基于MATLAB编程基础与典型应用书籍>,如涉及版权问题,请联系:156204968@qq.com. 出版社:人民邮电出版社, 页数:525. 本系列教程目前基于 ...
- matlab编程小结
最近在用matlab处理一些开源数据,需要用到matlab编程. 之前学习matlab编程不够系统和全面,对很多函数不太了解.数字图像处理和矩阵的一些常用处理居然不太熟悉,停留在学习的舒适区,有点不应 ...
- matlab编程风格
matlab编程风格 参考博客:https://blog.csdn.net/yuandm819/article/details/52453674?utm_medium=distribute.pc_re ...
- matlab编程选择语句,matlab编程控制语句
matlab编程控制语句 M文件 程序控制结构 函数文件 程序举例 程序调试,2.1 M文件 2.1.1 M文件概述 用MATLAB语言编写的程序,称为M文件.M文件可以根据调用方式的不同分为两类:命 ...
- Matlab:利用Matlab编程实现模拟分子布朗运动的动画展示
Matlab:利用Matlab编程实现模拟分子布朗运动的动画展示 目录 输出结果 实现代码 输出结果 实现代码 %Brownian motion clf; n=20; s=0.02; x = rand ...
- Algorithm之PrA:PrA之nLP非线性规划算法经典案例剖析+Matlab编程实现
Algorithm之PrA:PrA之nLP整数规划算法经典案例剖析+Matlab编程实现 目录 有约束非线性规划案例分析 1.投资决策问题 2.利用Matlab实现求解下列非线性规划 无约束极值问题 ...
- 互信息的数学解释以及matlab编程
1.互信息 互信息是信息论中的一个基本概念,通常用于描述两个系统间的统计相关性,或者是一个系统中所包含另一个系统中信息的多少. 引入互信息对图像进行配准是因为:不需要假定两幅图像灰度的对应关系,不需要 ...
最新文章
- UI设计培训技术分享:搞定萌萌哒可爱图标
- Java编程思想读书笔记--第21章并发
- pip包重新安装mysql_Ubuntu pip 安装 mysql-python包出错
- 制冷设备维保常见的10大故障!
- DeepFakes天敌来了!伯克利紧急研发“火眼金睛”防伪克星
- SAP 电商云 Spartacus UI added-to-cart 的端到端测试源代码解析
- 阿里技术|Java 的这些坑,你踩到了吗?
- python3.x教程_Python3.x基础教程1
- 关于 Java 9 你所需要知道的一切
- qdir 类似工具_实用工具——多窗口资源管理器qdir
- [java编程题]打印指定年指定月份的日历
- THINKPHP官方推荐第三方插件库
- 计算机辅助翻译与人工翻译论文,浅析计算机辅助翻译与纯人工翻译的利与弊
- 【GNSS】GREAT多频多系统GREAT-UPD开源代码-第一章 介绍
- Sublime Text2 好用技巧备忘
- python的if __name__ == “__main__“语法错误SyntaxError: invalid syntax
- AndroidStudio:Image not loaded Try to open it externally to fix format problem
- (附源码)Node.js图书管理小程序 毕业设计 250858
- HashMap初始化长度设置大小
- 关于录取志愿者培训人员名单的通知
热门文章
- 用Html制作网页简单例子
- 【算法java版09】:利用java实现对二进制数进行AMI编码
- matlab绘直方图,[转载]matlab中函数bar绘制直方图
- 问卷与量表数据分析(SPSS+AMOS)学习笔记(一) :问卷分析的流程
- 大众新能源电动车ID.3ID.4技术培训教材手册维修手册电路图
- 什么是网络处理器?网络处理器有哪些应用场景?
- 32位微型计算机原理接口,32位微机原理及接口技术
- “无线射频识别”可节省1/3物流成本
- 谷歌地球 Google Earth v7.3.2.5495 专业版
- [汇编] 最简单的汇编程序