c语言求直线方程
❶ c语言求直线方程最简式
intgcd(inta,intb){
intk;
if(a>b){
k=a;a=b;b=k;
}
while(1){
intk=a%b;
if(k==0)returnb;
a=b;
b=k;
}
}
voidsimplify(int*coef,int*out){
intg=gcd(coef[0],coef[1]);
out[0]=coef[0]/g;
g=coef[1]/g;
out[1]=-g;
out[2]=coef[2]*g;
if(out[0]==1)printf("x");
elseprintf("%dx",out[0]);
if(out[1]==-1)printf("-y");
elseprintf("%dy",out[1]);
printf("+%d=0
",out[2]);
}
intmain()
{
intn=0,i;
intcoef[3],out[3];
charstr[256];
scanf("%d",&n);
for(i=0;i<n;i++){
scanf("y=%d/%dx+%d",&coef[0],&coef[1],&coef[2]);
simplify(coef,out);
}
return0;
}
❷ C语言 解线性方程组
给，下面的代码已经编译运行确认，肯定好用了，试试吧:)
#include<conio.h>
#include<stdio.h>
#include<math.h>
#include<stdlib.h>
#define NUMBER 20
#define Esc 0x1b
#define Enter 0x0d
float A[NUMBER][NUMBER+1] ,ark;
int flag,n;
void exchange(int r,int k);
float max(int k);
void message();
int main()
{
float x[NUMBER]; /*此数组用于存放方程解*/
int k,i,j;
char celect;
system("cls");
printf("\n用Gauss列主元消元法解线性方程组");
printf("\n1.解方程组请按Enter.");
printf("\n2.退出程式请按Esc.");
celect=getch();
if(celect==Esc)
exit(0);
printf("\n 输入方程组的维数：n=");
scanf("%d",&n);
printf("\n现在输入系数矩阵A和向量b:");
for(i=1;i<=n;i++)
{
printf("\n请输入a%d1--a%d%d系数和向量b%d: \n",i,i,n,i);
/*实现将每一行中的系数和向量一次性输入，数之间用空格格开，输完后回车确定*/
for(j=1;j<=n+1;j++) /*将刚才输入的数存入数组*/
scanf("%f",&A[i][j]);
}
for(k=1;k<=n-1;k++)
{
ark=max(k);
if(ark==0) /*判断方程是否为线性方程，即是否合法*/
{
printf("\n此方程组不合法!");message();
}
else if(flag!=k)
exchange(flag,k);
for(i=k+1;i<=n;i++)
for(j=k+1;j<=n+1;j++)
A[i][j]=A[i][j]-A[k][j]*A[i][k]/A[k][k];
}
x[n]=A[n][n+1]/A[n][n];
for( k=n-1;k>=1;k--)
{
float me=0;
for(j=k+1;j<=n;j++)
{
me=me+A[k][j]*x[j];
}
x[k]=(A[k][n+1]-me)/A[k][k];
}
for(i=1;i<=n;i++)
{
printf("\nx%d=%f",i,x[i]);
}
message();
getch();
return 1;
}
void exchange(int r,int k) /*交换行的矩函数*/
{
int i;
for(i=1;i<=n+1;i++)
A[0][i]=A[r][i];
for(i=1;i<=n+1;i++)
A[r][i]=A[k][i];
for(i=1;i<=n+1;i++)
A[k][i]=A[0][i];
}
float max(int k) /*比校系数大小的函数*/
{
int i;
float temp=0;
for(i=k;i<=n;i++)
if(fabs(A[i][k])>temp)
{
temp=fabs(A[i][k]);
flag=i;
}
return temp;
}
void message() /*实现菜单选择的函数*/
{
printf("\n 继续运算按 Enter ,退出程式按 Esc!");
switch(getch())
{
case Enter: main();
case Esc: exit(0);
default:{printf("\n不合法的输入!");message();}
}
}
❸ C语言中用什么函数求方程参量
这里需要知道一个隐含的条件：ax+by+c=0是一条直线。两点确定一条直线，所以这里需要充分利用这两条来求直线。
表示直线的方式有好多种，ax+by+c=0只是其中之一。因为已经知道两点的坐标，我们采用两点式来表示直线。
令：(x0,y0)=(0, 0), (x1,y1)=(1,2)
则直线方程为：(y-y0)/(y0-y1)=(x-x0)/(x0-x1)
代入(y-0)/(0-1) = (x-0)/(0-2)
整理之后，得到直线方程：x - 2y = 0;
所以a=1, b=-2, c=0;方程直线为x - 2y.
❹ C语言编程，直线方程y=kx+3，要求输入y坐标以及k，求x坐标。
main()
{
基本语言我忘记了，用中文描述好了
定义Y，K，x的值为float；
scanf函数给Y,K变量复制；
X=(Y-3)/K；
printf函数输出X变量；
}
❺ c语言求直线方程
应该是
scanf("%lf %lf %lf %lf",&x1,&y1,&x2,&y2);
double 类型对应的输出格式为%lf，格式不匹配可能会出错。
不是"/"的问题，因为y2-y1和x2-x1均为double类型，因此这里不是整除。
❻ 用C语言编写 在平面上,穿过两个点(x1,y1)和(x2,y2)的直线是由下面的方程定义的:
p1和p2的坐标代入直线方程中，得到：
c1=a*x1+b*y1;
c2=a*xx2+by2;
如果c1和c2都小于c 或者c1和c2都大于c说明两个点位于直线的同一侧。
主要是数学问题，高中数学。
❼ 用c语言程序设计，求代码
提供两种思路，这里只给思路不给代码，因为网络知道现在似乎贴代码格式会乱，而且自己写代码更容易进步。
第一种暴力的方法：遍历所有格点，用勾股定理算出三条边长度，用海伦公式算出面积。
勾股定理这个大家都会。海伦公式可以通过三边长度算出三角形面积。
海伦公式大致内容：令 p = ( a + b + c ) / 2 ，则 S = sqrt [ p ( p - a ) ( p -b ) ( p - c ) ] 。
（之前题目看错了。。重新写一遍这个方法
暴力的方法在 9×9 的范围内当然跑得快，但是如果题目加强，让你在 10^7×10^7 的方格内找好点，你怎么找？
第二种方法：根据面积相等和底边长度关系，可以求出高的比值，然后设出其中一条高的长度，就能得到另一条高的长度。然后写出两个直线方程并联立，可以得到一个新的方程。这个方程的图象上的整点就是好点。
听起来有点绕，拿这道题做例子就好理解了。
设 △PAB 的高为 a （为了方便记作条件①），则 △PAC 的高为 2a （条件②）。
以方格左下角为原点建系，可以写出满足条件①的直线解析式 y = x + √2 a 。
同理满足条件②的直线解析式为 y = 12 - x - 2√2 a 。（这里 a 的系数可正可负，篇幅原因只算一种，另一种后面代入一遍就好了）
联立可得 P 坐标为 ( 6 - 3√2/2 a, 6 - √2/2 a ）。
根据坐标式子和题目性质可得 P 图象必定过 A（ 6, 6 ），并且图象是一条直线。
设 P : y = k ( x - 6 ) + 6 ，代入坐标解得 k = 1/3 。
所以 P 在直线 y = 1/3 x + 4 上。
同理，P 也在直线 y = 3 x - 12 上。
然后把所有横坐标代入就好啦，复杂度直接少了一阶。（事实上还可以继续优化，不过这样已经很优了）
❽ 雅克比迭代法求解线性方程组的C语言程序
void Solve ( double dCoef [] , double dY [] , unsigned int iOrder , double dErr)
{//用Jacobi迭代法解方程组, dCoef[]系数阵, Y[]向量, iOrder给出方程阶数, dErr给出精度
double res [Max]; //方程解
double res2[Max]; //保存上一阶方程解
if ( Max < iOrder )
{
printf ("最多支持%d阶方程组.", Max);
return;
}
for ( unsigned int i = 0 ; i < iOrder ; res2 [i++] = 0.0 ); //初始解向量 (0,0...)
while ( true )
{
bool bStopIterative = true;
for (unsigned int i = 0 ; i < iOrder ; ++i)
{
double dSum2 = 0;
for (unsigned int j = 0 ; j < iOrder ; j++)
{//求第二项
if ( j == i ) continue;
dSum2 += dCoef [i * iOrder + j] * res2 [j];
}
res[i] = 1/dCoef[i * iOrder + i] * ( dY[i] - dSum2 );
if ( abs ( res2[i] - res [i] ) > dErr )
bStopIterative = false;
}
if ( bStopIterative )
break;
for (unsigned int i = 0 ; i < iOrder ; i++ )
res2[ i ] = res[ i ];
}
//输出结果
for (unsigned int i = 0 ; i < iOrder ; i++)
printf ("x%d = %lf\n", i+1 , res[i]);
}
int main(int argc, char* argv[])
{
double a[] =
{
8,-3,2,
4,11,-1,
2,1,4
};
double b[3] =
{
20,33,12
};
Solve ( a , b , 15 , 1e-10);
getchar();
return 0;
}
❾ c语言解线性方程组的编程题 【做的好会多给分】
以下算法的适用条件：A的各阶主子式不为零
另外还可以采用
直接法：
消元法：Gauss-Jordan消元法，
分解法：Dolittle分解 （我用的是Courant分解法），追赶法，对称正定矩阵的LDL‘分解
----------
迭代法：
Jacobi迭代
Gauss-Seidel迭代
松弛迭代
-----------------
你上网可以搜索一下，或者看看数值计算方面的书
OK, 你看看这个, 另外还加了注释 ：
Courant分解算法：
aX = b, 作 A=LU, L是下三角矩阵， U是上三角矩阵
即L =
| L11
| L21 L22
| L31 L32 L33
| ..............
| Ln1 Ln2 ........Lnn
U =
| 1 U12 ..... U1n
| 空格 1 ..... U2n
| 空格 空格 ........
| 空格 空格 空格 空格 空格1
---------------------------------------------------
aX = b -----> LUX = b
记 UX = y,
由Ly = b得到
因为无法输出数学符号，以下采用[i, j]Ai 表示对Ai从i到j求和
yi = (bi - [j=1, i-1]Lij yj) / Lii i = 1, 2, ..., n
由UX = y得到
xi = yi - [j=i+1, n]uij xj j = n, n-1, ..., 2, 1
你在纸上验证一下就明白了
--------------------------------------------------------------
以下采用Courant分解 解 aX = b, 经检查，程序运行正确
这是运行结果：
--------------------------------------------------------------
Input n value(dim of Ax=b): 3
Now input the matrix a(i, j), i, j = 0, ..., 2:
1 2 1 -2 -1 -5 0 -1 6
Now input the matrix b(i), i = 0, ..., 2:
24 -63 50
Solve...x_i =
7.000000
4.000000
9.000000
--------------------------------------------------------------
#include "stdafx.h"
#include <stdio.h>
#include <math.h>
#define MAX_N 20
int main(int argc, char* argv[])
{
int n; // 未知数个数
int i, j, k;
static double a[MAX_N][MAX_N], b[MAX_N], x[MAX_N], y[MAX_N];
static double l[MAX_N][MAX_N], u[MAX_N][MAX_N];
printf("\nInput n value(dim of Ax=b): ");
scanf("%d", &n);
if(n >MAX_N)
{
printf("The input n is larger than MAX_N, please redefine the MAX_N.\n");
return 1;
}
if(n <= 0)
{
printf("Please input a number between 1 and %d.\n", MAX_N);
return 1;
}
// {{ 程序输入
printf("Now input the matrix a(i, j), i, j = 0, ..., %d:\n", n-1);
for (i=0; i<n; i++)
for (j=0; j<n; j++)
scanf("%lf", &a[i][j]);
printf("Now input the matrix b(i), i = 0, ..., %d:\n", n-1);
for(i=0; i<n; i++)
scanf("%lf", &b[i]);
// }} 程序输入
for(i=0; i<n; i++)
u[i][i] = 1; //
for(k=0; k<n; k++)
{
for(i=k; i<n; i++) // 计算L的第k列元素
{
l[i][k] = a[i][k];
for(j=0; j<=k-1; j++)
l[i][k] -= (l[i][j]*u[j][k]);
}
for(j=k+1; j<n; j++) //计算U的第k行元素
{
u[k][j] = a[k][j];
for(i=0; i<=k-1; i++)
u[k][j] -= (l[k][i]*u[i][j]);
u[k][j] /= l[k][k];
}
}
for(i=0; i<n; i++) // 解Ly = b
{
y[i] = b[i];
for(j=0; j<=i-1; j++)
y[i] -= (l[i][j]*y[j]);
y[i] /= l[i][i];
}
for(i=n-1; i>=0; i--) // 解UX = Y
{
x[i]=y[i];
for(j=i+1; j<n; j++)
x[i] -= (u[i][j]*x[j]);
}
printf("Solve...x_i = \n"); // 输出结果
for(i=0; i<n; i++)
printf("%f\n", x[i]);
return 0;
}
❿ c语言解y=a/x+b
用x'代替x，就转化为直线方程了，最小二乘法有公式，自己求偏导推一下也很方便，10组数据量很小，直接计算就OK了。。