IPB
ЛогинПароль:

> Прочтите прежде чем задавать вопрос!

1. Заголовок темы должен быть информативным. В противном случае тема удаляется ...
2. Все тексты программ должны помещаться в теги [code=pas] ... [/code], либо быть опубликованы на нашем PasteBin в режиме вечного хранения.
3. Прежде чем задавать вопрос, см. "FAQ", если там не нашли ответа, воспользуйтесь ПОИСКОМ, возможно такую задачу уже решали!
4. Не предлагайте свои решения на других языках, кроме Паскаля (исключение - только с согласия модератора).
5. НЕ используйте форум для личного общения, все что не относится к обсуждению темы - на PM!
6. Одна тема - один вопрос (задача)
7. Проверяйте программы перед тем, как разместить их на форуме!!!
8. Спрашивайте и отвечайте четко и по существу!!!

 
 Ответить  Открыть новую тему 
> Комплексные числа, нахождение косинуса
сообщение
Сообщение #1


Новичок
*

Группа: Пользователи
Сообщений: 23
Пол: Мужской
Реальное имя: valera

Репутация: -  0  +


Помогите закончить задачу

Пусть даны комплексное число z (пара вещественных чисел) и вещественное число Эпсилон > 0. Вычислить с точностью Эпсилон значение следующей комплексной функции:
cos z = z – z^2 /2! + z^4 /4! - … + (-1)^n z^2n /(2n)! + …
^ - это степень

program ComplexCos;

type
complex = record
re, im: real

end;
var
Eps: real;
z, res:Complex;

{Функция для определения модуля комплексного числа}
function AbsComplex(z: Complex): real;
begin
AbsComplex := sqrt(str(z.re) + sqr(z. im))
end;

{Процедура сложения двух комплексных чисел}
procedure AddComplex(z1,z2: Complex; var res:Complex);
begin
res.re := z1.re + z2.re;
res.im := z1.im + z2.im;
end;

{процедура умножения двух комплексных чисел}
procedure Multicomplex(z1,z2: Complex; var res:Complex);
begin
res.re := z1.re * z2.re - z1.im * z2.im;
res.im := z1.re * z2.im - z1.im * z2.re
end;

{Процедура деления комплексного числа на вещественное}
procedure DivCoplToReal(var z:Complex; c:real);
begin
if c < ??????????????
then writeln ('На ноль делить нельзя')
else
begin
z.re := z.re / c;
z.im := z.im / c
end
end;

{Процедура вычисления косинуса}
???????????????????????
???????????????????????

end;

begin
{Ввод комплексного числа z}
writeln('Введите комплексное число');
write ('Вещественная часть: ');
readln (z.re);
write ('Комплексная часть: ');
readln (z.im);

{Ввод точности вычеслений}
writeln ('Введите Eps - точность вычеслений' , '(Eps>0)');
repeat
write('Eps=');
readln (Eps);
if Eps <=0 then
writeln('Неправильно.Повторите ввод.');
until (Eps>0);

{Вычисление косинуса}
?????????????
?????????????

{Вывод результатов}
write ('Комплексное число равно', res.re:5:2);
if res.im>=0 then wtite ('+');
write(res.im:5:2, '*i');
readln
end.
 Оффлайн  Профиль  PM 
 К началу страницы 
+ Ответить 
сообщение
Сообщение #2


Новичок
*

Группа: Пользователи
Сообщений: 23
Пол: Мужской
Реальное имя: valera

Репутация: -  0  +


кто нибудь, ну помогите
 Оффлайн  Профиль  PM 
 К началу страницы 
+ Ответить 
сообщение
Сообщение #3


Знаток
****

Группа: Пользователи
Сообщений: 303
Пол: Мужской
Реальное имя: Роман

Репутация: -  2  +


В цикле сравнивай старое значение с новым до тех пор, пока модуль разницы не превысит эпсилон.


--------------------
Romiras HomeLab- материалы и статьи по разработке ПО, моделирование алгоритмов, обработка и анализ информации, нейронные сети, машинное зрение и прочее.
 Оффлайн  Профиль  PM 
 К началу страницы 
+ Ответить 
сообщение
Сообщение #4


Новичок
*

Группа: Пользователи
Сообщений: 23
Пол: Мужской
Реальное имя: valera

Репутация: -  0  +


если б ещё знал как
 Оффлайн  Профиль  PM 
 К началу страницы 
+ Ответить 
сообщение
Сообщение #5


Знаток
****

Группа: Пользователи
Сообщений: 303
Пол: Мужской
Реальное имя: Роман

