PDA

View Full Version : Calcolo delle radici in C++


Rhapsodyan [ITA]
29-01-2005, 20:18
Ciao a tutti :)
Ho un problema: non riesco a calcolare la radice "ennesima" di un numero...Ad esempio, che funzione posso utilizzare per calcolare la radice 350sima di un numero? (è solo un esempio, mi serve una regola generica da utilizzare con la media geometrica)
Ho pensato a questa formula, che in excel funziona:
numero^(1/350)...eppure in c++ non funziona :muro:
Come posso fare?

sirus
29-01-2005, 20:52
mai pensato alla funzione sqrt(...) nel'header MATH.H ;) ?!
PS esiste una variante che calcola qualsiasi radice nel campo dei double :D

Rhapsodyan [ITA]
29-01-2005, 21:00
l'sqrt fa la radice quadrata..io ho bisogno di radici molto più alte;)

Ubi-Wan
29-01-2005, 21:01
sei sicuro che pow(x,1/n) non funzioni??
io ho provato la radice 100-esima e va benissimo...

cmq puoi anche usare exp(log(x)*1/n)

Rhapsodyan [ITA]
29-01-2005, 21:18
Ora sembra funzionare..ma solo coi numeri bassi..eppure non capisco..ho dato double a tutti :muro:
E' la mia prima applicazione MCF,fino ad oggi lavoravo in console..
Allego qui il sorgente,scusate la lunghezza..

// Prova 3Dlg.cpp : implementation file
//

#include "stdafx.h"
#include "Prova 3.h"
#include "Prova 3Dlg.h"
#include "iostream.h"
#include "math.h"

#ifdef _DEBUG
#define new DEBUG_NEW
#undef THIS_FILE
static char THIS_FILE[] = __FILE__;
#endif

long int j=0,sommayarmonica=0,ctr_esponente=0;
double media,somma=0,xy=0,sommay=0,rapporto,frequenze[20000],ctr_esponente2=0,valori[20000],rapporti=0,moltiplica,xelevato;
/////////////////////////////////////////////////////////////////////////////
// CAboutDlg dialog used for App About

class CAboutDlg : public CDialog
{
public:
CAboutDlg();

// Dialog Data
//{{AFX_DATA(CAboutDlg)
enum { IDD = IDD_ABOUTBOX };
//}}AFX_DATA

// ClassWizard generated virtual function overrides
//{{AFX_VIRTUAL(CAboutDlg)
protected:
virtual void DoDataExchange(CDataExchange* pDX); // DDX/DDV support
//}}AFX_VIRTUAL

// Implementation
protected:
//{{AFX_MSG(CAboutDlg)
//}}AFX_MSG
DECLARE_MESSAGE_MAP()
};

CAboutDlg::CAboutDlg() : CDialog(CAboutDlg::IDD)
{
//{{AFX_DATA_INIT(CAboutDlg)
//}}AFX_DATA_INIT
}

void CAboutDlg::DoDataExchange(CDataExchange* pDX)
{
CDialog::DoDataExchange(pDX);
//{{AFX_DATA_MAP(CAboutDlg)
//}}AFX_DATA_MAP
}

BEGIN_MESSAGE_MAP(CAboutDlg, CDialog)
//{{AFX_MSG_MAP(CAboutDlg)
// No message handlers
//}}AFX_MSG_MAP
END_MESSAGE_MAP()

/////////////////////////////////////////////////////////////////////////////
// CProva3Dlg dialog

CProva3Dlg::CProva3Dlg(CWnd* pParent /*=NULL*/)
: CDialog(CProva3Dlg::IDD, pParent)
{
//{{AFX_DATA_INIT(CProva3Dlg)
m_sceltamedia = _T("Aritmetica");
m_frequenza = 0.0;
m_valore = 0.0;
m_media = 0.0;
//}}AFX_DATA_INIT
// Note that LoadIcon does not require a subsequent DestroyIcon in Win32
m_hIcon = AfxGetApp()->LoadIcon(IDR_MAINFRAME);
}

void CProva3Dlg::DoDataExchange(CDataExchange* pDX)
{
CDialog::DoDataExchange(pDX);
//{{AFX_DATA_MAP(CProva3Dlg)
DDX_CBString(pDX, IDC_tendina, m_sceltamedia);
DDX_Text(pDX, IDC_EDIT1, m_frequenza);
DDX_Text(pDX, IDC_EDIT2, m_valore);
DDX_Text(pDX, IDC_EDIT3, m_media);
//}}AFX_DATA_MAP
}

