dmi指标源码(高手改写个dmi选股公式？)
1. 高手改写个dmi选股公式？
N:=14;M:=6;MTR:= SUM(MAX(MAX(HIGH-LOW,ABS(HIGH-REF(CLOSE,1))),ABS(LOW-REF(CLOSE,1))),N)
;HD := HIGH-REF(HIGH,1)
;LD := REF(LOW,1)-LOW;DMP:= SUM(IF(HD>0 AND HD>LD,HD,0),N)
;DMM:= SUM(IF(LD>0 AND LD>HD,LD,0),N)
;PDI:=DMP*100/MTR;MDI:=DMM*100/MTR;ADX:=MA(ABS(MDI-PDI)/(MDI+PDI)*100,M)
;ADXR:=(ADX+REF(ADX,M))/2;XG:PDI>ADX AND ADX>ADXR AND CROSS(PDI,30);(已测试)
2. 32位浮点数乘法汇编算法？
写汇编的基本已经绝迹了现在,就老一辈的有些在用
#include <stdio.h>
#include <sys/time.h>
#define INIT_TIMER_VALIABLE
struct timeval tpstart,tpend;
float timeuse;
#define START_TIMER gettimeofday(&tpstart,NULL);
#define END_PRINTF_TIMER(name)
gettimeofday(&tpend,NULL);
timeuse=(tpend.tv_sec*1000*1000+tpend.tv_usec)-(tpstart.tv_sec*1000*1000+tpstart.tv_usec);
printf("func :%s:time use(us) %fn",name,timeuse);
float vfp_operate(float f1, float f2)
{
float sum=0;
__asm__ __volatile__(
"vmov s1, %1n"
"vmov s2, %2n"
"fmuls s0, s1, s2n"
"vmov %0, s0n"
:"=r"(sum)
:"r"(f1),"r"(f2)
);
return sum;
}
int main()
{
float f1,f2;
float result;
INIT_TIMER_VALIABLE
printf("input float data1:");
scanf("%f",&f1);
printf("input float data2:");
scanf("%f",&f2);
START_TIMER
result=f1*f2;
END_PRINTF_TIMER("use system function")
printf("result is %fn",result);
START_TIMER
result=vfp_operate(f1,f2);
END_PRINTF_TIMER("use vfp_operate")
printf("result is %fn",result);
}
/*
测试结果：
./vfp_helloworld
input float data1:0.125
input float data2:1.684
func :use system function:time use(us) 12.000000
result is 0.210500
func :use system function:time use(us) 5.000000
result is 0.210500
*/
测试使用编译选项 ：
源代码：
#include <stdio.h>
int main()
{
float f1,f2;
printf("input float data1:");
scanf("%f",&f1);
printf("input float data2:");
scanf("%f",&f2);
printf("float %f x %f =%fn", f1,f2,f1*f2);
}
编译选项：(只是多增加了-mfloat-abi=softfp 测试发现=sofltfp和=hard编译出的汇编一样)
arm-none-linux-gnueabi-gcc vfp_helloworld.c -S arm-none-linux-gnueabi-gcc vfp_helloworld.c -S -o vfp_helloworld.asm
-mfloat-abi=softfp -o vfp_helloworld.asm
.cpu arm10tdmi .cpu arm10tdmi
.eabi_attribute 27, 3 .fpu softvfp
.fpu vfp .eabi_attribute 20, 1
.eabi_attribute 20, 1 .eabi_attribute 21, 1
.eabi_attribute 21, 1 .eabi_attribute 23, 3
.eabi_attribute 23, 3 .eabi_attribute 24, 1
.eabi_attribute 24, 1 .eabi_attribute 25, 1
.eabi_attribute 25, 1 .eabi_attribute 26, 2
.eabi_attribute 26, 2 .eabi_attribute 30, 6
.eabi_attribute 30, 6 .eabi_attribute 18, 4
.eabi_attribute 18, 4 .file "vfp_helloworld.c"
.file "vfp_helloworld.c" .section .rodata
.section .rodata .align 2
.align 2 .LC0:
.LC0: .ascii "input float data1: 00"
.ascii "input float data1: 00" .align 2
.align 2 .LC1:
.LC1: .ascii "%f 00"
.ascii "%f 00" .align 2
.align 2 .LC2:
.LC2: .ascii "input float data2: 00"
.ascii "input float data2: 00" .global __aeabi_f2d
.align 2 .global __aeabi_fmul <<<------这里说明软件浮点运算。
.LC3: .align 2
.ascii "float %f x %f =%f 12 00" .LC3:
.text .ascii "float %f x %f =%f 12 00"
.align 2 .text
.global main .align 2
.type main, %function .global main
main: .type main, %function
.fnstart main:
.LFB2: .fnstart
@ args = 0, pretend = 0, frame = 16 .LFB2:
@ frame_needed = 1, uses_anonymous_args = 0 @ args = 0, pretend = 0, frame = 16
stmfd sp!, {fp, lr} @ frame_needed = 1, uses_anonymous_args = 0
.save {fp, lr} stmfd sp!, {r4, r5, r6, r7, r8, fp, lr}
.LCFI0: .save {r4, r5, r6, r7, r8, fp, lr}
.setfp fp, sp, #4 .LCFI0:
add fp, sp, #4 .setfp fp, sp, #24
.LCFI1: add fp, sp, #24
.pad #32 .LCFI1:
sub sp, sp, #32 .pad #36
.LCFI2: sub sp, sp, #36
ldr r0, .L3 .LCFI2:
bl printf ldr r0, .L3
sub r3, fp, #8 bl printf
ldr r0, .L3+4 sub r3, fp, #32
mov r1, r3 ldr r0, .L3+4
bl scanf mov r1, r3
ldr r0, .L3+8 bl scanf
bl printf ldr r0, .L3+8
sub r3, fp, #12 bl printf
ldr r0, .L3+4 sub r3, fp, #36
mov r1, r3 ldr r0, .L3+4
bl scanf mov r1, r3
flds s15, [fp, #-8] bl scanf
fcvtds d5, s15 ldr r3, [fp, #-32] @ float
flds s15, [fp, #-12] mov r0, r3
fcvtds d6, s15 bl __aeabi_f2d
flds s14, [fp, #-8] mov r5, r0
flds s15, [fp, #-12] mov r6, r1
fmuls s15, s14, s15 <<<-----直接使用硬件浮点指令 ldr r3, [fp, #-36] @ float
fcvtds d7, s15 mov r0, r3
fstd d6, [sp, #0] bl __aeabi_f2d
fstd d7, [sp, #8] mov r7, r0
ldr r0, .L3+12 mov r8, r1
fmrrd r2, r3, d5 ldr r3, [fp, #-32] @ float
bl printf ldr r2, [fp, #-36] @ float
sub sp, fp, #4 mov r0, r3
ldmfd sp!, {fp, pc} mov r1, r2
.L4: bl __aeabi_fmul <<<------这里调用软件浮点运算。
.align 2 mov r3, r0
.L3: mov r0, r3
.word .LC0 bl __aeabi_f2d
.word .LC1 mov r3, r0
.word .LC2 mov r4, r1
.word .LC3 stmia sp, {r7-r8}
.LFE2: str r3, [sp, #8]
.fnend str r4, [sp, #12]
.size main, .-main ldr r0, .L3+12
.ident "GCC: (Sourcery G++ Lite 2009q1-203) 4.3.3" mov r2, r5
.section .note.GNU-stack,"",%progbits mov r3, r6
bl printf
sub sp, fp, #24
ldmfd sp!, {r4, r5, r6, r7, r8, fp, pc}
.L4:
.align 2
.L3:
.word .LC0
.word .LC1
.word .LC2
.word .LC3
.LFE2:
.fnend
.size main, .-main
.ident "GCC: (Sourcery G++ Lite 2009q1-203) 4.3.3"
.section .note.GNU-stack,"",%progbits
摘自：RM的pdf文档 的说明
浮点运算的支持
ARM 处理器内核不包含浮点硬件。 必须使用以下两种方法之一，另行提供对浮点算法的支持：
在软件中，使用浮点库 fplib。 此库提供了执行浮点运算可以调用的函数，无需额外的硬件。 请参阅《库指南》 中第 4-2 页的软件浮点库 fplib。
在硬件中，使用含 VFP 硬件协处理器的 ARM 处理器内核来进行所需的浮点运算。 VFP 是执行 IEEE 浮点的协处理器体系结构，支持单精度和双精度，但不支持扩展精度。
Note
在实际编程中，VFP 中的浮点运算实际是组合使用硬件（执行常见的情况）和软件（处理不常见的情况和导致异常的情况）执行的。 请参阅VFP 支持。
Example 5.2 是一个用 C 执行浮点算法的函数，用以说明浮点算法的软件和硬件支持的不同。
Example 5.2. 浮点运算
float foo(float num1, float num2)
{
float temp, temp2;
temp = num1 + num2;
temp2 = num2 * num2;
return temp2-temp;
}
如果使用命令行选项 --cpu 5TE --fpu softvfp 编译Example 5.2 的 C 代码，则编译器生成的机器代码的反汇编如Example 5.3 所示。 在本示例中，在软件中通过调用库例程（如 __aeabi_fmul）来执行浮点算法。
Example 5.3. 软件中对浮点运算的支持
||foo|| PROC
PUSH {r4-r6, lr}
MOV r4, r1
BL __aeabi_fadd <<<<----直接向加
MOV r5, r0
MOV r1, r4
MOV r0, r4
BL __aeabi_fmul <<<<<----然后相乘
MOV r1, r5
POP {r4-r6, lr}
B __aeabi_fsub
ENDP
如果使用命令行选项 --fpu vfp 编译Example 5.2 的 C 代码，则编译器生成的机器代码的反汇编如Example 5.4 所示。 在本示例中，在硬件中通过浮点算法指令（如 VMUL.F32）来执行浮点算法。
Example 5.4. 硬件中对浮点运算的支持
||foo|| PROC
VADD.F32 s2, s0, s1
VMUL.F32 s0, s1, s1
VSUB.F32 s0, s0, s2
BX lr
ENDP
在实际编程中，使用硬件支持浮点算法的代码更为紧凑，并提供比在软件中执行浮点算法的代码更佳的性能。 但是，浮点算法的硬件支持需要 VFP 协处理器。
缺省情况下，如果有 VFP 协处理器，则会生成 VFP 指令。 如果没有 VFP 协处理器，则编译器会生成调用软件浮点库 fplib 的代码，用于执行浮点运算。fplib 是 C 库 RealView Development Suite 标准分发的组成部分。
3. dmi趋向指标公式？
股票走势的技术指标中，有趋势类、摆动类、能量类的指标，而DMI就是其中的一种经典的趋势指标，在实战中能够帮助投资者识别价格的趋势变化以及其中蕴含的买卖点位置，为了能够更加清晰的了解指标的应用技巧，我们从其源码公式进行学习。
首先我们来看一下具有顶底识别以及抄底逃顶作用的指标源码的编译：
【顶底DMI】
TR:=EMA(MAX(MAX(HIGH-LOW,ABS(HIGH-REF(CLOSE,1))),ABS(REF(CLOSE,1)-LOW)),7);
HD :=HIGH-REF(HIGH,1);
LD :=REF(LOW,1)-LOW;
DMP:=EMA(IF(HD>0 AND HD>LD,HD,0),7);
DMM:=EMA(IF(LD>0 AND LD>HD,LD,0),7);
PDI: DMP*100/TR;
MDI: DMM*100/TR;
ADX: EMA(ABS(MDI-PDI)/(MDI+PDI)*100,7);
ADXR:EMA(ADX,7);
抄底1:COUNT(PDI50,1) AND REF(ADX,1)>=REF(ADX,2) AND
ADX
逃顶1:=COUNT(PDI>MDI,3)=3 AND REF(ADX>50,1) AND REF(ADX,1)>=REF(ADX,2) AND
ADX
五十:50;
三十:30;
二十:20;
DRAWTEXT(抄底1,PDI,'抄底');
DRAWTEXT(逃顶1,ADX,'逃顶');
顶底dmi
从源码和走势图形中我们可以看出，其中指标线包含有PDI、ndi、adx等走势，具体的操作用法，请点击DMI指标，在这里主要讲解的是其源码公式的编译，以及简单的走势分析。
其次，我们可以通过其源码公式的展示，寻找指标线中的买入点位置：
【dmi买点】
{参数n14,m6}
TR:=EXPMEMA(MAX(MAX(HIGH-LOW,ABS(HIGH-REF(CLOSE,1))),ABS(REF(CLOSE,1)-LOW)),N);
HD :=HIGH-REF(HIGH,1);
LD :=REF(LOW,1)-LOW;
DMP:=EXPMEMA(IF(HD>0&&HD>LD,HD,0),N);
DMM:=EXPMEMA(IF(LD>0&&LD>HD,LD,0),N);
PDI: DMP*100/TR;
MDI: DMM*100/TR;
ADX: EXPMEMA(ABS(MDI-PDI)/(MDI+PDI)*100,M);
穿破短底:16;
线上发飙:30;
{2}预警:IF(ADX>90 AND PDI=0,30,0),COLORYELLOW,LINETHICK2;
买:IF(CROSS(PDI,MDI) AND (C >1.05*REF(C,1) OR REF(C,1)>1.05*REF(C,2)) AND (V>2*REF(V,1) OR REF(V,1)>2*REF(V,2))AND MA(V,5)>MA(V,60) ,20,0),COLORRED,LINETHICK2;
买2:IF(CROSS(PDI,MDI) AND COUNT(预警,10)>=1,10,0),COLORMAGENTA,LINETHICK3;
0;
{6}
RSV9:=(CLOSE-LLV(LOW,9))/(HHV(HIGH,9)-LLV(LOW,9))*100;
K1:=SMA(RSV9,3,1);
DK1:=K1-REF(K1,1),COLORSTICK;
出手就赢:IF(DK1>7 AND REF(DK1,1)<0 AND (DK1-REF(DK1,1))>9,50,0), STICK,COLORYELLOW;DRAWICON(出手就赢>0,50 ,1);
dmi买点
依据其指标线的走势形态，在股价的K线走势图中能够清楚的表现出，投资者应该买入以及卖出的获利点位置，希望投资者能够从中做好应对的操作策略。
【dmi组合】
VAR1:=(C-LLV(L,13))/(HHV(H,13)-LLV(L,13))*100;
牛:SMA(VAR1,7,1);
VAR2:=SMA(牛,5,1);
熊:IF(CROSS(牛,VAR2) AND CROSS(MA(C,5),MA(C,13) AND 牛>50),VAR2*2,VAR2);
DRAWBAND(牛,RGB(205,92,92),熊,RGB(28,134,238));
N:=14;M:=6;
TR:=EXPMEMA(MAX(MAX(HIGH-LOW,ABS(HIGH-REF(CLOSE,1))),ABS(REF(CLOSE,1)-LOW)),N);
HD :=HIGH-REF(HIGH,1);
LD :=REF(LOW,1)-LOW;
DMP:=EXPMEMA(IF(HD>0&&HD>LD,HD,0),N);
DMM:=EXPMEMA(IF(LD>0&&LD>HD,LD,0),N);
PDI:DMP*100/TR COLORMAGENTA;
MDI:=DMM*100/TR COLORGREEN;
ADX:=EXPMEMA(ABS(MDI-PDI)/(MDI+PDI)*100,M);
ADXR:=EXPMEMA(ADX,M);
STICKLINE(ADX>=ADXR,ADX,ADXR,3,0),COLOR800040;
STICKLINE(ADX
关注:STICKLINE(BARSLAST(CROSS(MDI,PDI))<5 AND MDI>PDI AND ADXADXR*0.97,ADX,ADXR,2,0) COLORRED;
10;90;
FVAR1:=REF(CLOSE,2);
会员专用: SMA(MAX(CLOSE-FVAR1,0),7,1)/SMA(ABS(CLOSE-FVAR1),7,1)*100,COLORYELLOW;
{一、PDI>MDI且趋势向上;二、ADX;ADXR在PDI及MDI中间金叉 }
走强:DRAWTEXT(CROSS(ADX,ADXR) AND PDI>MDI AND PDI>REF(PDI,1) AND ADX>REF(ADX,1) AND BETWEEN(ADX,MDI,PDI) AND ADXR>REF(ADXR,1),ADXR*0.98,'走强'),COLORYELLOW;
