RFpro.ru: Программирование на C / C++ (comp.soft.prog.cplus) : Рассылка : Subscribe.Ru

Подписаться Бесплатная «Золотая» новостная рассылка . Подписчиков 5.420 RSS

← Октябрь 2009 →
	1	2	3	4
5	6	7	8	9	10	11
12	13	14	15	16	17	18
19	20	21	22	23	24	25
26	27	28	29	30	31

За последние 60 дней ни разу не выходила

Сайт рассылки: http://prog.rusfaq.ru/C
Открыта: 14-03-2002

Автор

Калашников О.А.

Статистика

5.420 подписчиков
0 за неделю

← Все выпуски →

RFpro.ru: Программирование на C / C++

Хостинг портала RFpro.ru:
Московский хостер
Профессиональный платный хостинг на базе Windows 2008

РАССЫЛКИ ПОРТАЛА RFPRO.RU

Чемпионы рейтинга экспертов в этой рассылке

_Ayl_
Статус: Студент
Рейтинг: 1366
∙ повысить рейтинг »

Micren
Статус: Бакалавр
Рейтинг: 910
∙ повысить рейтинг »

Boriss
Статус: Академик
Рейтинг: 861
∙ повысить рейтинг »

/ КОМПЬЮТЕРЫ И ПО / Языки программирования / C/C++

Номер выпуска:	1402
Дата выхода:	11.10.2009, 21:30
Администратор рассылки:	Dr_Andrew, Модератор
Подписчиков / экспертов:	639 / 166
Вопросов / ответов:	1 / 1
IRC-канал по теме:	#C

Вопрос № 172943: Здравствуйте!Помогите реализовать с помощью интерполяционного многочлена лагранжа! Уплотнение в 2 раза таблицы функции ,заданной на равномерной сетке с оценкой погрешности результата. Все входные данные размещаются во входном текстовом файле,о...

Вопрос № 172943:

Здравствуйте!Помогите реализовать с помощью интерполяционного многочлена лагранжа!
Уплотнение в 2 раза таблицы функции ,заданной на равномерной сетке с оценкой погрешности результата.
Все входные данные размещаются во входном текстовом файле,он содержит:
Y-значения функции,расположенные по возрастанию аргумента.
N-число точек,в которых задана функция.
Выходной файл содержит:
eps- вычисленная погрешность результата
YY-уплотнёный вектор
I- код ошибки,равный 0,если ошибок нет и 1 если N<5.

Отправлен: 05.10.2009, 22:27
Вопрос задал: РАИ, Посетитель
Всего ответов: 1
Страница вопроса »

Отвечает Лысков Игорь Витальевич, Модератор :
Здравствуйте, РАИ.
Программа в приложении. Выводим и в файл, и на экран
Во входном файле вводим:
"Число точек n"
"x0" "шаг h"
"таблица из n+1 значений функции"
Примерный входной файл in.txt:

Код:

10
1,0 ,5
1, 2, 2,5 2,6 2,65
2,7 2,9 3,6 4,9 5,8 7,5

Выходной файл out.txt ("yy[i]" "погрешность Eps", в конце "код ошибки"):

Код:

1,000000 0,000000
1,442023 0,000323
2,000000 0,000000
2,342400 0,000051
2,500000 0,000000
2,565463
  0,000015
2,600000 0,000000
2,626934 0,000007
2,650000 0,000000
2,670858 0,000005
2,700000 0,000000
0