BEGIN_MESSAGE_MAP(CProva3Dlg, CDialog)
//{{AFX_MSG_MAP(CProva3Dlg)
ON_WM_SYSCOMMAND()
ON_WM_PAINT()
ON_WM_QUERYDRAGICON()
ON_BN_CLICKED(IDOK, OnMedia)
ON_BN_CLICKED(IDC_inserisci, OnInserisci)
ON_BN_CLICKED(IDazzera, Onazzera)
//}}AFX_MSG_MAP
END_MESSAGE_MAP()

/////////////////////////////////////////////////////////////////////////////
// CProva3Dlg message handlers

BOOL CProva3Dlg::OnInitDialog()
{
CDialog::OnInitDialog();

// Add "About..." menu item to system menu.

// IDM_ABOUTBOX must be in the system command range.
ASSERT((IDM_ABOUTBOX & 0xFFF0) == IDM_ABOUTBOX);
ASSERT(IDM_ABOUTBOX < 0xF000);

CMenu* pSysMenu = GetSystemMenu(FALSE);
if (pSysMenu != NULL)
{
CString strAboutMenu;
strAboutMenu.LoadString(IDS_ABOUTBOX);
if (!strAboutMenu.IsEmpty())
{
pSysMenu->AppendMenu(MF_SEPARATOR);
pSysMenu->AppendMenu(MF_STRING, IDM_ABOUTBOX, strAboutMenu);
}
}

// Set the icon for this dialog. The framework does this automatically
// when the application's main window is not a dialog
SetIcon(m_hIcon, TRUE); // Set big icon
SetIcon(m_hIcon, FALSE); // Set small icon

// TODO: Add extra initialization here

return TRUE; // return TRUE unless you set the focus to a control
}

void CProva3Dlg::OnSysCommand(UINT nID, LPARAM lParam)
{
if ((nID & 0xFFF0) == IDM_ABOUTBOX)
{
CAboutDlg dlgAbout;
dlgAbout.DoModal();
}
else
{
CDialog::OnSysCommand(nID, lParam);
}
}

// If you add a minimize button to your dialog, you will need the code below
// to draw the icon. For MFC applications using the document/view model,
// this is automatically done for you by the framework.

void CProva3Dlg::OnPaint()
{
if (IsIconic())
{
CPaintDC dc(this); // device context for painting

SendMessage(WM_ICONERASEBKGND, (WPARAM) dc.GetSafeHdc(), 0);

// Center icon in client rectangle
int cxIcon = GetSystemMetrics(SM_CXICON);
int cyIcon = GetSystemMetrics(SM_CYICON);
CRect rect;
GetClientRect(&rect);
int x = (rect.Width() - cxIcon + 1) / 2;
int y = (rect.Height() - cyIcon + 1) / 2;

// Draw the icon
dc.DrawIcon(x, y, m_hIcon);
}
else
{
CDialog::OnPaint();
}
}

// The system calls this to obtain the cursor to display while the user drags
// the minimized window.
HCURSOR CProva3Dlg::OnQueryDragIcon()
{
return (HCURSOR) m_hIcon;
}

void CProva3Dlg::OnMedia()
{
UpdateData(true);
if (m_sceltamedia=="Aritmetica")
{
somma=0;
sommay=0;

for (int i=0;i<j;i++)
{
xy=valori[i]*frequenze[i];
somma+=xy;
sommay+=frequenze[i];
}
m_media=somma/sommay;
if (sommay!=0)
UpdateData(false);
else
MessageBox("Impossibile calcolare la media: hai inserito i valori corretti?","Errore");
}
if (m_sceltamedia=="Armonica")
{
sommayarmonica=0;
rapporti=0;

for (int i=0;i<j;i++)
{
sommayarmonica+=frequenze[i];
rapporti+=frequenze[i]/valori[i];
}
m_media=sommayarmonica/rapporti;
if (rapporti!=0)
UpdateData(false);
else
MessageBox("Impossibile calcolare la media: hai inserito i valori corretti?","Errore");
}
if (m_sceltamedia=="Geometrica")
{
moltiplica=1;
sommay=0;

for (int i=0;i<j;i++)
{
xelevato=pow(valori[i],frequenze[i]);
moltiplica*=xelevato;
sommay+=frequenze[i];
}
m_media=pow(moltiplica,1/sommay);
if (sommay!=0)
UpdateData(false);
else
MessageBox("Impossibile calcolare la media: hai inserito i valori corretti?","Errore");

}
if (m_sceltamedia=="Quadratica")
{
somma=0;
sommay=0;

for (int i=0;i<j;i++)
{
xy=valori[i]*valori[i]*frequenze[i];
somma+=xy;
sommay+=frequenze[i];
}
m_media=sqrt(somma/sommay);
if (sommay!=0)
UpdateData(false);
else
MessageBox("Impossibile calcolare la media: hai inserito i valori corretti?","Errore");
}

}

