Координаты точек по окаймлению профиля

@Merhaba · Регистрация: 05.03.2011

Author24 — интернет-сервис помощи студентам

Доброго времени суток! Помогите Пожалуйста!
Данный код выводит на экран изображения руля Жуковского:

Fortran

type (windowconfig) wc
 
 d=GETWINDOWCONFIG(wc); p=wc.numxpixels; q=wc.numypixels
 step_px=2.0*Xmax/p; step_qy=2.0*Ymax/q
 
 !call number_pi (pi, pi_2, RAD)
 
    do i=0,p,1
        do j=0,q,1
        z=dcmplx(dble(+Xmin+i*step_px), dble(Ymin+j*step_qy))
       dzeta=conform_Joukowsky(z)/R_
        if (cdabs(dzeta) <=1 ) status=setpixel(i,j)
        end do
    end do

Мне нужно, чтобы программа выдавала координаты точек по окаймлению профиля! Как это сделать?

@Merhaba · 28.03.2014, 15:46 **[ТС]**

Krasme, это три вихря с заданной циркуляцией!

@Merhaba · 29.03.2014, 16:12 **[ТС]**

Krasme, как можно к профилю пририсовать прямую линию в чёрной точке, чтобы длину линии и угол её отклонения от горизонтального положения мы задавали сами?! Если Вам не сложно, помогите Пожалуйста мне это сделать!

@Krasme · 29.03.2014, 16:20

Merhaba, прям целый курсовик...
подумаю на досуге. может быть...

улетаю на профиле "Метла"

@Krasme · 29.03.2014, 19:36

Merhaba, держите

Кликните здесь для просмотра всего текста

Fortran

        status = setcolor(12)
    !do ib=iminb,imaxb
    !  status=setpixel(ib,jb)
    !end do
    !do jc=jminc,jmaxc
    !  status=setpixel(ic,jc)
    !end do
    call moveto(iminb,jb,xys)
    d=lineto(imaxb,jb);
    call moveto(ic,jminc,xys)
    d=lineto(ic,jmaxc);
    !=======================
    ! line with angle
    angle=-10*3.14/180;
    dlina=400;
    status = setcolor(0)
    call moveto(iminb,jb,xys)
    d=lineto(iminb+dlina*cos(angle),jb+dlina*sin(angle))
 
    !=======================

обратите внимание на закомментированные строки, ниже новый (более удобный) вариант.
положительный угол по часовой стрелке.

@Merhaba · 30.03.2014, 09:00 **[ТС]**

Krasme, Выдаёт:

Fortran

1	Error: There is no matching specific subroutine for this generic subroutine call. [MOVETO]

@Krasme · 30.03.2014, 09:05

Merhaba, вот весь текст, смотрите, чего не хватает

Кликните здесь для просмотра всего текста

Fortran

subroutine Paint_foil_Gukovsky (Xmax,Ymax,Xmin,Ymin)
 
use msflib
use foil_J
use param
 
!---------------------------------------------------------
INTERFACE
DOUBLE complex FUNCTION conform_Joukowsky(z) result(dzeta)
 
 
use foil_J
use point_
use param
 
complex(8),INTENT(in)::z 
END FUNCTION conform_Joukowsky
END INTERFACE
!---------------------------------------------------------
real(8),intent(in)::Xmin,Ymin, Xmax,Ymax
!-------------------------------------------------------------
logical(4)::d
complex(8)::dzeta, z
integer(4)::i,j,p,q
real(4)   ::step_px, step_qy
 
