题目
Mathematica问题
1.用牛顿迭代法求函数g(x)=(x - 5)^(-2/3)*(x - 17)^(5/3)的根,并画出其迭代过程的蛛网图.
2.当圆沿着给定的曲线滚动时,试建立相对于该圆位置固定的点的轨迹所对应的参数方程,并通过Mathematica模拟其形成过程和观察最终的曲线形状与特征,其中定点的位置为(1)在圆周外(2)在圆上(3)在圆内
3.随机产生90个正整数,(介于1-365之间),用这90个正整数代表一个班90个同学的生日),然后统计数据,观察是否有3人以上的人生日相同.当90个人中有3个人生日相同时,输出结果是“1”,否则为“0”.如此重复200次,计算这一事件发生的概率?要求模拟三次以上并记录下每次模拟计算的频率值?
要求:
需要mathematica源代码
关于第一个问题,方程明显有一个根x=17
第三个问题指的是有三个同学生日相同即可。
关于第二题,先建立圆上定点轨迹的参数方程,注意是【滚动】,像做摆线一样。把你给的程序修改后得到的是个滑动的结果。
1.用牛顿迭代法求函数g(x)=(x - 5)^(-2/3)*(x - 17)^(5/3)的根,并画出其迭代过程的蛛网图.
2.当圆沿着给定的曲线滚动时,试建立相对于该圆位置固定的点的轨迹所对应的参数方程,并通过Mathematica模拟其形成过程和观察最终的曲线形状与特征,其中定点的位置为(1)在圆周外(2)在圆上(3)在圆内
3.随机产生90个正整数,(介于1-365之间),用这90个正整数代表一个班90个同学的生日),然后统计数据,观察是否有3人以上的人生日相同.当90个人中有3个人生日相同时,输出结果是“1”,否则为“0”.如此重复200次,计算这一事件发生的概率?要求模拟三次以上并记录下每次模拟计算的频率值?
要求:
需要mathematica源代码
关于第一个问题,方程明显有一个根x=17
第三个问题指的是有三个同学生日相同即可。
关于第二题,先建立圆上定点轨迹的参数方程,注意是【滚动】,像做摆线一样。把你给的程序修改后得到的是个滑动的结果。
提问时间:2020-11-25
答案
第一个问题:
f[x_] := 2 x^2;
xlist = NestList[# - f[#]/f'[#] &,0.5,4]; ylist = f[xlist]; list1 = Transpose[{xlist,ylist}]; list2 =
Transpose[{Drop[xlist,1],Table[0,{Length[xlist] - 1}]}];
list = Riffle[list1,list2];
iterplot = ListLinePlot[list]; funplot = Plot[f[x],{x,0,1},PlotRange -> {0,1},PlotStyle -> {Black}]; Show[{funplot,iterplot}]
第二个问题:
singleTest := Module[{data,dataSorted,exist,howmanysame = 2},
data = RandomInteger[{1,365},90]; dataSorted = Sort[data];
If[Or @@ (Map[dataSorted[[#]] == dataSorted[[# + howmanysame]] &,
Range[Length[data] - howmanysame]]),exist = 1,exist = 0];
exist];
result = Table[singleTest,{200}];
In[203]:= Total[result]/200 // N
Out[203]= 0.505
第三个问题:
x[t_] := t; y[t_] := Sin[t]; r = 0.5;
Manipulate[
Show[
{ParametricPlot[{x[t],y[t]},{t,0,2 \[Pi]},
PlotRange -> {{0,2 \[Pi]},{-2,2}}],
Graphics[{
Circle[{x[a],y[a]} + r*{- y'[a],x'[a]}/Norm[{x'[a],y'[a]}],r],
Point[{x[a],y[a]} + r*{- y'[a],x'[a]}/Norm[{x'[a],y'[a]}]]}]}
],
{a,0,2 \[Pi]}]
第二个和第三个粘贴反了.一些意见,1里面方程似乎没有根;3里面到底是三个还是三个以上说的不清楚.这个让人很不爽.
2里面我只做了圆心,算是hint吧,其他可以自己搞定的
那么1,3解决了
关于2你说的无摩擦滚动=滑动+绕圆心一定速率的转动.依我看只有用弧长参数代码才容易写.不过就是x,y各增加一项,琢磨琢磨吧:)
f[x_] := 2 x^2;
xlist = NestList[# - f[#]/f'[#] &,0.5,4]; ylist = f[xlist]; list1 = Transpose[{xlist,ylist}]; list2 =
Transpose[{Drop[xlist,1],Table[0,{Length[xlist] - 1}]}];
list = Riffle[list1,list2];
iterplot = ListLinePlot[list]; funplot = Plot[f[x],{x,0,1},PlotRange -> {0,1},PlotStyle -> {Black}]; Show[{funplot,iterplot}]
第二个问题:
singleTest := Module[{data,dataSorted,exist,howmanysame = 2},
data = RandomInteger[{1,365},90]; dataSorted = Sort[data];
If[Or @@ (Map[dataSorted[[#]] == dataSorted[[# + howmanysame]] &,
Range[Length[data] - howmanysame]]),exist = 1,exist = 0];
exist];
result = Table[singleTest,{200}];
In[203]:= Total[result]/200 // N
Out[203]= 0.505
第三个问题:
x[t_] := t; y[t_] := Sin[t]; r = 0.5;
Manipulate[
Show[
{ParametricPlot[{x[t],y[t]},{t,0,2 \[Pi]},
PlotRange -> {{0,2 \[Pi]},{-2,2}}],
Graphics[{
Circle[{x[a],y[a]} + r*{- y'[a],x'[a]}/Norm[{x'[a],y'[a]}],r],
Point[{x[a],y[a]} + r*{- y'[a],x'[a]}/Norm[{x'[a],y'[a]}]]}]}
],
{a,0,2 \[Pi]}]
第二个和第三个粘贴反了.一些意见,1里面方程似乎没有根;3里面到底是三个还是三个以上说的不清楚.这个让人很不爽.
2里面我只做了圆心,算是hint吧,其他可以自己搞定的
那么1,3解决了
关于2你说的无摩擦滚动=滑动+绕圆心一定速率的转动.依我看只有用弧长参数代码才容易写.不过就是x,y各增加一项,琢磨琢磨吧:)
举一反三
已知函数f(x)=x,g(x)=alnx,a∈R.若曲线y=f(x)与曲线y=g(x)相交,且在交点处有相同的切线,求a的值和该切线方程.
我想写一篇关于奥巴马的演讲的文章,写哪一篇好呢?为什么好
奥巴马演讲不用看稿子.为什么中国领导演讲要看?
想找英语初三上学期的首字母填空练习……
英语翻译
最新试题
- 1一根长2m的硬杆AB水平旋转,A端挂30N的物体,B端挂50N的物体,要使杠杆平衡,支点应距A端多远?(画简图)
- 2360/(18*4)怎么简算
- 3time to go to the park改为同义句
- 4请大家帮我指导一下这个公式是怎么推出的 【物质的量浓度与溶解度S的换算】 S=C.V.M.100/(1000VP-C.V.M
- 5在分式1/f=1/f1+1/f2中,f2≠-f2,则f=
- 6有人知道英文吗,你们知道Main翻译成中问是什么吗
- 7我喜欢倾听优美的舞蹈和动听的音乐,修改病句
- 8长方形的面积公式是怎么来的?为什么把宽乘长=面积?
- 9be worn out from
- 10当a=(),b=()时,下面两个等式成立.a÷b=12余15 a+b=353
热门考点
- 1早餐 午餐 晚餐的英文
- 2新宰的畜禽肉,过一段时间再煮,肉更鲜嫩,原因是由于存在 核糖体?溶酶体?线粒体?中心体?
- 3越南语翻译:anh oi a la ai vay a noi that cho minh hieu vay co duoc ko
- 40.6=( ):10=1/2=( )%=( )÷40=6:( )=( )÷60
- 5形容高兴的受不了的成语、没有明确说明间接透露话语意思的成语?
- 62分之x-3≥3分之2x-5求正整数x
- 7英语翻译
- 84.现欲探究一固体混合物A的成分,已知其中可能含有Mg(OH)2、CuCl2、NH4NO3、 (NH4)2SO4、 CuO五种物质中的两种或多种.按图所示进行实验,出现的现象如图所示(设过程中所有可能
- 9用爱是和可以组一个句子
- 10按照语法结构特点,世界上的语言一般可以分成哪几类