void OnUpdateEdit2()
{
}
void CProva3Dlg::OnInserisci()
{
UpdateData(true);

if (m_frequenza!=0 && m_valore!=0)
{
frequenze[j]=m_frequenza;
valori[j]=m_valore;
m_frequenza=0; //azzeramento per visualizzazione
m_valore=0;
UpdateData(false);
j++;
}
else
MessageBox("Errore: reinserisci i dati.","Errore");
}

void CProva3Dlg::OnCancel()
{

CDialog::OnCancel();
}

void CProva3Dlg::Onazzera()
{
j=0;
media=0;
somma=0;
xy=0;
sommay=0;
m_valore=0;
m_frequenza=0;
m_media=0;
UpdateData(false);

}

Rhapsodyan [ITA]
29-01-2005, 22:48
Aggiungo anche il programma per chi volesse provarlo:O

supermario
30-01-2005, 10:49
mancano i tuoi header

Rhapsodyan [ITA]
30-01-2005, 11:03
mancano i tuoi header
Non serve, ora funziona,grazie :)
Però funziona solo con i numeri non troppo alti, dopodichè sballa i risultati...ho dato alle variabili ecc tipo double, forse dovrei dargli long?

supermario
30-01-2005, 11:13
boh ho scaricato il programma ma all'inizio mi chiedeva 2 dll........le ho scaricate e installate ma ora mi da errore nella mfc42d

Rhapsodyan [ITA]
30-01-2005, 11:21
Originariamente inviato da supermario
boh ho scaricato il programma ma all'inizio mi chiedeva 2 dll........le ho scaricate e installate ma ora mi da errore nella mfc42d
Cazz, ancora questo problema...l'ho passato ad un amico che ha c++6 (come il mio) e gli funziona senza problemi...l'ho passato ad un altro che ha il .NET e come te gli mancano delle dll...Lui dice che è perchè l'ho compilato in modalità debug, ma io non sò manco come si fa...:muro:

supermario
30-01-2005, 11:24
Originariamente inviato da Rhapsodyan [ITA]
Cazz, ancora questo problema...l'ho passato ad un amico che ha c++6 (come il mio) e gli funziona senza problemi...l'ho passato ad un altro che ha il .NET e come te gli mancano delle dll...Lui dice che è perchè l'ho compilato in modalità debug, ma io non sò manco come si fa...:muro:


appunto io ho il .net

magati se mi passi tutti gli header e il sorgente lo posso compilare su .net e vediamo se va:)

Rhapsodyan [ITA]
30-01-2005, 14:13
Certo,se mi dai l'email te lo mando, sono quasi 2MB.
Grazie per la disponibilità :)

supermario
30-01-2005, 16:07
hai PVT;)

Rhapsodyan [ITA]
30-01-2005, 16:15
Inviato:D

supermario
30-01-2005, 16:18
arrivato :D


creo subito;)

supermario
30-01-2005, 16:29
umh mi dava 3 errori e li ho risolti...ma ne sn usciti altri 3 interni al programma e nn so come risolverli

ho messo nei source

Prova 3.cpp
Prova 3dlg.cpp
StdAfx.cpp

e negli header i rispettivi .h

nel resource

Prova 3.rc

ma nn va:cry:

Rhapsodyan [ITA]
30-01-2005, 16:34
Errori? Strano, a me dà solo un warning ma funziona...incompatibilità?

supermario
30-01-2005, 16:45
nn saprei.........ecco gli errori


http://img111.exs.cx/img111/7338/err0mb.jpg

Rhapsodyan [ITA]
30-01-2005, 17:40
Non ne ho idea:( :cry:

Rhapsodyan [ITA]
31-01-2005, 20:27
Nessuno sà aiutarmi?