im trying to solve this system of 10 equations

2 views (last 30 days)
clc
clear all
M=10
M = 10
deltax=2*pi/M
deltax = 0.6283
syms f1 f2 f3 f4 f5 f6 f7 f8 f9 f10
U1=f1*exp(1i*deltax)+f2*exp(1i*2*deltax)+f3*exp(1i*3*deltax)+f4*exp(1i*4*deltax)+f5*exp(1i*5*deltax)+f6*exp(1i*6*deltax)+f7*exp(1i*7*deltax)+f8*exp(1i*8*deltax)+f9*exp(1i*9*deltax)+f10*exp(1i*10*deltax)==cos(deltax)
U1 = 
U2=f1*exp(2*1i*deltax)+f2*exp(2*1i*2*deltax)+f3*exp(2*1i*3*deltax)+f4*exp(2*1i*4*deltax)+f5*exp(2*1i*5*deltax)+f6*exp(2*1i*6*deltax)+f7*exp(2*1i*7*deltax)+f8*exp(2*1i*8*deltax)+f9*exp(2*1i*9*deltax)+f10*exp(2*1i*10*deltax)==cos(2*deltax)
U2 = 
U3=f1*exp(3*1i*deltax)+f2*exp(3*1i*2*deltax)+f3*exp(3*1i*3*deltax)+f4*exp(3*1i*4*deltax)+f5*exp(3*1i*5*deltax)+f6*exp(3*1i*6*deltax)+f7*exp(3*1i*7*deltax)+f8*exp(3*1i*8*deltax)+f9*exp(3*1i*9*deltax)+f10*exp(3*1i*10*deltax)==cos(3*deltax)
U3 = 
U4=f1*exp(4*1i*deltax)+f2*exp(4*1i*2*deltax)+f3*exp(4*1i*3*deltax)+f4*exp(4*1i*4*deltax)+f5*exp(4*1i*5*deltax)+f6*exp(4*1i*6*deltax)+f7*exp(4*1i*7*deltax)+f8*exp(4*1i*8*deltax)+f9*exp(4*1i*9*deltax)+f10*exp(4*1i*10*deltax)==cos(4*deltax)
U4 = 
U5=f1*exp(5*1i*deltax)+f2*exp(5*1i*2*deltax)+f3*exp(5*1i*3*deltax)+f4*exp(5*1i*4*deltax)+f5*exp(5*1i*5*deltax)+f6*exp(5*1i*6*deltax)+f7*exp(5*1i*7*deltax)+f8*exp(5*1i*8*deltax)+f9*exp(5*1i*9*deltax)+f10*exp(5*1i*10*deltax)==cos(5*deltax)
U5 = 
U6=f1*exp(6*1i*deltax)+f2*exp(6*1i*2*deltax)+f3*exp(6*1i*3*deltax)+f4*exp(6*1i*4*deltax)+f5*exp(6*1i*5*deltax)+f6*exp(6*1i*6*deltax)+f7*exp(6*1i*7*deltax)+f8*exp(6*1i*8*deltax)+f9*exp(6*1i*9*deltax)+f10*exp(6*1i*10*deltax)==cos(6*deltax)
U6 = 
U7=f1*exp(7*1i*deltax)+f2*exp(7*1i*2*deltax)+f3*exp(7*1i*3*deltax)+f4*exp(7*1i*4*deltax)+f5*exp(7*1i*5*deltax)+f6*exp(7*1i*6*deltax)+f7*exp(7*1i*7*deltax)+f8*exp(7*1i*8*deltax)+f9*exp(7*1i*9*deltax)+f10*exp(7*1i*10*deltax)==cos(7*deltax)
U7 = 
U8=f1*exp(8*1i*deltax)+f2*exp(8*1i*2*deltax)+f3*exp(8*1i*3*deltax)+f4*exp(8*1i*4*deltax)+f5*exp(8*1i*5*deltax)+f6*exp(8*1i*6*deltax)+f7*exp(8*1i*7*deltax)+f8*exp(8*1i*8*deltax)+f9*exp(8*1i*9*deltax)+f10*exp(8*1i*10*deltax)==cos(8*deltax)
U8 = 
U9=f1*exp(9*1i*deltax)+f2*exp(9*1i*2*deltax)+f3*exp(9*1i*3*deltax)+f4*exp(9*1i*4*deltax)+f5*exp(9*1i*5*deltax)+f6*exp(9*1i*6*deltax)+f7*exp(9*1i*7*deltax)+f8*exp(9*1i*8*deltax)+f9*exp(9*1i*9*deltax)+f10*exp(9*1i*10*deltax)==cos(9*deltax)
U9 = 
U10=f1*exp(10*1i*deltax)+f2*exp(10*1i*2*deltax)+f3*exp(10*1i*3*deltax)+f4*exp(10*1i*4*deltax)+f5*exp(10*1i*5*deltax)+f6*exp(10*1i*6*deltax)+f7*exp(10*1i*7*deltax)+f8*exp(10*1i*8*deltax)+f9*exp(10*1i*9*deltax)+f10*exp(10*1i*10*deltax)==cos(10*deltax)
U10 = 
linsolve(U1,U2,U3,U4,U5,U6,U7,U8,U9,U10,f1,f2,f3,f4,f5,f6,f7,f8,f9,f10)
Error using sym/linsolve
Too many input arguments.
im trying to solve this system of equations
any help would be great