type (windowconfig) wc
type( XYCOORD ) xys
 d=GETWINDOWCONFIG(wc); p=wc.numxpixels; q=wc.numypixels
 step_px=2.0*Xmax/p; step_qy=2.0*Ymax/q
 
 call number_pi (pi, pi_2, RAD)
 dj=0;di=0; imax=0;imin=p;
    do i=0,p,1
        jmax=0;jmin=q;
        do j=0,q,1
         z=dcmplx(dble(Xmin+i*step_px), dble(Ymin+j*step_qy))
         dzeta=conform_Joukowsky(z)/R_
         if (cdabs(dzeta) <=1.and.cdabs(dzeta)>0.99 ) then
            status=setpixel(i,j)
            !--------------------
            if (jmax<j) jmax=j;             if (jmin>j) jmin=j;
            if (dj<jmax-jmin) then
               dj=jmax-jmin;
               jmaxc=jmax;   jminc=jmin;  ic=i;
            end if
            !--------------------
            if (imax<i) imax=i;             if (imin>i) imin=i;
            if (di<imax-imin) then
               di=imax-imin;
               imaxb=imax;   iminb=imin;  jb=j;
            end if
         end if
        end do
    end do
    !-----------------
   write (222,*) dj,ic,jmaxc,jminc
   write (222,*) di,jb,imaxb,iminb
    status = setcolor(12)
    !do ib=iminb,imaxb
    !  status=setpixel(ib,jb)
    !end do
    !do jc=jminc,jmaxc
    !  status=setpixel(ic,jc)
    !end do
    call moveto(iminb,jb,xys)
    d=lineto(imaxb,jb);
    call moveto(ic,jminc,xys)
    d=lineto(ic,jmaxc);
    !=======================
    ! line with angle
    angle=-10*3.14/180;
    dlina=400;
    status = setcolor(0)
    call moveto(iminb,jb,xys)
    d=lineto(iminb+dlina*cos(angle),jb+dlina*sin(angle))
 
    !=======================
return
end subroutine Paint_foil_Gukovsky

@Merhaba · 30.03.2014, 09:42 **[ТС]**

Krasme, попробую написать программу обтекания профиля с такой штукой...

Добавлено через 2 минуты
Krasme, как-нить можно узнавать координаты точек этой прямой?

@Krasme · 30.03.2014, 09:45

Сообщение от Merhaba

как-нить можно узнавать координаты точек этой прямой?

а эта формула ни о чем не говорит?

Fortran

1	d=lineto(iminb+dlinacos(angle),jb+dlinasin(angle))

@Merhaba · 30.03.2014, 09:49 **[ТС]**

Krasme, отсюда мы получим координаты точек конца прямой!

@Krasme · 30.03.2014, 09:53

Merhaba, ну зачем так буквально читать? это наводка..
а код можно записать след.образом:

Fortran

do dl=1,dlina
idl=iminb+dl*cos(angle)
jdl=jb+dl*sin(angle)
end do

написано навскидку, без проверки.. сами, сами, сами!

@Merhaba · 30.03.2014, 21:18 **[ТС]**

Krasme, как можно сделать так, чтобы прямая появлялась и исчезала по нашему требованию?
Мне кажется, что тут нужна логическая переменная.
Пытался написать такой код, но программа с ним не работает:

Fortran

logical :: flap 
if flap =.TRUE. then
call moveto(iminb,jb,xys)
    d=lineto(iminb+dlina*cos(angle),jb+dlina*sin(angle))
end if

@Krasme · 30.03.2014, 21:23

Fortran

logical :: flap =.true.
...
if (flap) then
call moveto(iminb,jb,xys)
    d=lineto(iminb+dlina*cos(angle),jb+dlina*sin(angle))
end if

@Merhaba · 12.04.2014, 21:22 **[ТС]**

Krasme, как можно найти координаты точки "прикрепления" прямой линии к профилю в системе координат wxycoord ?

@Krasme · 12.04.2014, 21:28

Merhaba,

Сообщение от Merhaba

координаты точки "прикрепления" прямой линии к профилю

не понимаю...

@Merhaba · 12.04.2014, 21:47 **[ТС]**

Krasme, координату черной точки в системе координат wxycoord

Добавлено через 7 минут
Krasme, call moveto(iminb,jb,xys)

@Krasme · 12.04.2014, 22:18

Merhaba, на цап-царап..

Fortran

1	call getwindowcoord(iminb,jminc,wxys)

Не по теме:

никогда не занималась различными системами координат, потому ни о каких преобразованиях понятия не имею :mda:

@Merhaba · 17.04.2014, 17:39 **[ТС]**

Krasme, А можно как-нить вместо прямой линии, от точки (iminb,jb), нарисовать сигмент профиля, выделенный черным цветом и управлять углом его отклонения?

@Krasme · 17.04.2014, 18:09

Сообщение от Merhaba

А можно как-нить вместо прямой линии, от точки (iminb,jb), нарисовать сигмент профиля, выделенный черным цветом

думаю, что можно...

Сообщение от Merhaba

управлять углом его отклонения

не понимаю, что за угол отклонения..