Репутация: -  2  +


Операторов цикла не знаешь? Учи матчасть.


--------------------
Romiras HomeLab- материалы и статьи по разработке ПО, моделирование алгоритмов, обработка и анализ информации, нейронные сети, машинное зрение и прочее.
 Оффлайн  Профиль  PM 
 К началу страницы 
+ Ответить 
сообщение
Сообщение #6


Новичок
*

Группа: Пользователи
Сообщений: 23
Пол: Мужской
Реальное имя: valera

Репутация: -  0  +


вроде всё!
program ComplexCos;
type
complex = record
re, im: real
end;
var
Eps: real;
z, res:Complex;
{Функция для определения модуля комплексного числа}
function AbsComplex(z: Complex): real;
begin
AbsComplex := sqrt(sqr(z.re) + sqr(z. im))
end;
{Процедура сложения двух комплексных чисел}
procedure AddComplex(z1,z2: Complex; var res:Complex);
begin
res.re := z1.re + z2.re;
res.im := z1.im + z2.im;
end;
{процедура умножения двух комплексных чисел}
procedure Multicomplex(z1,z2: Complex; var res:Complex);
begin
res.re := z1.re * z2.re - z1.im * z2.im;
res.im := z1.re * z2.im - z1.im * z2.re
end;
{Процедура деления комплексного числа на вещественное}
procedure MultCoplToReal(var z:Complex; c:real);
begin
z.re := z.re * c;
z.im := z.im * c
end;
{Процедура вычисления косинуса}
Procedure CosZ(z,a: Complex; e: real; var res: Complex);
var {a: real;}
i: integer;
begin
{write('z=');readln(z);}
{ write('e=');readln(e); }
Multicomplex(z,z,a);
MultCoplToReal(a, -1/2);
i:=1; res:=1+a;
repeat
i:=i+1;
a:=a*(-1)*(sqr(z.re) + sqr(z. im))/(2*i)/(2*i-1);
res:=res+a;
until abs(a)<e;
writeln('cos z = ',res:10:7,' ',cos(z)10.gif7);
readln;
end;
begin
{Ввод комплексного числа z}
writeln('Введите комплексное число');
write ('Вещественная часть: ');
readln (z.re);
write ('Комплексная часть: ');
readln (z.im);
{Ввод точности вычеслений}
writeln ('Введите Eps - точность вычеслений' , '(Eps>0)');
repeat
write('Eps=');
readln (Eps);
if Eps <=0 then
writeln('Неправильно.Повторите ввод.');
until (Eps>0);
{Вычисление косинуса}
{CosZ(z,Eps,res);}
{Вывод результатов}
write ('Комплексное число равно', res.re:5:2);
if res.im>=0 then write ('+');
write(res.im:5:2, '*i');
readln
end.


Результат
Введите комплексное число
Вещественная часть:2
Комплексная часть:1
Введите Eps - точность вычеслений (Eps>0)
Eps=1
Комплексное число равно 0,00+ 0,00*i1

Чёт по моему не то

здесь
Комплексное число равно 0,00+ 0,00*i1
после i она просит ещё раз ввод и потом ставит это число рядом с собой получается 0,00+ 0,00*i1, поставить 2 получается 0,00+ 0,00*i2
и почему всё время получается 0,00+ 0,00*i
 Оффлайн  Профиль  PM 
 К началу страницы 
+ Ответить 
сообщение
Сообщение #7


Гость






А что, возведение комплексного числа в квадрат теперь делается так:
sqr(z.re) + sqr(z. im) ??? И с каких пор? blink.gif
 К началу страницы 
+ Ответить 
сообщение
Сообщение #8


Новичок
*

Группа: Пользователи
Сообщений: 23
Пол: Мужской
Реальное имя: valera

Репутация: -  0  +


если поставить просто sqr (Z) дала ошибку:
функция SQR с параметрами указанный типов не найдена


что за ошибка, вроде всё проверил
 Оффлайн  Профиль  PM 
 К началу страницы 
+ Ответить 
сообщение
Сообщение #9


Гость






Правильно, потому что Complex - НЕ встроенный тип Турбо Паскаля... А поэтому делаем по-другому: допустим, имеем комплексное число Z = (a + b*i)... Возведем его в квадрат:

sqr(z) = z * z = (a + b*i) * (a + b*i) = a^2 + a*b*i + a*b*i + (b*i)^2 = a^2 + 2*a*b*i - b^2 = (a^2 – b^2) + 2*a*b*i

А что у тебя получилось?
 К началу страницы 
