Обчислює арккосинус, використовуючи метод CORDIC. Практична робота з Програмування. Робота № 500255

Обчислює арккосинус, використовуючи метод CORDIC

Інформація про навчальний заклад

ВУЗ:

Національний університет Львівська політехніка

Інститут:

Не вказано

Факультет:

Не вказано

Кафедра:

Не вказано

Інформація про роботу

Рік:

2015

Тип роботи:

Практична робота

Предмет:

Програмування

Завантажити

Частина тексту файла (без зображень, графіків і формул):

АРККОСИНУС function theta = arccos_cordic ( t, n ) %*****************************************************************************80 % %% ARCCOS_CORDIC returns the arccosine of an angle using the CORDIC method. % % Licensing: % % This code is distributed under the GNU LGPL license. % % Modified: % % 15 June 2007 % % Author: % % John Burkardt % % Reference: % % Jean-Michel Muller, % Elementary Functions: Algorithms and Implementation, % Second Edition, % Birkhaeuser, 2006, % ISBN13: 978-0-8176-4372-0, % LC: QA331.M866. % % Parameters: % % Input, real T, the cosine of an angle. -1 <= T <= 1. % % Input, integer N, the number of iterations to take. % A value of 10 is low. Good accuracy is achieved with 20 or more % iterations. % % Output, real THETA, an angle whose cosine is T. % % Local Parameters: % % Local, real ANGLES(60) = arctan ( (1/2)^(0:59) ); % angles = [ ... 7.8539816339744830962E-01, ... 4.6364760900080611621E-01, ... 2.4497866312686415417E-01, ... 1.2435499454676143503E-01, ... 6.2418809995957348474E-02, ... 3.1239833430268276254E-02, ... 1.5623728620476830803E-02, ... 7.8123410601011112965E-03, ... 3.9062301319669718276E-03, ... 1.9531225164788186851E-03, ... 9.7656218955931943040E-04, ... 4.8828121119489827547E-04, ... 2.4414062014936176402E-04, ... 1.2207031189367020424E-04, ... 6.1035156174208775022E-05, ... 3.0517578115526096862E-05, ... 1.5258789061315762107E-05, ... 7.6293945311019702634E-06, ... 3.8146972656064962829E-06, ... 1.9073486328101870354E-06, ... 9.5367431640596087942E-07, ... 4.7683715820308885993E-07, ... 2.3841857910155798249E-07, ... 1.1920928955078068531E-07, ... 5.9604644775390554414E-08, ... 2.9802322387695303677E-08, ... 1.4901161193847655147E-08, ... 7.4505805969238279871E-09, ... 3.7252902984619140453E-09, ... 1.8626451492309570291E-09, ... 9.3132257461547851536E-10, ... 4.6566128730773925778E-10, ... 2.3283064365386962890E-10, ... 1.1641532182693481445E-10, ... 5.8207660913467407226E-11, ... 2.9103830456733703613E-11, ... 1.4551915228366851807E-11, ... 7.2759576141834259033E-12, ... 3.6379788070917129517E-12, ... 1.8189894035458564758E-12, ... 9.0949470177292823792E-13, ... 4.5474735088646411896E-13, ... 2.2737367544323205948E-13, ... 1.1368683772161602974E-13, ... 5.6843418860808014870E-14, ... 2.8421709430404007435E-14, ... 1.4210854715202003717E-14, ... 7.1054273576010018587E-15, ... 3.5527136788005009294E-15, ... 1.7763568394002504647E-15, ... 8.8817841970012523234E-16, ... 4.4408920985006261617E-16, ... 2.2204460492503130808E-16, ... 1.1102230246251565404E-16, ... 5.5511151231257827021E-17, ... 2.7755575615628913511E-17, ... 1.3877787807814456755E-17, ... 6.9388939039072283776E-18, ... 3.4694469519536141888E-18, ... 1.7347234759768070944E-18 ]; if ( 1.0 < abs ( t ) ) fprintf ( 1, '\n' ); fprintf ( 1, 'ARCCOS_CORDIC - Fatal error!\n' ); fprintf ( 1, ' Input argument 1 < |T|.\n' ); error ( 'ARCCOS_CORDIC - Fatal error!' ); end theta = 0.0; z = [ 1.0, 0.0 ]'; poweroftwo = 1.0; for j = 1 : n if ( z(2) < 0.0 ) sign_z2 = -1.0; else sign_z2 = +1.0; end if ( t <= z(1) ) sigma = + sign_z2; else sigma = - sign_z2; end if ( j <= length ( angles ) ) angle = angles(j); else angle = angle / 2.0; end r = [ 1.0, -sigma * poweroftwo; ... sigma * poweroftwo, 1.0 ]; z = r * r * z; theta = theta + 2.0 * sigma * angle; t = t + t * poweroftwo * poweroftwo; poweroftwo = poweroftwo / 2.0; end return end

Практична робота Програмування

2015

Коментарі

Ви не можете залишити коментар. Для цього, будь ласка, або зареєструйтесь.

Ділись своїми роботами та отримуй миттєві бонуси!

поділитись