我只看了你的第二個(gè)
主要從事網(wǎng)頁設(shè)計(jì)、PC網(wǎng)站建設(shè)(電腦版網(wǎng)站建設(shè))、wap網(wǎng)站建設(shè)(手機(jī)版網(wǎng)站建設(shè))、自適應(yīng)網(wǎng)站建設(shè)、程序開發(fā)、微網(wǎng)站、小程序開發(fā)等,憑借多年來在互聯(lián)網(wǎng)的打拼,我們?cè)诨ヂ?lián)網(wǎng)網(wǎng)站建設(shè)行業(yè)積累了豐富的網(wǎng)站設(shè)計(jì)、成都做網(wǎng)站、網(wǎng)絡(luò)營(yíng)銷經(jīng)驗(yàn),集策劃、開發(fā)、設(shè)計(jì)、營(yíng)銷、管理等多方位專業(yè)化運(yùn)作于一體,具備承接不同規(guī)模與類型的建設(shè)項(xiàng)目的能力。
你定義的f函數(shù)只有一個(gè)形參x,所以在main函數(shù)中輸入的你所認(rèn)為的實(shí)參abcd 的值無法轉(zhuǎn)給f函數(shù) 所以運(yùn)行時(shí)會(huì)出錯(cuò) 只要賦給abcd確定的值或增加f函數(shù)的形參
scanf("%f,%f,%f,%f",d,d,d,d);這句也錯(cuò)了應(yīng)該是scanf("%f,%f,%f,%f",a,b,c,d);
S*cos(α+0.01410087457*S*cosα)=K
令b=0.01410087457,u=α+b*S*cosα,
則S=(u-α)/(bcosα),原方程變?yōu)?(u-α)cosu=bKcosα
令f(u)=(u-α)cosu,則f'(u)=cosu-(u-α)sinu
∵f(u)=(u-α)cosu=bKcosα=常數(shù),∴f'(u)=cosu-(u-α)sinu=0
聯(lián)立上述兩個(gè)方程,得 sin2u+bKcosα*sinu-1=0
∵△=(bKcosα)^2+4≥4,∴√△≥2
∴sinu=(-bKcosα+√△)/2 (另一根≤-1舍棄)
∴u=arcsin[(-bKcosα+√△)/2]
(若你的解不是唯一,則應(yīng)取u=nπ±arcsin[(-bKcosα+√△)/2],n∈Z)
∴S=(u-α)/(bcosα)={arcsin[(-bKcosα+√△)/2]-α}/(bcosα)
(其中△=(bKcosα)^2+4,b=0.01410087457;其余α,K均為已知)
用這個(gè)公式就可以直接算出你要求的S的值了
#includestdio.h
#includemath.h
void?main()
{
float?x0,x1,x2,fx0,fx1,fx2;
do
{
printf("enter?x1??x2:");
scanf("%f,%f",x1,x2); //輸入在哪個(gè)區(qū)間尋找方程的根
fx1=x1*((2*x1-4)*x1+3)-6; //計(jì)算在兩個(gè)端點(diǎn),函數(shù)的值
fx2=x2*((2*x2-4)*x2+3)-6;
}while(fx1*fx20); //必須要函數(shù)在兩個(gè)端點(diǎn),函數(shù)的值異號(hào),才能保證區(qū)間內(nèi)有根
do
{
x0=(x1+x2)/2; //x0為區(qū)間的中點(diǎn)
fx0=x0*((2*x0-4)*x0+3)-6; //中點(diǎn)的函數(shù)值
if((fx0*fx1)0) //將中點(diǎn)修改為函數(shù)值與中點(diǎn)函數(shù)值同號(hào)的端點(diǎn)
{
x2=x0;
fx2=fx0;
}
else
{
x1=x0;
fx2=fx0;
}
}while(fabs(fx0)=1e-5); //直到中點(diǎn)的函數(shù)值與0足夠接近
printf("x=%6.2f\n",x0);
}
可以使用二分法,弦截法,牛頓法,以前寫的,希望對(duì)LZ有幫助。
①弦截法:
:
用弦截法求方程F(x)=X^3-5X^2+16X-80=0的根。
方法如下:
(1)取兩個(gè)不同的點(diǎn)x1,x2,如果f(x1)和f(x2)符號(hào)相反,則(x1,x2)區(qū)間內(nèi)必有一個(gè)根。如果f(x1)和f(x2)同符號(hào),則應(yīng)改變x1,x2,直到f(x1)和f(x2)異號(hào)為止。注意x1和x2的值不應(yīng)差太大,以保證(x1,x2)區(qū)間內(nèi)只有一根。
(2)連接(x1,f(x1))和(x2,f(x2))兩點(diǎn),此線交x軸于x,x的坐標(biāo)可以用下式求出:x=(x1*f(x2)-x2*f(x1))/f(x2)-f(x1),
再從x求出f(x);
(3)若f(x)和f(x1)同符號(hào),則根必在(x,x2)內(nèi),此時(shí)將x作為新的x1;如果f(x)和f(x2)同符號(hào),則表示根在(x1,x),此時(shí)將x作為新的x2。
(4)重復(fù)步驟(2)和(3),知道|f(x)|e,e是一個(gè)很小的數(shù),例如10^-6;
程序如下:
#includestdio.h
#includemath.h
float f(float x)
{
float y;
y=((x-5.0)*x+16.0)*x-80.0;
return y;
}
float xpoint(float x1,float x2)
{//求(x1,f(x1)和(x2,f(x2))連線與x軸的交點(diǎn)
float y;
y=(x1*f(x2)-x2*f(x1))/(f(x2)-f(x1));
return y;
}//xpoint函數(shù)
float root(float x1,float x2)
{
float x,y,y1;
y1=f(x1);
do
{
x=xpoint(x1,x2);
y=f(x);
if(y*y10)
{
y1=y;
x1=x;
}
else x2=x;
}while (fabs(y)=0.00001);
return x;
}//root函數(shù)
void main()
{
float x1,x2,f1,f2,x;
do
{
printf("input x1,x2:\n");
scanf("%f%f",x1,x2);
f1=f(x1);
f2=f(x2);
} while(f1*f2=0);
x=root(x1,x2);
printf("a root of equation is %8.4f\n",x);
getch();
}
②:牛頓迭代法:
牛頓迭代法
牛頓迭代法又稱牛頓切線法,它采用以下方法求根:先任意設(shè)定一個(gè)與真實(shí)的根接近的值x0作為第一個(gè)近似根,由x0求出f(x0),過(x0,f(x0))點(diǎn)做f(x)的切線,交x軸于x1,把它作為第二次近似根,再由x1求出f(x1),再過(x1,f(x1))點(diǎn)做f(x)的切線,交x軸于x2,再求出f(x2),再作切線……如此繼續(xù)下去,直到足夠接近真正的x為止。
其中f'(X0)是函數(shù)在X0處的斜率,也就是在X0處的導(dǎo)數(shù)。
代碼如下:
#includestdio.h
#includemath.h
float f(float a,float b,float c,float d,float x)
{
float f;
f=((a*x+b)*x+c)*x+d;
return f;
}
float f1(float a,float b,float c,float x)
{
float f;
f=(x*3*a+2*b)*x+c;
return f;
}
float root(float a,float b,float c,float d)
{
float x0,x1=1;
do
{
x0=x1;
x1=x0-f(a,b,c,d,x0)/f1(a,b,c,x0);
}while(fabs(x1-x0)=1e-6);
return x0;
}
void main()
{
float a,b,c,d,x;
printf("input four float numbers:\n");
scanf("%f%f%f%f",a,b,c,d);
x=root(a,b,c,d);
printf("%.1fX^3+%.1fX^2+%.1fX+%.1f=0 its root near x=1.5 is :%.4f\n",a,b,c,d,x);
getch();
}
求sin的:參考下 #includestdio.h void main() { double x,a,b,sum=0; printf("請(qǐng)輸入x的弧度值:\n"); scanf("%lf",x); int i,j,count=0; for(i=1;;i+=2) { count++; a=b=1; for(j=1;j=i;j++) { a*=x; b*=(double)j; } if(a/b0.0000001) break; else { if(count%2==0) sum-=a/b; else sum+=a/b; } } printf("%lf\n",sum); }
要利用根與根之差的絕對(duì)值=1這個(gè)條件。
下面是剛學(xué)的時(shí)候?qū)懙模憧梢钥纯?/p>
#includestdio.h
#includemath.h
int a,b,c,d;
int count=0;
double x[3]={-200,-200,-200};
void fun(double x1,double x2)
{
double x3,y1,y2,y3;
int i=0;
while (x2-x110e-6)
{
x3=(x1+x2)/2;
y1=a*x1*x1*x1+b*x1*x1+c*x1+d;
y2=a*x2*x2*x2+b*x2*x2+c*x2+d;
y1=a*x1*x1*x1+b*x1*x1+c*x1+d;
y3=a*x3*x3*x3+b*x3*x3+c*x3+d;
if (y1==0)
{
for (i=0;i3;i++)
{
if (x[i]==x1)
{
return;
}
if (x[i]==-200)
{
count++;
x[i]=x1;
return;
}
}
}
else if (y2==0)
{
for (i=0;i3;i++)
{
if (x[i]==x2)
{
return;
}
if (x[i]==-200)
{
count++;
x[i]=x2;
return;
}
}
}
else if (y3==0)
{
for (i=0;i3;i++)
{
if (x[i]==x3)
{
return;
}
if (x[i]==-200)
{
count++;
x[i]=x3;
return;
}
}
}
else if (y1*y30)
{
x2=x3;
}
else
{
x1=x3;
}
}
for (i=0;i3;i++)
{
if (x[i]==x1)
{
return;
}
if (x[i]==-200)
{
count++;
x[i]=x1;
return;
}
}
}
int main()
{
double x1,x2,y1,y2,x3,x4;
scanf("%d,%d,%d,%d",a,b,c,d);
for (x1=-100,x2=x1+1;x1100;x1+=1,x2+=1)
{
y1=a*x1*x1*x1+b*x1*x1+c*x1+d;
y2=a*x2*x2*x2+b*x2*x2+c*x2+d;
x3=x1;
x4=x2;
if (y1*y2=0)
{
fun(x3,x4);
if (count==3)
{
break;
}
}
}
for(int i=0;i3;i++)
{
printf("%.2lf ",x[i]);
}
printf("\n");
return 0;
}
分享標(biāo)題:c語言解一元三角函數(shù)方程 c語言實(shí)現(xiàn)三角函數(shù)
鏈接URL:http://aaarwkj.com/article20/dooogco.html
成都網(wǎng)站建設(shè)公司_創(chuàng)新互聯(lián),為您提供微信公眾號(hào)、品牌網(wǎng)站建設(shè)、商城網(wǎng)站、關(guān)鍵詞優(yōu)化、網(wǎng)站設(shè)計(jì)公司、移動(dòng)網(wǎng)站建設(shè)
聲明:本網(wǎng)站發(fā)布的內(nèi)容(圖片、視頻和文字)以用戶投稿、用戶轉(zhuǎn)載內(nèi)容為主,如果涉及侵權(quán)請(qǐng)盡快告知,我們將會(huì)在第一時(shí)間刪除。文章觀點(diǎn)不代表本網(wǎng)站立場(chǎng),如需處理請(qǐng)聯(lián)系客服。電話:028-86922220;郵箱:631063699@qq.com。內(nèi)容未經(jīng)允許不得轉(zhuǎn)載,或轉(zhuǎn)載時(shí)需注明來源: 創(chuàng)新互聯(lián)