+ Ответить 
сообщение
Сообщение #10


Новичок
*

Группа: Пользователи
Сообщений: 23
Пол: Мужской
Реальное имя: valera

Репутация: -  0  +


мне её надо где-то вписать в код sqr(z)?
немножко не понял вопроса

Сообщение отредактировано: valeral -
 Оффлайн  Профиль  PM 
 К началу страницы 
+ Ответить 
сообщение
Сообщение #11


Новичок
*

Группа: Пользователи
Сообщений: 23
Пол: Мужской
Реальное имя: valera

Репутация: -  0  +


Volvo обьясни пожалуйста
 Оффлайн  Профиль  PM 
 К началу страницы 
+ Ответить 
сообщение
Сообщение #12


Гость






Значит, так... Вот рабочая программа для вычисления значения Cos(z) комплексного аргумента разложением в ряд:
Program ZnachCos;

type
complex = record
re, im: real;
end;

procedure GetComplex(var cplx: Complex);
begin
write('(re) part: '); readln(cplx.re);
write('(im) part: '); readln(cplx.im);
end;
procedure PutComplex(var cplx: Complex);
const sign: array[boolean] of char = ('+', '-');
begin
write('(', cplx.re:7:4, sign[cplx.im < 0], abs(cplx.im):7:4, '*i');
end;

function AbsComplex(z: Complex): real;
begin
AbsComplex := sqrt(sqr(z.re) + sqr(z.im))
end;

procedure AddComplex(z1,z2: Complex; var res:Complex);
begin
res.re := z1.re + z2.re;
res.im := z1.im + z2.im;
end;

procedure Multicomplex(z1,z2: Complex; var res:Complex);
begin
res.re := z1.re * z2.re - z1.im * z2.im;
res.im := z1.re * z2.im + z1.im * z2.re;
end;
procedure ScaleComplex(var z:Complex; c:real);
begin
z.re := z.re * c;
z.im := z.im * c
end;

Var
n,k:Integer;
b,e,p,Eps:real;
z, a, s, sq: complex;
Begin
writeln('complex Z = '); GetComplex(z);
Write('real Eps = '); Readln(Eps);

n := 2;
a.re := 1; a.im := 0;
s.re := 1; s.im := 0;
MultiComplex(z, z, sq);
While AbsComplex(a) > Eps do Begin
MultiComplex(a, sq, a);
ScaleComplex(a, - 1 / (pred(n) * n));
AddComplex(s, a, s);
Inc(n, 2)
End;
Writeln('COS(z) = '); PutComplex(s);
End.


Теперь - разбор твоих ошибок (может, еще кому интересно будет)... Первая твоя ошибка - в формуле. Правильное разложение косинуса:
cos z = 1 – z^2 /2! + z^4 /4! - … ,
а не z как у тебя...

Вторая - в реализации операции умножения комплексных чисел... Ты поставил знак "-" вместо нужного плюса:
procedure Multicomplex(z1,z2: Complex; var res:Complex);
begin
res.re := z1.re * z2.re - z1.im * z2.im;
res.im := z1.re * z2.im + z1.im * z2.re; { <--- Здесь у тебя был МИНУС }
end;


Ну, и третье - смотри свою реализацию в repeat ... until, и мою в while. Как говорится, найди 10 отличий...

Проверялось на Z = (2, 1)
при Eps = 0.001
Ответ совпадает с тем, что получено на С++ с использованием <complex.h>
 К началу страницы 
+ Ответить 
сообщение
Сообщение #13


Новичок
*

Группа: Пользователи
Сообщений: 39
Пол: Мужской
Реальное имя: leahov

Репутация: -  0  +


здесь выдало ошибку
const sign: array[boolean] of char = ('+', '-');

[ boolean ] - ожидалось выражение
 Оффлайн  Профиль  PM 
 К началу страницы 
+ Ответить 
сообщение
Сообщение #14


Новичок
*

Группа: Пользователи
Сообщений: 39
Пол: Мужской
Реальное имя: leahov

Репутация: -  0  +


звиняйте в турбо паскале всё работает! Всём спасибо! Volvo - отдельный благодарчик
 Оффлайн  Профиль  PM 
 К началу страницы 
+ Ответить 

 Ответить  Открыть новую тему 
1 чел. читают эту тему (гостей: 1, скрытых пользователей: 0)
Пользователей: 0

 




- Текстовая версия 21.08.2017 9:44
Хостинг предоставлен компанией "Веб Сервис Центр" при поддержке компании "ДокЛаб"