Диссертация (1137458), страница 22

Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

Текст из файла (страница 22)

Великая депрессия 1929-1933гг. (падение индекса Доу-Джонса на90% за 2,5 года) – мировой экономический кризис,1703. «Черный понедельник» 19.10.1987г в США (падение индекса ДоуДжонса на 22,6% за 1 день), вызванный распространением автоматических программ для трейдинга,4. Dotcom кризис в 2000г. в США (индекс высокотехнологичныхкомпаний Nasdaq упал на 82% за 2000-2002гг.) – спекулятивныйпузырь цен, надутый большими ожиданиями относительно развития интернет-бизнеса,5. Великая рецессия 2008-2012гг.

(падение индекса S&P500 на 52%за 1,5 года) – мировой экономический кризис, начавшийся с кризиса в финансовом секторе США.Во время биржевых кризисов большая часть трейдеров теряетденьги. Заработать на биржевом кризисе, конечно же, можно – дляэтого нужно продать акции на пике, непосредственно накануне обвалацен, а также встать по возможности в короткую позицию, т.е. продатьзаемные акции и вернуть их уже по сильно сниженной цене посленаступления кризиса.

При знании момента начала биржевого крахаприбыль от использования такой стратегии огромна, однако для ееиспользования необходимо точно знать момент краха, т.е. быть инсайдером. Подавляющее большинство трейдеров не имеют доступа ктакой информации и во время кризиса теряет значительные средства.Приложение 2. Комплекс программ для исследованиявозможностей участников фондовой биржи по работес финансовыми инструментамиИсходный код для базовой модели из подраздела 4.2.1library(fBasics)rm(list=ls(all=TRUE))# загружаем данные по ценамData <- read.table("C:/temp/data_n.txt", header=T)# выставление границ вероятностей р в этой серии экспериментов171n1<-0.47n2<-0.57# число экспериментов в серииnn<-100# переменные для критериев успешности стратегийav_wealth<-matrix(0,nrow=(n2-n1)*100+1,ncol=nn)better_wealth<-matrix(0,nrow=(n2-n1)*100+1,ncol=nn)bankrupts<-matrix(0,nrow=(n2-n1)*100+1, ncol=nn)# число тактов и число агентов в экспериментахtt<-2500n<-100# величина кредитного рычага (если маржинальна торговля запрещена, то# leverage=0)leverage<-2# достаем данныеprice<- matrix(0,nrow = tt, ncol=1)R<- matrix(0,nrow = tt, ncol=1)date<- Data[F:(tt),1 ]price<- Data[F:(tt),2 ]R<- Data[F:(tt),3 ]month<- Data[F:(tt),4 ]#начальное благосостояниеcash0<-10000stocks0<-0wealth0<-cash0+stocks0*price[1]#матрицы благосостояний агентовcash<-mat.or.vec(n,2)stocks<-mat.or.vec(n,2)wealth<-mat.or.vec(n,tt)#матрицы решений агентовs<-matrix(0, nrow = n, ncol=1)vol<-matrix(0, nrow = n, ncol=1)for (ii in 1:nn) {i1<-0for (pr in n1*100:n2*100) {i1<-i1+1p_min<-n1p_max<-n2rm(wealth)wealth<-mat.or.vec(n,tt)cash[,1]<-cash0stocks [,1]<-stocks0wealth[,1]<-wealth0172#матрица банкротствb<-matrix(0,nrow=1,ncol=tt)# матрица со случайными числами для розыгрыша вероятностиu1 <- runif(n*tt)z1<- matrix(u1, nrow = n, ncol=tt, byrow=TRUE)u2 <- runif(n*tt)z2 <- matrix(u2, nrow = n, ncol=tt, byrow=TRUE)# вектор вероятностей верного определения направления изменения цен исписок, в который записываются значения вероятности у тех, кто обанкротилсяp<-runif(n, min=p_min, max=p_max)m<-c()k<-1# счетчик# Блок 1.

