Список форумов Форум космопорта Форум космопорта
Космос и все, что с ним связано...
 
 FAQFAQ   ПоискПоиск   ПользователиПользователи   ГруппыГруппы   РегистрацияРегистрация 
 ПрофильПрофиль   Войти и проверить личные сообщенияВойти и проверить личные сообщения   ВходВход 

Периметр эллипса
На страницу Пред.  1, 2, 3  След.
 
Начать новую тему   Ответить на тему    Список форумов Форум космопорта -> Поговорим?
Предыдущая тема :: Следующая тема  
Автор Сообщение
Гость






СообщениеДобавлено: 14 Июн 2005 [22:48]    Заголовок сообщения: Ответить с цитатой

и еще. было несколько некорректно называть вашу новую величину средним эллиптическим. Средняя величина подразумевет наличие какого-то ряда, будь то арифметическая, геометрическая прогрессии или гармонический ряд. Среднее вводится как метод определения членая An, при известных An-1 и An+1. А эллиптического ряда не вводили. поэтому это средний радиус. и поэтому ваше утверждение о том. что средние величины при почти одинкаовых a и b теряют между собой различие неверно, так как средняя эллиптическая не удовлетворяет определению в обычном смысле.

С уважением Андрюша Мехедов
Вернуться к началу
Непризнанный Гений
Ветеран


Зарегистрирован: 27.02.2005
Сообщения: 606

СообщениеДобавлено: 15 Июн 2005 [03:20]    Заголовок сообщения: Ответить с цитатой

Вы Андрей прям как великий Леонардо!

Ваша программа не компилится Sad

Код:

In file included from /usr/include/c++/3.2.2/backward/iostream.h:31,
                 from ellips.cxx:1:
