司老师第九讲QQ讨论部分文档格式.docx
- 文档编号:20308480
- 上传时间:2023-01-21
- 格式:DOCX
- 页数:13
- 大小:492.72KB
司老师第九讲QQ讨论部分文档格式.docx
《司老师第九讲QQ讨论部分文档格式.docx》由会员分享,可在线阅读,更多相关《司老师第九讲QQ讨论部分文档格式.docx(13页珍藏版)》请在冰豆网上搜索。
%计算每月的还款差额
xk=A0/N+(N+1)*A0*r/N/2
%计算每月还款额
33:
11
p=[0.08
0.9
0.02;
0.08
0.24
0.7
0.06];
a=[p'
-eye(3);
ones(1,3)];
%构造非线性方程组的系数阵
b=[zeros(3,1);
1];
%构造非线性方程组的常数项列
plimit=a\b
37:
02
[vec,val]=eigs(p'
1);
vec=vec/sum(vec);
[plimit,vec]
38:
44
[1]许哲,李号雷,基于灰色模型和Bootstrap理论的大规模定制质量控制方法研究,数学的实践与认识,第42卷第21期,121~127,2012年11月。
47:
27
cumsum(x)
33
就是累加运算
海航司守奎(94339146)20:
00:
08
x0=[2.5320,2.6470,2.6290,2.5840,2.6090,2.6010,2.5280,2.5630,2.6540,2.6190];
n=length(x0);
me=quantile(x0,0.5)
%计算中位数
[h,p,stat]=runstest(x0,me)
%进行游程检验
x1=cumsum(x0);
%求累加序列
zk=(x1(1:
end-1)+x1(2:
end))/2
%求累加序列的均值序列
B=[-zk'
ones(size(zk'
))];
yn=x0(2:
end)'
;
ab=B\yn
%拟合参数a,b
x=dsolve('
Dx+a*x=b'
'
x(0)=x00'
)
%求微分方程的符号解
x=subs(x,{'
a'
b'
x00'
},{ab
(1),ab
(2),x0
(1)});
xx=vpa(x,6)
yuce=subs(x,'
t'
[0:
n+5]);
%求累加序列的预测值
yuce0=[x0
(1),diff(yuce)]
%求原始数据的预测值
c=std(yuce0(1:
n))/std(x0)
%求后验差比值c
nyuce=yuce0(n+1:
end)
%提取6个新的预测值
nyb=[x0,
nyuce];
%构造新的样本数据
nnyb=reshape(nyb,[4,4])
mu=mean(nnyb)
%分别求4个子样本的均值
jc=range(nnyb)
%分别求4个子样本的极差
xlswrite('
hb.xls'
[nnyb;
mu;
jc])
%把数据写到Excel文件中,便于做表使用
b=rand(16,250);
%产生16行250列的随机数矩阵
h=floor(n*b)+1;
%把每一列的随机数映射为相应的编号
mx=nyb;
bb=mx(h);
bb=reshape(bb,[4,1000]);
mmu=mean(bb);
mjc=range(bb);
%计算1000个子样本的均值和极差
smu=sort(mmu);
sjc=sort(mjc);
%把均值和极差按照从小到大的次序排列
alpha=0.01;
k1=floor(1000*alpha/2),
k2=floor(1000*(1-alpha/2))
mqj=[smu(k1),
smu(k2)]
%显示均值的置信区间
jqj=[sjc(k1),
sjc(k2)]
%显示极差的置信区间
subplot(1,2,1),
plot(mu,'
*-'
),
hold
on
plot([1,4],[mqj
(1),mqj
(1)]),
plot([1,4],[mqj
(2),mqj
(2)])
ylabel('
样本均值'
)
subplot(1,2,2),
plot(jc,'
plot([1,4],[jqj
(1),jqj
(1)]),
plot([1,4],[jqj
(2),jqj
(2)])
极差'
yuce0=double(yuce0);
14:
03
82.99217.1948.4
88.09321.3969.6
99.99625.19710.4
105.39429.09711.4
117.710034.010012.2
131.010140.010114.2
148.210544.010415.8
161.811249.010917.9
174.211251.011119.6
184.711253.011120.8
185.310848.610318.9
175.999529920.5
166.810648.810219.2
156.810549.798.818.6
148.99350.696.517.9
136.59647.997.917.6
49
d=load('
huigui2.txt'
);
xs=[ones(16,1),d(:
[1:
4])];
cs=xs\d(:
5)
[cs2,cs2int,r,rint,stats]=regress(d(:
5),xs)
[cs,cs2]
36:
35
[1,3])];
[cs1,cs1int,r,rint,stats]=regress(d(:
40:
04
a=d(:
4]);
b=d(:
5);
rstool(a,b)
45:
30
a=load('
huigui.txt'
AA=a(:
end-1]);
B=a(:
end);
rstool(AA,B)
n=size(a,1);
xs=[ones(n,1),AA,AA(:
1).*AA(:
2),AA(:
3),AA(:
4),AA(:
2).*AA(:
3).*AA(:
4),AA.^2]
cs=xs\B
%为了验证完全二次式的回归答案
57:
huigui3.txt'
x0=a(:
y0=a(:
end)
stepwise(x0,y0,[1:
8])
%初始时使用全部8个自变量
海航司守奎(94339146)21:
02:
28
stepwise(x0,y0)
%初始时不使用任何变量
04:
下面开始答疑
201304001(903528428)21:
05:
46
马尔科预测的适用范围是什么?
微分霸天虎(1355097263)21:
07:
“匿名函数”格式是什么?
预测极限状态
34
函数格式:
function
返回值=函数名(输入参数)
10
匿名函数格式:
返回值=@(输入参数)
函数表达式
201307031(820014430)21:
10:
32
逐步回归
向前选择,向后删除,选择哪个方法更好?
56
何时用“匿名函数”?
岭回归的含义是什么?
1]具体是怎样的内容?
16:
25
18:
如果有多个自变量,但自变量间可能相关,这用主成分回归或者偏最小二乘回归比较好,对吗?
22:
下次课讲?
47
建议老师尽量多分析程序来源
48
最后一次课讲随机模拟
23:
00
什么意思?
18
讲讲预测方法各自的特点和适用范围
24:
19
就是matlab程序的构建
21
谢谢老师!
27:
40
28:
12
Description
Y
=
wblpdf(X,A,B)
computesthe
Weibull
at
each
of
the
values
in
X
usingthe
corresponding
scale
parameter,
A
and
shapeparameter,
B.
X,
A,and
B
can
be
vectors,
matrices,
or
multidimensionalarrays
that
all
have
same
size.
scalar
input
is
expanded
toa
constant
array
size
as
other
inputs.
The
parametersin
must
positive.
is
Some
references
refer
to
distribution
with
a
singleparameter.
This
corresponds
wblpdf
1.
Examples
37
29:
画图的Matlab程序如下
%数学与Matlab参数:
b=tau;
a=beta^(-1/b);
subplot(121)
fx=@(x,tau,beta)tau*beta*x.^(tau-1).*exp(-beta*x.^tau);
fplot(@(x)fx(x,2,0.5),[0,5],'
.'
text(1.6,0.5,'
$\beta=0.5,
\tau=2$'
Interpreter'
Latex'
xlabel('
$x$'
$f(x)$'
title('
数学上的参数表示'
subplot(122)
fx2=@(x,a,b)wblpdf(x,a,b);
fplot(@(x)fx2(x,0.5^(-1/2),2),[0,5],'
str=['
a='
num2str(0.5^(-1/2)),'
b=2'
];
text(1.6,0.5,str)
Matlab的参数表示'
30:
201307014--程小乐(646092267)09:
这是是不是不能在matlab
7.0上运行啊
201307019(347667369)09:
19:
为啥
呃,那倒真有可能
17
工具箱里没有。
。
201307004(82985573)09:
49:
07
用2013a吧
我们机房安装的是7.0的
内民大-刘京(1374413450)16:
25:
net=newff(minmax(P),[8,5],{'
logsig'
},'
traingdx'
这个BP网络里面的
minmax(P)是什么意思?
201307019(347667369)16:
32:
45
minmax不是最小值和最大值的意思么。
35:
xm=minmax(x)%求x的最小值和最大值
15
这步变化是怎么得到的?
p是个啥?
P是我定义的矩阵
矩阵?
对
59
呃,对矩阵做minmax我就不清楚了,要不help一下看看
201307029(785144557)16:
54
是不是第一行最小值0,最大值1?
第二行也是?
41:
我DOC了一下。
这个例子。
完全看不明白。
26
按行返回结果?
doc
哪个
按列返回结果。
42
DOC
minmax
42:
那就不知道了
哦
bp是什么
201307005(1445891025)16:
39
minmax(A)即是求解A矩阵中每个行向量中最小值和最大值
58
bp里面用的不是mapminmax吗
内民大-刘京(1374413450)17:
16
哦哦哦哦哦。
我去。
23
这么简单
就是没想明白。
谢谢
那mapminmax是啥意思?
201307017(122192130)17:
标准化?
201307005(1445891025)18:
mapminmax是归一化处理,防止bp网络里面如果有些参量过大可能影响bp网络的结果。
.
31:
53
88
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 老师 第九 QQ 讨论 部分