Приложение:
#include <locale> #include <iostream> #include <stdio.h> #include <math.h> using namespace std; //Функция вводит данные из файла //Параметры: //file - имя файла //n - адрес количества точек //x0 - адрес начального значения x0 //h - адрес шага //y - адрес массива табличных значений функции int GetData(char * file, int * n, float * x0, float * h, float ** y) { FILE *fin = fopen( file, "r" ); //открываем файл на чтение int i; if( fin == NULL ) //проверка на наличие файла { wcout<<L"Файл "; cout<<file; wcout<<L" не найден"<<endl; return 0; } else { //вводим данные, одновременно wcout<<L"Исходные данные:"<<endl;//выводим на экран cout<<"n = "; if (fscanf(fin, "%d", n )) //количество { cout<<*n; *y = (float*)malloc((*n+1)*sizeof(float)); //запросим память } else { wcout<<L"ошибка ввода"<<endl; return 0; //при ошибке прекращаем ввод } cout<<" x0 = "; if (fscanf(fin, "%f", x0 )) //x0 cout<<*x0; else { wcout<<L"ошибка ввода"<<endl; return 0; } cout<<" h = "; if (fscanf(fin, "%f", h )) //шаг cout<<*h<<endl; else { wcout<<L"ошибка ввода"<<endl; return 0; } for(i=0; i<=*n ; i++) //вводим n+1 значений функции { cout<<"y["<<i<<"] = "; //одновременно выводим if (fscanf( fin, "%f", &(*y)[i] )) cout<<(*y)[i]<<endl; else { wcout<<L"ошибка ввода"<<endl; return 0; } } fclose( fin ); //закрываем файл return 1; //все ок } } //Функция выводит данные в файл //Параметры: //file - имя файла //n - адрес количества точек //yy - адрес массива интерполированных значений функции int SetData(char * file, int n, float * yy, float * eps) { FILE *fout = fopen( file, "w+" );//открываем файл на запись int i; if( fout == NULL ) //проверяем на ошибку создания { wcout<<L"Файл "; cout<<file; wcout<<L" не создан"<<endl; return 0; } else { wcout<<L"Результат:"<<endl; //выводим на экран и в файл for (i=0; i<=n; i++) //n+1 значений { //в одной строке элемент yy и погрешность fprintf(fout, "%f %f\n", yy[i], eps[i]); cout<<"yy["<<i<<"] = "<<yy[i]<<"\teps = "<<eps[i]<<endl; } if (n < 5) //проверим условие n < 5 { fprintf(fout, "1\n"); wcout<<L"Код ошибки: 1"<<endl; } else { fprintf(fout, "0\n"); wcout<<L"Код ошибки: 0"<<endl; } fclose( fout ); return 1; } } //расчет Li(x) = П(j=0,n;j!=i)[(x-x0-j*h)/(h*(i-j)] //Параметры: //x - аргумент, для которого расчитываем интерполяцию //i - индекс коэффициента Li(x) //n - число точек //x0- начальное значение аргумента //h - шаг float CalcL(float x, int i, int n, float x0, float h) { float l; //тут будем накапливать произведение int j; for(l=1,j=0; j<=n; j++) { if (j != i) { l *= x - x0 - j * h; l /= h * (i - j); } } return l; } //расчет F(x) = Сум(i=0,n)[y(i)*Li(x)] //Параметры: //x - аргумент, для которого расчитываем интерполяцию //y - адрес массива табличных значений функции //n - число точек //x0- начальное значение аргумента //h - шаг float CalcF(float x, float * y, int n, float x0, float h) { float f; //тут будем накапливать сумму int i; for(f=0,i=0; i<=n; i++) f += y[i] * CalcL(x, i, n, x0, h); return f; } //расчет погрешности Eps = ((П(i=0,n)[x-x0-i*h])*max(y))/n! //Параметры: //x - аргумент, для которого расчитываем интерполяцию //y - адрес массива табличных значений функции //n - число точек //x0- начальное значение аргумента //h - шаг float CalcEps(float x, float * y, int n, float x0, float h) { float eps, M, fact; int i; for(M=fabs(y[0]),fact=1.,eps=1.,i=0; i<n+1; i++) { eps *= x - x0 - i * h; if (M < fabs(y[i])) M = fabs(y[i]); if (i > 1) fact *= i; } return fabs(eps * M / fact); } int main() { int i,n; float x0; float h; float *y = 0; float *yy = 0; float *eps = 0; locale::global(locale("russian_russia.866")); //чтобы писалось по-русски if (GetData("in.txt", &n, &x0, &h, &y)) { //запросим память yy = (float*)malloc((2*n+1)*sizeof(float)); eps = (float*)malloc((2*n+1)*sizeof(float)); //считаем результат интерполяции и погрешность for(i=0; i<21; i++) { yy[i] = CalcF(x0+i*h/2, y, n, x0, h); eps[i] = CalcEps(x0+i*h/2, y, n, x0, h); } //выводим результат SetData("out.txt", n, yy, eps); //свободим память free(eps); free(yy); free(y); return 0; } if (y) //может быть выделена free(y); return 1; }
-----
Удачи!

Ответ отправил: Лысков Игорь Витальевич, Модератор
Ответ отправлен: 10.10.2009, 06:16
Украина, Кировоград
Тел.: +380957525051
ICQ # 234137952
Mail.ru-агент: igorlyskov@mail.ru
Абонент Skype: igorlyskov

Как сказать этому эксперту "спасибо"?

Отправить SMS #thank 255244 на номер 1151 (Россия) | Еще номера »