/usr/include/c++/3.2.2/backward/backward_warning.h:32:2: warning: #warning This file includes at least one deprecated or antiquated header. Please consider using one of the 32 headers found in section 17.4.1.2 of the C++ standard. Examples include substituting the <X> header for the <X.h> header for C++ includes, or <sstream> instead of the deprecated header <strstream.h>. To disable this warning use -Wno-deprecated.
ellips.cxx: In function `double Intf(double)':
ellips.cxx:35: no matching function for call to `f(double, double&)'
ellips.cxx:10: candidates are: double f(double)
ellips.cxx:28:                 double f(double, double, double)
ellips.cxx:35: no matching function for call to `f(double&, double&)'
ellips.cxx:10: candidates are: double f(double)
ellips.cxx:28:                 double f(double, double, double)
Вернуться к началу
Посмотреть профиль Отправить личное сообщение
Непризнанный Гений
Ветеран


Зарегистрирован: 27.02.2005
Сообщения: 606

СообщениеДобавлено: 15 Июн 2005 [03:54]    Заголовок сообщения: Ответить с цитатой

В функции f(x,y) что есть y?

x это t, как я понимаю.

А y - это отношение a к b?
Вернуться к началу
Посмотреть профиль Отправить личное сообщение
Непризнанный Гений
Ветеран


Зарегистрирован: 27.02.2005
Сообщения: 606

СообщениеДобавлено: 15 Июн 2005 [04:08]    Заголовок сообщения: Ответить с цитатой

Андрей, мне пришлось подправить Вашу программу, чтобы она заработала.

Теперь она такая :-

cat ellips.cxx
Код:

#include <iostream.h>
#include <fstream.h>
#include <math.h>

#define N 200000
#define Pi 3.141592653
#define X 1.534928536

double Intf(double y);
double f(double x);

using namespace std;

int main() {
  fstream file;
  file.open("dat.txt",ios::out);
  double x, i, s;
  for(int n=0;n<=100;n++) {
    x=(double)0.05*n;
    cout << "I[0,Pi/2] f(" << x << ") dx =";
    i=Intf(x);
    s=0.5*Pi*pow(0.5*(1+pow(x,X)), 1.0/X);
    file << x << " " << i << " " << s << " " << 0.5*Pi*sqrt(x+4*(1-x)*(1-x)/(Pi*Pi)) << endl;
    cout << i << "\t" << s << endl;
   
  }
  return 0;
}

double f(double x, double y)
{
return sqrt(sin(x)*sin(x)+y*y*cos(x)*cos(x));//!!!!!!!!!!!!!!!!!!!!!!!!!!!!
}

double Intf(double y) {
  double x=0, i=0, dx=0.5*Pi/N;
  for(int j=0;j<=N;j++) {
    i+=(f(x+dx,y)+f(x,y))*dx/2;
    x+=dx;
  }
  return i;
}


cat Makefile
Код:

ellips: ellips.o
   gcc -s -lstdc++ ellips.o -o ellips

ellips.o: ellips.cxx
   gcc -c ellips.cxx -o ellips.o


Как видите, я увеличил N в два раза и использовал y как отношение a к b в функции f. Думаю это разумное предположение. Если неправ, то извиняюсь.

Вот результаты работы программы:-

cat dat.txt
Код:

0 1.00001 1 1
0.05 1.00486 1.00655 1.01285
0.1 1.016 1.01891 1.02798
0.15 1.03153 1.03509 1.04528
0.2 1.05051 1.0543 1.06465
0.25 1.07231 1.07606 1.08598
0.3 1.09649 1.10001 1.10915
0.35 1.12269 1.12588 1.13406
0.4 1.15066 1.15347 1.16059
0.45 1.18018 1.18258 1.18863
0.5 1.21106 1.21306 1.21807
0.55 1.24317 1.24478 1.24883
0.6 1.27636 1.27763 1.2808
0.65 1.31053 1.31149 1.31389
0.7 1.3456 1.34629 1.34803
0.75 1.38148 1.38195 1.38313
0.8 1.41809 1.41838 1.41913
0.85 1.45538 1.45554 1.45595
0.9 1.4933 1.49336 1.49354
0.95 1.53179 1.5318 1.53184
1 1.5708 1.5708 1.5708
1.05 1.61031 1.61032 1.61036
1.1 1.65028 1.65034 1.6505
1.15 1.69067 1.6908 1.69116
1.2 1.73146 1.73169 1.73231
1.25 1.77261 1.77298 1.77391
1.3 1.81412 1.81463 1.81594
1.35 1.85594 1.85663 1.85836
1.4 1.89808 1.89895 1.90115
1.45 1.9405 1.94158 1.94428
1.5 1.98319 1.98449 1.98774
1.55 2.02613 2.02767 2.03149
1.6 2.06932 2.0711 2.07553
1.65 2.11274 2.11477 2.11984
1.7 2.15637 2.15866 2.16439
1.75 2.20021 2.20277 2.20917
1.8 2.24424 2.24708 2.25418
1.85 2.28845 2.29158 2.29939
1.9 2.33284 2.33626 2.34479
1.95 2.3774 2.38111 2.39038
2 2.42212 2.42612 2.43614
2.05 2.46699 2.47129 2.48207
2.1 2.512 2.5166 2.52815
2.15 2.55716 2.56206 2.57438
2.2 2.60244 2.60765 2.62074
2.25 2.64786 2.65337 2.66724
2.3 2.69339 2.6992 2.71386
2.35 2.73904 2.74516 2.7606
2.4 2.78481 2.79122 2.80745
2.45 2.83067 2.8374 2.85441
2.5 2.87665 2.88367 2.90147
2.55 2.92272 2.93004 2.94862
2.6 2.96888 2.97651 2.99587
2.65 3.01514 3.02306 3.04321
2.7 3.06149 3.0697 3.09063
2.75 3.10792 3.11643 3.13813
2.8 3.15443 3.16323 3.18571
2.85 3.20102 3.21011 3.23336
2.9 3.24768 3.25706 3.28108
2.95 3.29442 3.30409 3.32886
3 3.34123 3.35118 3.37671
3.05 3.38811 3.39834 3.42463
3.1 3.43505 3.44556 3.4726
3.15 3.48206 3.49284 3.52063
3.2 3.52913 3.54018 3.56871
3.25 3.57626 3.58758 3.61684
3.3 3.62344 3.63503 3.66503
3.35 3.67069 3.68253 3.71326
3.4 3.71798 3.73009 3.76154
3.45 3.76533 3.7777 3.80986
3.5 3.81273 3.82535 3.85823
3.55 3.86018 3.87305 3.90663
3.6 3.90768 3.9208 3.95508
3.65 3.95522 3.96859 4.00356
3.7 4.00281 4.01642 4.05208
3.75 4.05044 4.0643 4.10064
3.8 4.09812 4.11221 4.14923
3.85 4.14583 4.16016 4.19786
3.9 4.19359 4.20815 4.24651
3.95 4.24139 4.25617 4.2952
4 4.28922 4.30423 4.34392
4.05 4.33709 4.35233 4.39266
4.1 4.385 4.40045 4.44143
4.15 4.43294 4.44861 4.49024
4.2 4.48091 4.49681 4.53906
4.25 4.52892 4.54503 4.58791
4.3 4.57696 4.59328 4.63679
4.35 4.62503 4.64156 4.68569
4.4 4.67314 4.68987 4.73461
4.45 4.72127 4.7382 4.78356
4.5 4.76943 4.78657 4.83253
4.55 4.81763 4.83495 4.88151
4.6 4.86584 4.88337 4.93052
4.65 4.91409 4.93181 4.97955
4.7 4.96236 4.98027 5.0286
4.75 5.01066 5.02876 5.07766
4.8 5.05898 5.07727 5.12675
4.85 5.10733 5.1258 5.17585
4.9 5.1557 5.17435 5.22497
4.95 5.2041 5.22292 5.2741
5 5.25252 5.27152 5.32325


Длинна, полученная по формуле Салтера различается в сотых от длинны, полученной численным интегрированием.

Это я N увеличил всего в два раза... А если увеличу в три раза - что будет? Но я не буду проверять - может я неверно понял программу. А пока формула Салтера показывает превосходные результаты.

P.S.

Сейчас уже так пишут по ANSI стандарту :-

Код:

#include <iostream>
#include <fstream>
#include <math.h>

Стандартные библиотеки c++ подключаются без расширения .h
Кроме того, все глобальные переменные и константы вынесены в namespace std.
Вернуться к началу
Посмотреть профиль Отправить личное сообщение
Непризнанный Гений
Ветеран


Зарегистрирован: 27.02.2005
Сообщения: 606

СообщениеДобавлено: 15 Июн 2005 [05:48]    Заголовок сообщения: Ответить с цитатой

Anonymous писал(а):

по формуле салтера:6.2830585


Я понял, в чём Вы ошибались.

Это не по формуле Салтера. Это просто 2*Pi.

Так получилось, потому что Вы написали 1/X вместо 1.0/X. Ваш компилятор C++ преобразовал это к (int)1/X и получился 0. А pow(...., 0) дал еденицу.

Заметьте, что я Вашу вторую формулу не правил, но результат получился почти точным, потому что я пользовался gcc. Но по стандарту Вы сами должны константы целые приводить к вещественным, так как тип арифметической операции умножения или деления определяется по типу первого операнда.

то есть целая константа помноженная или поделённая на что-либо даст целое число в Вашем компиляторе.

Что за компилятор? BC 3.0 какой-нибудь?

Я использовал gcc 3.2.2 и GNU/Linux. Считаю, что для численных расчётов это наиболее удачное решение.
Вернуться к началу
Посмотреть профиль Отправить личное сообщение
Гость






СообщениеДобавлено: 15 Июн 2005 [10:01]    Заголовок сообщения: Ответить с цитатой

И ещё ошибочка: for(int j=0;j<=N;j++) надо писать как for(int j=0;j<N;j++)

Sad
Вернуться к началу
Salter
Ветеран


Зарегистрирован: 15.02.2005
Сообщения: 323
Откуда: Соликамск

СообщениеДобавлено: 15 Июн 2005 [19:16]    Заголовок сообщения: Ответить с цитатой

Андрюша писал: "Я вывел из элементарных соображений свою формулу длины периметра эллипса:
L=2*Pi*(a*b+(2*(a-b)/Pi)^2)
Теперь посчитаем численным методом длину дуги эллипса a=1 b=2 по формуле ...
.... ... .... ... ....
А теперь проведем расчет по моей формуле и по формуле салтера:
по моей формуле: 7.4483836
по формуле Салтера: 6.2830585
".

Уважаемый Андрюша, судя по результату Ваших расчетов, мне кажется Вы не научились делать расчет Средней эллиптической (Rср.элл) и длину периметра эллипса (Lэлл)по мой формуле, т.к. при заданных значениях a=1 b=2
у Вас получился ошибочный результат

Если считать по моей формуле, то получится следующее (разжевываю):
Моя формула Rср.элл = ([а^x + b^x]/2)^(1/x)
Также известно, что x = ln2/ln(Pi/2) = 1,53... = 1,5349285356613752020529480451829

Показываю пошаговый расчет
Тогда а^x = 1^1,53... = 1
Тогда b^x = 2^1,53... = 2,9... = 2,8977407650824266294195437386109
Тогда [а^x + b^x] = 1 + 2,9... = 3,9... = 3,8977407650824266294195437386109
Тогда ([а^x + b^x]/2) = 3,9.../2 = 1,9... = 1,9488703825412133147097718693054
Тогда ([а^x + b^x]/2)^(1/x) = (1,9...)^(1/1,53...) = (1,9...)^( 0,65...) = 1,54... = 1,5445162216614120315541811981067

Rср.элл = 1,5445... = 1,5445162216614120315541811981067
Lэлл = 2 * Pi * Rср.элл. = 2 * 3,14... * 1,5445... =
Lэлл = 9,7044816306435134323363799509294
Таким образом, длина периметра эллипса по формуле Салтера будет равна:
Lэлл = 9,7...
а, не 6.2830585 , как утверждаете Вы.

Что интересно, но Вы не смогли даже правильно сделать расчет по своей формуле
L=2*Pi*(a*b+(2*(a-b)/Pi)^2)
где, очевидно, за средний радиус эллипса Вы пытаетесь выдать (О! Ужас!)
Rср.эл. = (a*b+(2*(a-b)/Pi)^2)
т.е. произведение радиусов (a*b) да еще и с каким-то немыслимым плюсом.
Наверняка, Вы тут где-то напутали (или действительно за лохов людей принимаете). Посмотрите, какой расчет получается по Вашей формуле в действительности. Привожу пошаговый расчет:
(2*(a-b)/Pi) = 2*(1-2)/3,14... = 2 * (-1)/3,14... = - 0,6... = 0,63661977236758134307553505349006
(2*(a-b)/Pi)^2 = (- 0,6...)^2 = 0,4... = 0,40528473456935108577551785283891
Rср.эл. = (a*b+(2*(a-b)/Pi)^2) = 1 * 2 + 0,4... = 2,4... =
= 2,4052847345693510857755178528389
L = 2*Pi*(a*b+(2*(a-b)/Pi)^2) = 2 * 3,14...* 2,4... =
= 15,112849703829498326152713747078
Таким образом, длина периметра эллипса по Вашей, уважаемый Андрюша, формуле будет равна:
L = 15,1..
что на много больше величины (7.4483836), на которую ссылались Вы.
Во всем у Вас сквозит торопливость и небрежность в тексте, в результате чего Вы делаете поспешные выводы, кишащие ошибками. Математика этого не любит.
Исходите из того, что я был бы рад найти в своей формуле ошибку. Если бы я этого не хотел, я не стал бы ничего говорить о функции, о ее пределе, о модели геометрической интерпретации перехода к пределу, о секретах вывода формул Салтера. Не расскажи я всего этого, Вы сейчас вообще бились бы головой о стену.
Об явных ошибках в программировании у Вас, уважаемый Андрюша, и у Непризнанного Гения я говорить ничего не буду. Сами разбирайтесь.

Непризнанный Гений, получив величину (5.32325) написал: "Длинна, полученная по формуле Салтера различается в сотых от длинны, полученной численным интегрированием.".

Тут я не понял, на основании чего он делает такой вывод.
Вы обои, уважаемые, если беретесь программировать, то будьте добры указывать обозначение буквенное и цифровое полученных результатов ясно и не двусмысленно, чтобы даже школьник мог понять о каких входных и выходных данных в решении задачи идет речь. Вы должны пояснить откуда Вами взялся входной параметр х = 1.534928536
Я то, величину (х) получил из своей формулы периметра эллипса при переходе к пределу функции при b -> 0

x = ln2/ln(π/2) = 1,53...
= 1,5349285356613752020529480451829
а у вас то, откедова эта величина объявилась?
Тут Вам нужно давать объяснения, раз вы взялись устраивать проверку работы моей формулы.

У меня складывается впечатление, что вами обоими до конца так и не осмыслена еще пока сущность открытия мною новой красивой Средней, которая наилучшим образом решает задачу определения длины периметра эллипса в отличие от всех прочих известных средних величин. Для лучшего уяснения этой сути сравните на примере а = 1; b = 2 попытки определить длину периметра эллипса:

Средний радиус - как Среднее геометрическое:
Rср = (a*b)^0,5 = 2^0,5 = 1,414... ; Lэлл = 2 * Pi * Rср = 2 * 3,14...* 1,414 = 8,886

Средний радиус - как Среднее арифметическое:
Rср = (a+b)/2 = 3/2 = 1,5 ; Lэлл = 2 * Pi * Rср = 2 * 3,14...* 1,5 = 9,4248

Средний радиус - как Среднее квадратичное:
Rср = [(a^2 + b^2)/2]^0,5 = [(1 + 4)/2]^0,5=2,5^0,5 = 1,581... ;
Lэлл = 2 * Pi * Rср = 2 * 3,14...* 1,581...= 9,935

Средний радиус - как Среднее эллиптическое (по Салтеру):
Rср = Rср.элл = ([а^1,53.. + b^1,53..]/2)^(1/1,53..)= 1,54...;
Lэлл = 2 * Pi * Rср = 2 * 3,14...* 1,54... = 9,7044

Вот теперь и пробуйте дать более точный расчет и среднего радиуса эллипса - Rср.элл. (т.е. радиуса такого круга, длина окружности которого равнялась бы длине периметра данного эллипса) и соответствующей этому радиусу - длины периметра эллипса - Lэлл. А если не можете, то и не беритесь.
В моих расчетах все предельно ясно. А в ваших - только понапрасну время терять, сплошной бардак.
Уважаемый Непризнанный Гений, я Вам очень благодарен за поддержку моей позиции и моей формулы. Простите, что мне пришлось покритиковать и Вас. Если в чем я тут ошибся, то поправьте меня. Я был бы рад этому.

С уважением, Salter.
Вернуться к началу
Посмотреть профиль Отправить личное сообщение Отправить e-mail Посетить сайт автора
Гость






СообщениеДобавлено: 15 Июн 2005 [20:31]    Заголовок сообщения: Ответить с цитатой

Слышь сальтер на все вопросы по порядку:
1. Я в моей формуле забыл одну штуку дописать-мог заметить по размерности. На самом деле я забыл там корень дописать.Формула выглядит так: L=2*Pi*sqrt(a*b+(2*(a-b)/Pi)^2)
2.Программа считает одновременно и интеграл, и твою ф-лу, и мою ф-лу-поэтому на входе получается этот 1.53..

Зацени пока точность моей формулы-тоже неплохо.

Завтра закончу исследование и напишу окончательный результат.
Вернуться к началу
Непризнанный Гений
Ветеран


Зарегистрирован: 27.02.2005
Сообщения: 606

СообщениеДобавлено: 16 Июн 2005 [03:33]    Заголовок сообщения: Ответить с цитатой

2 Salter:-

Это программа Андрюшы. Я лишь привёл её в работоспособный вид. Я не исправляю всех ошибок в чужих программах - только критические. Предпочитаю делать всё по минимуму с чужими программами. В программировании иначе нельзя.

В файле dat.txt, генерируемом программой, есть три столбца:
<b> <результат численного интегрирования> <результат по формуле Салтера> <результат по формуле Андрюшы>.

Максимальная ошибка достигается для a=1 b=5.
Код:

5 5.25252 5.27152 5.32325


Для случая a=1 и b=2 программа даёт такой результат.
Код:

2 2.42212 2.42612 2.43614


Собственно, это меня и удивило - по программе Андрея получалось совсем не то, о чём он писал. А всё из-за 1/X вместо 1.0/X Cool
Вернуться к началу
Посмотреть профиль Отправить личное сообщение
Salter
Ветеран


Зарегистрирован: 15.02.2005
Сообщения: 323
Откуда: Соликамск

СообщениеДобавлено: 16 Июн 2005 [04:42]    Заголовок сообщения: Ответить с цитатой

Спасибо уважаемый Непризнанный Гений за разъяснение. Я так и понял.

to Андрюша.
Уважаемый Андрюша,
точность Вашей формулы по определению Среднего радиуса эллипса и его периметра, действительно неплохая.
Теперь нужно проверить ее (формулы) поведение, как функции при устремлении (a/b) к бесконечности, т.е. при b -> 0 , и тем самым уточнить правомерность применения принятого Вами показателя степени (у Вас этот показатель пока равен числу 2).
Для этого смело берите предел функции (х) при b -> 0

Lim {2*Pi*[(a*b+(2*(a-b)/Pi)^x)]^(1/x)} = 4 * a
b -> 0


Найдете чему равен х , тогда человечество получит, действительно, еще одну (достаточно точную) формулу для расчета длины периметра эллипса и его среднего радиуса (если, конечно, Вы умеете решать уравнения с Пределами). Smile
Может быть тогда, вы сможете оценить роль и значение мозговой атаки на космическом форуме Kuasara по достоинству и примите в ней активное участие
Вернуться к началу
Посмотреть профиль Отправить личное сообщение Отправить e-mail Посетить сайт автора
Непризнанный Гений
Ветеран


Зарегистрирован: 27.02.2005
Сообщения: 606

СообщениеДобавлено: 16 Июн 2005 [05:02]    Заголовок сообщения: Ответить с цитатой

А можно ещё x по-другому искать - численно.

Пусть f(a,b,x) - это наша параметрическая формула длинны эллипса.

1. Задать набор точек (a,b) с весовыми коэффициентами r(a, b)
2. Просчитать для них точную длинну эллипса l(a, b)
3. Найти x при котором следующая сумма даёт минимум:

Код:

   Sum(по всем (a,b); r(a,b)*abs(f(a,b,x) - l(a,b)))


Тогда уравнения с пределами не надо решать - компьютер сам всё сделает.
Вернуться к началу
Посмотреть профиль Отправить личное сообщение
Salter
Ветеран


Зарегистрирован: 15.02.2005
Сообщения: 323
Откуда: Соликамск

СообщениеДобавлено: 16 Июн 2005 [07:04]    Заголовок сообщения: Ответить с цитатой

Давайте не будем, уважаемый Непризнанный Гений, пренебрегать аналитикой и уходить в сторону от вопроса.
Возьмите сначала тот Предел функции, который я показал, а потом решим, что нужно делать (как поступить) с формулой Андрюши.

С уважением, Salter.
Вернуться к началу
Посмотреть профиль Отправить личное сообщение Отправить e-mail Посетить сайт автора
Непризнанный Гений
Ветеран


Зарегистрирован: 27.02.2005
Сообщения: 606

СообщениеДобавлено: 16 Июн 2005 [09:48]    Заголовок сообщения: Ответить с цитатой

Так он не зависит от x.
Вернуться к началу
Посмотреть профиль Отправить личное сообщение
Salter
Ветеран


Зарегистрирован: 15.02.2005
Сообщения: 323
Откуда: Соликамск

СообщениеДобавлено: 16 Июн 2005 [16:09]    Заголовок сообщения: Ответить с цитатой

Раз предел при любом значении х один и тот же, то это означает только одно - функция в формуле Андрюши образует расходящийся ряд, т.е. при любом значении Х (в т.ч. и при х = 2) мы будем получать все что угодно, только не длину периметра эллипса. В противном случае длина периметра эллипса по формуле Андрюши у конкретной пары чисел а и b при х = 2 была бы точно такой же, как и при х = 1 , как и при х = 3 и т.д.
Поэтому число 2 столь же случайно, как и числа 1,52 или 4,59, дающие при некоторых значениях a/b достаточно близкие значения длины периметра эллипса.
Поэтому пользоваться формулой Андрюши я бы никому не рекомендовал.
Факир был пьян - фокус не удался!
В отличие от подобного рода "формул" Средняя эллиптическая (по Салтеру) переживет века. Smile


С уважением, Salter.
Вернуться к началу
Посмотреть профиль Отправить личное сообщение Отправить e-mail Посетить сайт автора
Гость






СообщениеДобавлено: 16 Июн 2005 [19:17]    Заголовок сообщения: Ответить с цитатой

Численное исследование закончено.
Разбиение при взятии интеграла-1000000 точек.
Число пи и число сальтера-точность-все разряды калькулятора.

Посчитал интеграл, число сальтера и по моей формуле
Построил график погрешностей(абсолютных)
И что же получилось?
Не спорю, что формула Сальтера хорошая аппроксимация, но не более того. Расхождение в четвертом знаке стабильно, тогда как
при уменьшении точности числа пи на 10(!) знаков(впрочем как у числа х). Так что как нам утверждал Сальтер-погрешность его формулы именно в том, что она аппроксимация, а не из-за числа пи или х. Я уже несколько раз писал что если бы формула была абсолютно точна, то она была бы решением эллитического интеграла, но сальтер это игнорирует.
Также при уменьшении числа точек разбиения точность интеграла не меняется-ставил опыт с прогой.
Итог: Формула сальтера-хорошая аппроксимация, она не является точной. Расхождение стабильно в четвертом знаке.
А моя формула тоже неплохая и попроще-максимум погрешности 2%, для физики пойдет на ура, так как гораздо проще сальтеровской.

Если сальтер мне не верит на слово, пусть напишет на маил tahtamouse1@mail.ru письмо со своим маилом, я ему скину прогу и построенные мною графики.

С уважением Андрюша Мехедов
Вернуться к началу
Гость






СообщениеДобавлено: 17 Июн 2005 [04:02]    Заголовок сообщения: Ответить с цитатой

2 Андрюша

А какова максимальная и минимальная погрешность формулы Салтера по Вашей программе?
Вернуться к началу
zenixt
Ветеран


Зарегистрирован: 14.04.2005
Сообщения: 1278
Откуда: Россия

СообщениеДобавлено: 19 Июн 2005 [11:18]    Заголовок сообщения: Ответить с цитатой

Андрюша говорил: "рассмотрим формулу Q=I^2*R*t
Если мы сейчас сюды вместо t подставим -t-то будет эквивалентно, если заменить Q на -Q или R на -R с точки зрения эксперимента мы знаем что невозможно подводя к проводу тепло получить акой уровень напряжения-конечно некоторая разность потенциалов появится, но она будет обусловлена тепловым шумом электронов Джонсона, а про материалы, которые вырабатывали ток и были по структуре такие же как обычные резисторы я не слышал(возможно слышали вы).

Так вот такое вот доказательство.

Я конечно догадываюсь чего вы хотите добиться этим эффектом, и наверное добились бы, но сначала преодолейте те трудности которые перед вами поставил хотя бы я."
Андрюша. Много вы от меня хотите. Человеческие возможности ограничены. Должна быть специализация. А я решаю только те задачи, которые МНЕ интересны.
Вернуться к началу
Посмотреть профиль Отправить личное сообщение Отправить e-mail Посетить сайт автора
Андрюша
Новичок


Зарегистрирован: 08.06.2005
Сообщения: 44
Откуда: долгопа

СообщениеДобавлено: 23 Июн 2005 [13:33]    Заголовок сообщения: Ответить с цитатой

сальтер, ты куда пропаал?
Вернуться к началу
Посмотреть профиль Отправить личное сообщение Отправить e-mail Посетить сайт автора
Андрюша
Новичок


Зарегистрирован: 08.06.2005
Сообщения: 44
Откуда: долгопа

СообщениеДобавлено: 23 Июн 2005 [13:45]    Заголовок сообщения: Ответить с цитатой

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

dS/dt>0 для замкнутых систем! вот в чем сиила, а не в том, чтобы время повернуть в обратную сторону. И притом с макросистемами, даже если ты время повернещь в обратную сторону, ты обратного процесса не получишь-как например если ты смешиваешь диффузионно два газа, то если время назад пустить, то они не разделятся все равно. Или в крайнем случае(пускай даже у тебя в распоряжении есть работающий причинно-следственный принцип) если газы разделятся, потом они снова смешаются!

то есть есть подкласс необратимых процессов, которые даже поворотом времени назад не пустишь в обратную сторону.
Кстатит к ним относится горение-при горении излучаетсяв пространство информация. Если ты смеешь рассуждать про микромир,
то должен знать, что все элементарные частицы дуальны. И вот как раз эту дуальность компенсирует информация. Поступает к частице информация-допустим мы на нее смотрим-и она превращается в частицу. А когда не смотрим-она превращается в волну-электрон может интерферировать сам с собой! Так что информация объективно существует, и в процессе горения она излучается. Кстати потрясная связь энтропии и информации, а получена она кстати из формальной термодинамики и статистики, как раз не случайно, а слуит для компенсации дуализма.

И кстати ты велел мне найти ошибку в твоем расчете ТБШ. Я нашел и написал, но ты че то не реагируешь.
Если ты
Вернуться к началу
Посмотреть профиль Отправить личное сообщение Отправить e-mail Посетить сайт автора
Salter
Ветеран


Зарегистрирован: 15.02.2005
Сообщения: 323
Откуда: Соликамск

СообщениеДобавлено: 24 Июн 2005 [01:51]    Заголовок сообщения: Ответить с цитатой

Андрюша писал(а):
сальтер, ты куда пропаал?


Я здеся. Very Happy
А что такое?
Вернуться к началу
Посмотреть профиль Отправить личное сообщение Отправить e-mail Посетить сайт автора
Показать сообщения:   
Начать новую тему   Ответить на тему    Список форумов Форум космопорта -> Поговорим? Часовой пояс: GMT + 2
На страницу Пред.  1, 2, 3  След.
Страница 2 из 3

 
Перейти:  
Вы не можете начинать темы
Вы не можете отвечать на сообщения
Вы не можете редактировать свои сообщения
Вы не можете удалять свои сообщения
Вы не можете голосовать в опросах