Answers (1)

Walter Roberson
Walter Roberson on 6 Feb 2023
M=10
M = 10
deltax=2*sym(pi)/M
deltax = 
syms f1 f2 f3 f4 f5 f6 f7 f8 f9 f10
U1=f1*exp(1i*deltax)+f2*exp(1i*2*deltax)+f3*exp(1i*3*deltax)+f4*exp(1i*4*deltax)+f5*exp(1i*5*deltax)+f6*exp(1i*6*deltax)+f7*exp(1i*7*deltax)+f8*exp(1i*8*deltax)+f9*exp(1i*9*deltax)+f10*exp(1i*10*deltax)==cos(deltax);
U2=f1*exp(2*1i*deltax)+f2*exp(2*1i*2*deltax)+f3*exp(2*1i*3*deltax)+f4*exp(2*1i*4*deltax)+f5*exp(2*1i*5*deltax)+f6*exp(2*1i*6*deltax)+f7*exp(2*1i*7*deltax)+f8*exp(2*1i*8*deltax)+f9*exp(2*1i*9*deltax)+f10*exp(2*1i*10*deltax)==cos(2*deltax);
U3=f1*exp(3*1i*deltax)+f2*exp(3*1i*2*deltax)+f3*exp(3*1i*3*deltax)+f4*exp(3*1i*4*deltax)+f5*exp(3*1i*5*deltax)+f6*exp(3*1i*6*deltax)+f7*exp(3*1i*7*deltax)+f8*exp(3*1i*8*deltax)+f9*exp(3*1i*9*deltax)+f10*exp(3*1i*10*deltax)==cos(3*deltax);
U4=f1*exp(4*1i*deltax)+f2*exp(4*1i*2*deltax)+f3*exp(4*1i*3*deltax)+f4*exp(4*1i*4*deltax)+f5*exp(4*1i*5*deltax)+f6*exp(4*1i*6*deltax)+f7*exp(4*1i*7*deltax)+f8*exp(4*1i*8*deltax)+f9*exp(4*1i*9*deltax)+f10*exp(4*1i*10*deltax)==cos(4*deltax);
U5=f1*exp(5*1i*deltax)+f2*exp(5*1i*2*deltax)+f3*exp(5*1i*3*deltax)+f4*exp(5*1i*4*deltax)+f5*exp(5*1i*5*deltax)+f6*exp(5*1i*6*deltax)+f7*exp(5*1i*7*deltax)+f8*exp(5*1i*8*deltax)+f9*exp(5*1i*9*deltax)+f10*exp(5*1i*10*deltax)==cos(5*deltax);
U6=f1*exp(6*1i*deltax)+f2*exp(6*1i*2*deltax)+f3*exp(6*1i*3*deltax)+f4*exp(6*1i*4*deltax)+f5*exp(6*1i*5*deltax)+f6*exp(6*1i*6*deltax)+f7*exp(6*1i*7*deltax)+f8*exp(6*1i*8*deltax)+f9*exp(6*1i*9*deltax)+f10*exp(6*1i*10*deltax)==cos(6*deltax);
U7=f1*exp(7*1i*deltax)+f2*exp(7*1i*2*deltax)+f3*exp(7*1i*3*deltax)+f4*exp(7*1i*4*deltax)+f5*exp(7*1i*5*deltax)+f6*exp(7*1i*6*deltax)+f7*exp(7*1i*7*deltax)+f8*exp(7*1i*8*deltax)+f9*exp(7*1i*9*deltax)+f10*exp(7*1i*10*deltax)==cos(7*deltax);
U8=f1*exp(8*1i*deltax)+f2*exp(8*1i*2*deltax)+f3*exp(8*1i*3*deltax)+f4*exp(8*1i*4*deltax)+f5*exp(8*1i*5*deltax)+f6*exp(8*1i*6*deltax)+f7*exp(8*1i*7*deltax)+f8*exp(8*1i*8*deltax)+f9*exp(8*1i*9*deltax)+f10*exp(8*1i*10*deltax)==cos(8*deltax);
U9=f1*exp(9*1i*deltax)+f2*exp(9*1i*2*deltax)+f3*exp(9*1i*3*deltax)+f4*exp(9*1i*4*deltax)+f5*exp(9*1i*5*deltax)+f6*exp(9*1i*6*deltax)+f7*exp(9*1i*7*deltax)+f8*exp(9*1i*8*deltax)+f9*exp(9*1i*9*deltax)+f10*exp(9*1i*10*deltax)==cos(9*deltax);
U10=f1*exp(10*1i*deltax)+f2*exp(10*1i*2*deltax)+f3*exp(10*1i*3*deltax)+f4*exp(10*1i*4*deltax)+f5*exp(10*1i*5*deltax)+f6*exp(10*1i*6*deltax)+f7*exp(10*1i*7*deltax)+f8*exp(10*1i*8*deltax)+f9*exp(10*1i*9*deltax)+f10*exp(10*1i*10*deltax)==cos(10*deltax);
eqns = [U1;U2;U3;U4;U5;U6;U7;U8;U9;U10];
vars = [f1;f2;f3;f4;f5;f6;f7;f8;f9;f10];
[A,b] = equationsToMatrix(eqns, vars);
sol = A\b;
vars(:) == sol
ans = 

Products


Release

R2022b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!