Отправить WebMoney:

Вам помогли? Пожалуйста, поблагодарите эксперта за это!

Оценить выпуск »
Нам очень важно Ваше мнение об этом выпуске рассылки!

Задать вопрос экспертам этой рассылки »

Скажите "спасибо" эксперту, который помог Вам!

Отправьте СМС-сообщение с тестом #thank НОМЕР_ОТВЕТА
на короткий номер 1151 (Россия)
Номер ответа и конкретный текст СМС указан внизу каждого ответа.

Полный список номеров »

* Стоимость одного СМС-сообщения от 7.15 руб. и зависит от оператора сотовой связи. (полный список тарифов)
** При ошибочном вводе номера ответа или текста #thank услуга считается оказанной, денежные средства не возвращаются.
*** Сумма выплаты эксперту-автору ответа расчитывается из суммы перечислений на портал от биллинговой компании.

© 2001-2009, Портал RFpro.ru, Россия
Авторское право: ООО "Мастер-Эксперт Про"
Автор: Калашников О.А. | Программирование: Гладенюк А.Г.
Хостинг: Компания "Московский хостер"
Версия системы: 2009.6.9 от 25.09.2009

В избранное

{#template MAIN} <div id="loginForm" style="display:none;" class="subscriberu_popup"> <div class="popup_register"> {#include js_tmpl_auth_reg_tab} {#if $P.login_register_tab == 1} <form class="authentication-form" method="post" action="/MEMBERLOGIN_authen_cred"> <dl class="rg_block_options"> <dt id="js_tap_panel_auth"> <h1>Войти на сайт</h1> {* {#include js_tmpl_auth_reg_button} *} {#include js_tmpl_auth_reg_action} <hr class="logreg_line noPhones"> <div class="logreg_descr noPhones"><p>{#include js_tmpl_auth_reg_descr} </p></div> <div class="logreg_advice noPhones"> Если вы еще не с нами, то начните с <a href="#" onclick="rgNav('js_tab_reg');return false;" class="dashed" data-func="registr">регистрации</a> </div> <br><br> <a class="dashed auth-enter" href="/manage/author/"><b>Вход для авторов</b></a> </dt> </dl> </form> {#/if} {#if $P.login_register_tab == 2} <div class="rg_block_options"> <div id="js_tap_panel_auth"> <h1>Регистрация</h1> <div class="social_reg"> {* <div class="rg_description">{#include js_tmpl_soc_auth_reg_descr}</div> *} {#include js_tmpl_auth_reg_soc} <div class="rg_soc_auth_agree">{#include js_tmpl_auth_reg_agree}</div> </div> <div class="subscribe_reg"> {* <div class="rg_description"> #include js_tmpl_auth_reg_descr </div> *} {#include js_tmpl_auth_reg_action} </div> {* {#include js_tmpl_auth_reg_button} *} <div class="clr"> </div> <hr class="logreg_line noPhones"> <div class="logreg_descr noPhones">{#include js_tmpl_auth_reg_descr} {#include js_tmpl_soc_auth_reg_descr} </div> </div> </div> {#/if} </div> {* <div class="gray_bg register_shadow"></div> *} </div> {#/template MAIN} {#template js_tmpl_auth_reg_tab} <ul class="rg_nav"> <li id="js_tab_auth" class="{#if $P.login_register_tab == 1} rg_active_nav {#/if} rg_first_nav"><a onclick="rgNav('js_tab_auth');return false;" href="">Вход на сайт</a></li> <li id="js_tab_reg" class="{#if $P.login_register_tab == 2} rg_active_nav {#/if}"><a onclick="rgNav('js_tab_reg');return false;" href="">Регистрация </a></li> </ul> <span onclick="hidebo();" class="rg_closed"> </span> {#/template js_tmpl_auth_reg_tab} {#template js_tmpl_auth_reg_action} {#if $P.login_register_tab == 1} {#include js_tmpl_auth_reg_soc} {#/if} <div class="rg_forms"> <input type="hidden" id="login_register_destination" value="{$P.login_register_destination}"/> {#if $P.login_register_tab == 1} <div class="rg_for_input"> <span class="rg_text_inner">E-mail или код подписчика</span> <input id="credential_0" class="js_keydown_selector rg_input_text" data-js_submit="no" data-js_next_input_name="credential_1" name="" type="text" /> </div> <div class="rg_for_input"> <span class="rg_text_inner">Пароль</span> <input id="credential_1" class="js_keydown_selector rg_input_text" data-js_submit="yes" data-js_action="js_loginFormBut" name="" type="password" onkeyup="showAttention(this,!!window.event.shiftKey)" /> <span class="pswd_attention" id="attention_pswd"> <span class="icon_attention"></span> <span class="pswd_attention-text" id="attention-text_pswd1">Русская раскладка клавиатуры!</span> <span class="pswd_attention-text" id="attention-text_pswd2">У вас включен Caps Lock!</span> <span class="pswd_attention-text" id="attention-text_pswd3">У вас включен Caps Lock и русская раскладка клавиатуры!</span> </span> </div> <div class="rg_for_input input-alien"> <span class="chk noPhones"><input id="chk_alien" name="" type="checkbox" /></span><label for="chk_alien" class="noPhones"> Чужой компьютер</label> <a class="forgot_pass" href="/member/totalrecall">Забыли пароль?</a> </div> <div class="rg_for_input"> <em id="auth_msg" class="reg_error"></em> <input id="lf_typeauthid" value="email" type="hidden"> <input type="submit" class="button button-red logreg_submit" id="js_loginFormBut" value="Войти">  <div class="loading loading-cover" style="display: none;"><div class="loader"></div></div> </div> {#/if} {#if $P.login_register_tab == 2} <div class="rg_for_input"> <span class="rg_text_inner">E-mail</span> <input id="arfemail" class="js_keydown_selector rg_input_text" name="" type="text" data-js_submit="yes" data-js_action="js_regFormBut"/> </div> <div class="rg_for_input rg_set_lineh rg_for_input_wide"> <label class="js_tap_panel_checkbox"> <span class="chk"><input name="" id='js_tap_panel_checkbox_terms' type="checkbox" data-js_submit="yes" /></span> Я ознакомился и согласен с <a class="link_txd logreg_accLink" href="/faq/vereinbarung.html">условиями сервиса Subscribe.ru</a> </label> <br /> <label class="js_tap_panel_checkbox"> <span class="chk"><input name="" id='js_tap_panel_checkbox_personal' type="checkbox" data-js_submit="yes" /></span> Нажимая на кнопку "Готово!", я даю <a class="link_txd logreg_accLink" href="/faq/persverordnung.html">согласие на обработку персональных данных</a> </label> </div> {* <div style="float: left;position: absolute;left: 11em;"> <img src="http://www.kupivip.ru/images/vip/logo.png?1604" style="width: 86px; vertical-align: middle;display: block;"> </div> <div class="rg_for_input rg_set_lineh"> <label class="js_tap_panel_checkbox"><input name="" id="js_tap_panel_checkbox_kupivip" type="checkbox" data-js_submit="yes"> Я хочу получать новости о скидках на одежду</label> </div> *} <div class="rg_for_input"> <em id="reg_msg" class="reg_error rg_for_input_wide"></em> <em id="reg_msg2" class="reg_error rg_for_input_wide"></em> <input id="rf_typeauthid" value="email" type="hidden"> <a class="button button-red logreg_submit" id="js_regFormBut" href="#">Готово!</a> <div class="loading loading-cover" style="display: none;"><div class="loader"></div></div> </div> {#/if} </div> {#/template js_tmpl_auth_reg_action} {#template js_tmpl_auth_reg_agree} <div class="rg_for_input rg_set_lineh rg_for_input_wide"> <label class="js_tap_panel_checkbox"> <span class="chk"><input name="" id='js_tap_panel_checkbox_terms_reg' type="checkbox" data-js_submit="yes" /></span> Я ознакомился и согласен с <a class="link_txd logreg_accLink" href="/faq/vereinbarung.html">условиями сервиса Subscribe.ru</a></label> <em id="reg_msg_soc" class="reg_error rg_for_input_wide"></em> </div> {#/template js_tmpl_auth_reg_agree} {#template js_tmpl_auth_reg_button} <div class="rg_butons_socials"> {#if $P.login_register_tab == 1} <a class="rg_btn_soc rg_bs_01 js_tap_panel_selector" action="auth_email" href="#"><span><i></i>Email</span></a> <a class="rg_btn_soc rg_bs_01 js_tap_panel_selector" action="auth_openid" href="#"><span><i></i>OpenID</span></a> <a class="rg_btn_soc rg_bs_02 js_tap_panel_selector" action="auth_vkontakte" href="#"><span><i></i>Вконтакте</span></a> <a class="rg_btn_soc rg_bs_02 js_tap_panel_selector" action="auth_mailru" href="#"><span><i></i>Mail.Ru</span></a> {#/if} {#if $P.login_register_tab == 2} <a class="rg_btn_soc rg_bs_01 js_tap_panel_selector" action="reg_email" href="#"><span><i></i>Email</span></a> <a class="rg_btn_soc rg_bs_01 js_tap_panel_selector" action="reg_openid" href="#"><span><i></i>OpenID</span></a> <a class="rg_btn_soc rg_bs_02 js_tap_panel_selector" action="reg_vkontakte" href="#"><span><i></i>Вконтакте</span></a> <a class="rg_btn_soc rg_bs_02 js_tap_panel_selector" action="reg_mailru" href="#"><span><i></i>Mail.Ru</span></a> {#/if} </div> {#/template js_tmpl_auth_reg_button} {#template js_tmpl_auth_reg_descr} {#if $P.login_register_tab == 1} Для оформления подписки на выбранную рассылку, работы с интересующей вас группой или доступа в нужный вам раздел, просим авторизоваться на Subscribe.ru {#/if} {#if $P.login_register_tab == 2} Для регистрации укажите ваш e-mail адрес. Адрес должен быть действующим, на него сразу после регистрации будет отправлено письмо с инструкциями и кодом подтверждения. {#/if} {#/template js_tmpl_auth_reg_descr} {#template js_tmpl_soc_auth_reg_descr} Или зарегистрируйтесь через социальную сеть. {#/template js_tmpl_soc_auth_reg_descr} {#template js_tmpl_auth_reg_soc} <div class="rg_soc"> {#if $P.login_register_tab == 1} <a onclick="return _checkSocConfirm(event)" href="https://oauth.vk.com/authorize?client_id=3954260&scope=wall,offline,photos,groups,video,audio,email&redirect_uri={location.protocol+'//'+location.host}/member/login/vk/&response_type=code&v=5.15" class="login_register_vk_button"> <span class="login_register_vk_icon"></span> </a> {#/if} {#if $P.login_register_tab == 2} <a onclick="return _checkSocConfirm(event)" href="https://oauth.vk.com/authorize?client_id=3954260&scope=wall,offline,photos,groups,video,audio,email&redirect_uri={location.protocol+'//'+location.host}/member/join/vk&response_type=code&v=5.15" class="login_register_vk_button"> <span class="login_register_vk_icon"></span> </a> {#/if} </div> {#/template js_tmpl_auth_reg_soc}

{#template MAIN} <div id="loginForm" style="display:none;" class="subscriberu_popup"> <div class="popup_register"> {#include js_tmpl_auth_reg_tab} <dl class="rg_block_options"> <dt id="js_tap_panel_auth"> <p class="rg_description">{#include js_tmpl_auth_reg_descr}</p> <div class="clr"> </div> {#include js_tmpl_auth_reg_action} <div class="clr"> </div> </dt> </dl> </div>  </div> {#/template MAIN} {#template js_tmpl_auth_reg_tab} <ul class="rg_nav"> <li id="js_tab_reg" class="rg_active_nav rg_first_nav"><a href="" onclick="return false;" >Регистрация</a></li> </ul> <span onclick="hidebo();" class="rg_closed"> </span> {#/template js_tmpl_auth_reg_tab} {#template js_tmpl_auth_reg_descr} <strong>Пожалуйста, подтвердите ваш адрес.</strong><br><br>Вам отправлено письмо для подтверждения вашего адреса {$P.register_confirm_mail}.<br>Для подтверждения адреса перейдите по ссылке из этого письма. {#/template js_tmpl_auth_reg_descr} {#template js_tmpl_auth_reg_action} <div class="rg_forms confirm_code_from_letter"> <div class="rg_for_input"> <span class="rg_inp_descr" style="width:15em;">Или введите код из письма:</span> <input type="text" value="" id="confirm_code" name="" data-js_submit="yes" data-js_action="js_confirmFormBut" class="js_keydown_selector rg_input_text_conf" > </div> <div class="rg_for_input"><label>Не пришло письмо? <b>Пожалуйста, проверьте папку Спам</b><br /> (папку для нежелательной почты).</label><br /> <a href="" onclick="ajax_recall_code();return false" >Вышлите мне письмо еще раз!</a></div> <div class="rg_for_input"> <em class="reg_error" id="confirm_msg"></em> <a href="#" class="button button-red" id="js_confirmFormBut">Готово</a> <div class="loading loading-cover" style="display: none;"><div class="loader"></div></div> <br> </div> </div> {#/template js_tmpl_auth_reg_action}