Принятие решений агентамиfor (t in 1:(tt-1)){for (i in 1:n){# те, кто обанкротился, не участвуют в дальнейшей торговлеif (p[i]==0) {wealth[i,t]<-0;stocks[i,1]<-0;cash[i,1]<-0} else {if ((z1[i,t]<p[i])& (price[t+1]-price[t]>0)) { s[i,1] <- 1}else if ((z1[i,t]<p[i])& (price[t+1]-price[t]<=0)) { s[i,1] <- -1}else if ((z1[i,t]>=p[i])& (price[t+1]-price[t]>0)) { s[i,1] <- -1}else if ((z1[i,t]>=p[i])& (price[t+1]-price[t]<=0)) {s[i,1] <- 1}# выставление объема заявкиif ((s[i,1]==-1) & (stocks[i,1]>0)& (R[t]==0)) {if (leverage==0) vol[i,1] <- z2[i,t]*stocks[i,1] else vol[i,1] <z2[i,t]*leverage*wealth[i,t]/price[t]} else if ((s[i,1]==-1) & (stocks[i,1]>0)& (R[t]==1)) {vol[i,1] <- z2[i,t]*stocks[i,1]} else if ((s[i,1]==-1) & (stocks[i,1]<=0) & (R[t]==0)) {vol[i,1] <- z2[i,t]*leverage*wealth[i,t]/price[t]} else if ((s[i,1]==-1) & (stocks[i,1]<=0) & (R[t]==1)) {vol[i,1] <- 0} else if ((s[i,1]==1)&(R[t]==1)) {if (cash[i,1]>=0) vol[i,1] <- z2[i,t]*cash[i,1]/price[t] elsevol[i,1]<-0} else if ((s[i,1]==1)& (cash[i,1]>=0) &(R[t]==0)) {if (leverage==0) vol[i,1]<-z2[i,t]*cash[i,1]/price[t] else vol[i,1]<z2[i,t]*(wealth[i,t]*leverage)/price[t]}else if ((s[i,1]==1)& (cash[i,1]<=0) &(R[t]==0)) {if (leverage==0) vol[i,1]<-0 else vol[i,1]<z2[i,t]*(wealth[i,t]*leverage)/price[t]}173# пересчет параметров после сделокif ((s[i,1]>=0) & (vol[i,1]*price[t]<=cash[i,1])) {cash[i,2]<-cash[i,1]-vol[i,1]*price[t]; stocks[i,2]<stocks[i,1]+vol[i,1]}else if ((s[i,1]>=0) & (vol[i,1]*price[t]>cash[i,1])){cash[i,2] <- 0; stocks[i,2] <- stocks[i,1]+vol[i,1](vol[i,1]*price[t]-cash[i,1])/price[t+1]} else if ((s[i,1]<0)& (stocks[i,1]>0)& (stocks[i,1]>=vol[i,1])) {cash[i,2] <- cash[i,1]+vol[i,1]*price[t]; stocks[i,2] <- stocks[i,1]vol[i,1]} else if ((s[i,1]<0)& (stocks[i,1]>0)& (stocks[i,1]<vol[i,1])) {cash[i,2] <- cash[i,1]+stocks[i,1]*price[t]+ (vol[i,1]stocks[i,1])*(price[t]-price[t+1]); stocks[i,2] <- 0} else if ((s[i,1]<0)& (stocks[i,1]<=0)) {cash[i,2] <- cash[i,1]+vol[i,1]*(price[t]-price[t+1]); stocks[i,2]<- 0}#пересчет благосостояния и анализ банкротствwealth[i,t+1]<-cash[i,2]+stocks[i,2]*price[t+1]cash[i,1]<-cash[i,2]stocks[i,1]<-stocks[i,2]if (wealth[i,t+1]<=wealth0/2) {b[t]<-b[t]+1m[k]<-p[i]p[i]<-0k<-k+1}}}}num<-0for (i in 1:n) if (wealth[i,tt]>=wealth0) num<-num+1better_wealth[i1,ii]<-numbankrupts[i1,ii]<-sum(b)if (bankrupts[i1,ii]<n)av_wealth[i1,ii]<-sum(wealth[1:(n),tt])/(n-bankrupts[i1,ii])else av_wealth[i1,ii]<-0}}av_wealthbetter_wealthbankruptsv<-matrix(10, nrow = n2-n1+1, ncol=1)h<-0; for (j in n1:n2) {h<-h+1; v[h]<-j/100}x11()par(mfrow = c(1, 3))plot(v,rowMeans(av_wealth), xlab="time", ylab="money", main="averagewealth", col="black", type = "l")174plot(v, rowMeans(better_wealth), xlab="time", ylab="percentage",main="fraction of agent with wealth greater than initial",col="black", type = "l")plot(v, rowMeans(bankrupts), xlab="time", ylab="number of agents",main="average number of bankrupts", col="black", type = "l")Исходный код для модели «Лидеры и последователи» из подраздела 4.2.2library(fBasics)rm(list=ls(all=TRUE))# загружаем данные по ценамData <- read.table("C:/temp/data_n.txt", header=T)# выставление границ вероятностей р в этой серии экспериментовn1<-0.4n2<-0.6# число экспериментов в серииnn<-100# переменные для критериев успешности стратегийav_wealth<-matrix(0,nrow=(n2-n1)*100+1,ncol=nn)better_wealth<-matrix(0,nrow=(n2-n1)*100+1,ncol=nn)bankrupts<-matrix(0,nrow=(n2-n1)*100+1, ncol=nn)# число тактов и число агентов в экспериментахtt<-2500n<-100# величина кредитного рычага (если маржинальна торговля запрещена, то# leverage=0)leverage<-2# достаем данныеdate<- matrix(0,nrow = tt, ncol=1)price<- matrix(0,nrow = tt, ncol=1)R<- matrix(0,nrow = tt, ncol=1)date<- Data[F:(tt),1 ]price<- Data[F:(tt),2 ]R<- Data[F:(tt),3 ]month<- Data[F:(tt),4 ]#начальное благосостояниеcash0<-10000stocks0<-0wealth0<-cash0+stocks0*price[1]#матрицы благосостояний агентовcash<-mat.or.vec(n,2)stocks<-mat.or.vec(n,2)wealth<-mat.or.vec(n,tt)#матрицы решений агентовs<-matrix(0, nrow = n, ncol=1)vol<-matrix(0, nrow = n, ncol=1)175ii<-0for (ii in 1:nn) {print(ii)i1<-0for (pr in n1:n2) {i1<-i1+1p_min<-n1p_max<-n2s<-matrix(0, nrow = n, ncol=tt)vol<- matrix(0, nrow = n, ncol=1)cash<-mat.or.vec(n,2)stocks<-mat.or.vec(n,2)wealth<-mat.or.vec(n,tt)cash[,1]<-cash0stocks [,1]<- stocks0wealth[,1]<-wealth0# матрица со случайными числами для розыгрыша вероятностиu1 <- runif(n*tt)z1 <- matrix(u1, nrow = n, ncol=tt, byrow=TRUE)u2 <- runif(n*tt)z2 <- matrix(u2, nrow = n, ncol=tt, byrow=TRUE)p<-runif(n, min=p_min, max=p_max)m2<-c(0)numbers2<-c()k2<-1# счетчик# Блок 1.

Принятие решений агентамиfor (t in 1:(tt-1)){for (i in 1:n){# те, кто обанкротилися не участвуют в дальнейшей торговлеif (p[i]==0) {wealth[i,t]<-0;stocks[i,1]<-0;cash[i,1]<-0} else {# принятие решений лидерамиif (i<=(n/2)) {if ((z1[i,t]<p[i]) & (price[t+1]- price[t]>0)) { s[i,t] <- 1}else if ((z1[i,t]<p[i]) & (price[t+1]-price[t]<=0)) { s[i,t] <-1}else if ((z1[i,t]>=p[i]) & (price[t+1]-price[t]>0)) { s[i,t] <-1}else if ((z1[i,t]>=p[i]) & (price[t+1]-price[t]<=0)) {s[i,t] <1}# принятие решений последователями} else if (i>(n/2)) {if (t==1) s[i,t]<-0else s[i,t]<-s[i-n/2,t-1]176}}}# выставление объема заявки последователямиfor (i in (n/2+1):n){if ((s[i,t]==-1) & (stocks[i,1]>0) & (R[t]==0)) {if (leverage==0) vol[i,1] <- z2[i,t]*stocks[i,1] else vol[i,1] <z2[i,t]*leverage*wealth[i,t]/price[t]} else if ((s[i,t]==-1) & (stocks[i,1]>0) & (R[t]==1)) {vol[i,1] <- z2[i,t]*stocks[i,1]} else if ((s[i,t]==-1) & (stocks[i,1]<=0) & (R[t]==0)) {vol[i,1] <- z2[i,t]*leverage*wealth[i,t]/price[t]} else if ((s[i,t]==-1) & (stocks[i,1]<=0) & (R[t]==1)) {vol[i,1] <- 0} else if ((s[i,t]==1) & (R[t]==1)) {if (cash[i,1]>=0) vol[i,1] <- z2[i,t]*cash[i,1]/price[t] elsevol[i,1]<-0} else if ((s[i,t]==1) & (cash[i,1]>=0) & (R[t]==0) ) {if (leverage==0) vol[i,1]<-z2[i,t]*cash[i,1]/price[t] else vol[i,1]<z2[i,t]*(wealth[i,t]*leverage)/price[t]}else if ((s[i,t]==1) & (cash[i,1]<=0) & (R[t]==0) ) {if (leverage==0) vol[i,1]<-0 else vol[i,1]<z2[i,t]*(wealth[i,t]*leverage)/price[t]}# пересчет параметров после сделокif ((s[i,t]>=0) & (vol[i,1]*price[t]<=cash[i,1])) {cash[i,2]<-cash[i,1]-vol[i,1]*price[t]; stocks[i,2]<stocks[i,1]+vol[i,1]} else if ((s[i,t]>=0) & (vol[i,1]*price[t]>cash[i,1])) {cash[i,2] <- 0; stocks[i,2] <- stocks[i,1]+vol[i,1](vol[i,1]*price[t]-cash[i,1])/price[t+1]} else if ((s[i,t]<0) & (stocks[i,1]>0)& (stocks[i,1]>=vol[i,1])) {cash[i,2] <- cash[i,1]+vol[i,1]*price[t]; stocks[i,2] <- stocks[i,1]vol[i,1]} else if ((s[i,t]<0) & (stocks[i,1]>0)& (stocks[i,1]<vol[i,1])) {cash[i,2] <- cash[i,1]+stocks[i,1]*price[t] + (vol[i,1]stocks[i,1])*(price[t]-price[t+1]); stocks[i,2] <- 0} else if ((s[i,t]<0) & (stocks[i,1]<=0)) {cash[i,2] <- cash[i,1]+vol[i,1]*(price[t]-price[t+1]); stocks[i,2]<- 0}# пересчет благосостояния и анализ банкротстваwealth[i,t+1]<-cash[i,2]+stocks[i,2]*price[t+1]cash[i,1]<-cash[i,2]stocks[i,1]<-stocks[i,2]if ((wealth[i,t+1]<=wealth0/2)&(p[i]!=0)){numbers2[k2]<-ik2<-k2+1p[i]<-0}}}177num1<-0for (i in (n/2+1):n) {if (wealth[i,tt]>=wealth0)num1<-num1+1}better_wealth[i1,ii]<-num1bankrupts[i1,ii]<-k2-1num2<-0for (i in 1:n) {if (wealth[i,tt]>0) {av_wealth[i1,ii]<-av_wealth[i1,ii]+wealth[i,tt]num2<-num2+1}}av_wealth[i1,ii]<-av_wealth[i1,ii]/num2}}# создание графиков критериевv<-matrix(10, nrow = n2-n1+1, ncol=1)h<-0; for (j in n1:n2) {h<-h+1; v[h]<-j/100}x11()par(mfrow = c(1, 3))plot(v,rowMeans(av_wealth,na.rm=TRUE), xlab="time", ylab="money",main="average wealth", col="black", type = "l",ylim=c(min(av_wealth,na.rm=TRUE),max(av_wealth,na.rm=TRUE)))for (i in 1:nn) points(v,av_wealth[,i], col="grey", type = "p")lines(v,rowMeans(av_wealth,na.rm=TRUE), col="black", type = "l")plot(v, rowMeans(better_wealth), xlab="time", ylab="percentage",main="fraction of agent with wealth greater than initial",col="black", type = "l",ylim=c(min(better_wealth,na.rm=TRUE),max(better_wealth,na.rm=TRUE)))for (i in 1:nn) points(v,better_wealth[,i], col="grey", type = "p")lines(v,rowMeans(better_wealth,na.rm=TRUE), col="black", type = "l")plot(v, rowMeans(bankrupts), xlab="time", ylab="number of agents",main="average number of bankrupts", col="black", type = "l",ylim=c(min(bankrupts,na.rm=TRUE),max(bankrupts,na.rm=TRUE)))for (i in 1:nn) points(v, bankrupts[,i], col="grey", type = "p")lines(v,rowMeans(bankrupts,na.rm=TRUE), col="black", type = "l")Исходный код для модели «Искатели Черных лебедей» из подраздела 4.2.3library(fBasics)rm(list=ls(all=TRUE))# загружаем данные по ценамData <- read.table("C:/temp/data_n.txt", header=T)# выставление границ вероятностей р в этой серии экспериментовn1<-0.4n2<-0.6# число экспериментов в серии178nn<-100# критерии успешности стратегийav_wealth<-matrix(0,nrow=(n2-n1)+1,ncol=nn)better_wealth<-matrix(0,nrow=(n2-n1)+1,ncol=nn)bankrupts<-matrix(0,nrow=(n2-n1)+1,ncol=nn)# границы вероятностей определения направления изменения цен для обеихгрупп участников в период стабильной экономики и параметр дельта, отвечающий за изменение этой вероятности в кризисные периоды – у обычныхвероятность уменьшится на дельту, а у искателей Черных лебедей, наоборот, увеличитсяp_min1<-0.4p_max1<-0.5p_min2<-0.5p_max2<-0.6delta<-0.4# величина кредитного рычага (если маржинальна торговля запрещена, то# leverage=0)leverage<-5# достаем данные из матрицыdate <- matrix(0,nrow = tt, ncol=1)price <- matrix(0,nrow = tt, ncol=1)R <- matrix(0,nrow = tt, ncol=1)date <- Data[1:(tt),1 ]price <- Data[1:(tt),2 ]R <- Data[1:(tt),3 ]#начальное благосостояниеcash0<-10000stocks0<-0wealth0<-cash0+stocks0*price[1]#матрицы благосостояний агентовcash<-mat.or.vec(n,tt)stocks<-mat.or.vec(n,tt)wealth<-mat.or.vec(n,tt)#матрицы решений агентовs<-matrix(10, nrow = n, ncol=tt)vol<- matrix(10, nrow = n, ncol=tt)for (ii in 1:nn) {i1<-0for (pr in n1*100:n2*100) {i1<-i1+1# матрица банкротствb<-matrix(0,nrow=2,ncol=tt)cash[,1]<-cash0stocks [,1]<- stocks0wealth[,1]<-wealth0179# матрица со случайными числами для розыгрыша вероятностиu1 <- runif(n*tt)z1 <- matrix(u1, nrow = n, ncol=tt, byrow=TRUE)u2 <- runif(n*tt)z2 <- matrix(u2, nrow = n, ncol=tt, byrow=TRUE)# вектор вероятностей верного определения направления изменения цендля обычных трейдеров и икателей Черных лебедейp<-matrix(0, nrow = n, ncol = 1)p1<-runif(n/2, min=p_min1, max=p_max1)p2<-runif(n/2, min=p_min2, max=p_max2)for (i in 1:(n/2)) p[i]<-p1[i]for (i in (n/2+1):n) p[i]<-p2[i-n/2]m1<-c()m2<-c()k1<-1;k2<-1# счетчик# Блок 1.

Характеристики

Список файлов диссертации