@Merhaba · 17.04.2014, 18:22 **[ТС]**

Krasme, Для линии мы сами задаём угол отклонения - angle

Fortran

1 2	call moveto(iminb,jb,xys) d=lineto(iminb+dlinacos(angle),jb+dlinasin(angle))

Добавлено через 1 минуту
Krasme, у Вас есть какие-нить мысли, как это можно сделать?

@Krasme · 17.04.2014, 18:28

Merhaba, навскидку - дополнительный двумерный массив.. но это будет на основе полученного профиля, но чтобы отдельно, да еще и угол менять.. не понимаю как.. та линия хоть строится самостоятельно.. работа уже не ограничивается одной графикой, надо вникать в полную задачу.

мне не очень-то хочется

@Merhaba 3 / 3 / 5 Регистрация: 05.03.2011 Сообщений: 270
	30.03.2014, 09:42 [ТС]	27
	Krasme, попробую написать программу обтекания профиля с такой штукой... Добавлено через 2 минуты Krasme, как-нить можно узнавать координаты точек этой прямой? 0

@Merhaba 3 / 3 / 5 Регистрация: 05.03.2011 Сообщений: 270
	12.04.2014, 21:22 [ТС]	33
	Krasme, как можно найти координаты точки "прикрепления" прямой линии к профилю в системе координат wxycoord ? 0

@Krasme 6830 / 4890 / 2065 Регистрация: 02.02.2014 Сообщений: 13,048
	12.04.2014, 21:28	34
	Merhaba, Сообщение от Merhaba координаты точки "прикрепления" прямой линии к профилю не понимаю... 0

@Merhaba 3 / 3 / 5 Регистрация: 05.03.2011 Сообщений: 270
	12.04.2014, 21:47 [ТС]	35
	Krasme, координату черной точки в системе координат wxycoord Добавлено через 7 минут Krasme, call moveto(iminb,jb,xys) 0

@Merhaba 3 / 3 / 5 Регистрация: 05.03.2011 Сообщений: 270
	17.04.2014, 17:39 [ТС]	37
	Krasme, А можно как-нить вместо прямой линии, от точки (iminb,jb), нарисовать сигмент профиля, выделенный черным цветом и управлять углом его отклонения? 0

@Merhaba 3 / 3 / 5 Регистрация: 05.03.2011 Сообщений: 270
	28.03.2014, 15:46 [ТС]	21
	Krasme, это три вихря с заданной циркуляцией! 0

@Merhaba 3 / 3 / 5 Регистрация: 05.03.2011 Сообщений: 270
	29.03.2014, 16:12 [ТС]	22
	Krasme, как можно к профилю пририсовать прямую линию в чёрной точке, чтобы длину линии и угол её отклонения от горизонтального положения мы задавали сами?! Если Вам не сложно, помогите Пожалуйста мне это сделать! 0

@Krasme 6830 / 4890 / 2065 Регистрация: 02.02.2014 Сообщений: 13,048
	29.03.2014, 16:20	23
	Сообщение было отмечено Merhaba как решение Решение Merhaba, прям целый курсовик... подумаю на досуге. может быть... улетаю на профиле "Метла" 1

@Merhaba 3 / 3 / 5 Регистрация: 05.03.2011 Сообщений: 270
	30.03.2014, 09:49 [ТС]	29
	Krasme, отсюда мы получим координаты точек конца прямой! 0

@Krasme 6830 / 4890 / 2065 Регистрация: 02.02.2014 Сообщений: 13,048
	17.04.2014, 18:09	38
	Сообщение от Merhaba А можно как-нить вместо прямой линии, от точки (iminb,jb), нарисовать сигмент профиля, выделенный черным цветом думаю, что можно... Сообщение от Merhaba управлять углом его отклонения не понимаю, что за угол отклонения.. 0

@Krasme 6830 / 4890 / 2065 Регистрация: 02.02.2014 Сообщений: 13,048
	17.04.2014, 18:28	40
	Merhaba, навскидку - дополнительный двумерный массив.. но это будет на основе полученного профиля, но чтобы отдельно, да еще и угол менять.. не понимаю как.. та линия хоть строится самостоятельно.. работа уже не ограничивается одной графикой, надо вникать в полную задачу. мне не очень-то хочется 1

Координаты точек по окаймлению профиля

Решение