matlab程式**
繪製[}e^sin(\\sqrtt)', 'altimg': '', 'w': '250', 'h': '54'}]的函式圖象fv
clear;
t=0:0.02:10;
f1=10/sqrt(7.75).*exp(-1.5*t);
f2=sin(sqrt(7.75).*t);
y=f1.*f2;
plot(t,y,'-k',t,y,'ok');
xlabel('t');ylabel('y(t) ');title('函式影象')
axis([-2 10 -0.5 2])
拉氏變換
clear;
clc;
syms s t fs1 fs2 fs3 ft1 ft2 ft3;
l=1,c=0.1,r=[1.5 3 5];
h1=1/(l*c*s^2+r(1)*c*s+1);
h2=1/(l*c*s^2+r(2)*c*s+1);
h3=1/(l*c*s^2+r(3)*c*s+1);
fs1=h1*(1/s);
fs2=h2*(1/s);
fs3=h3*(1/s);
ft1=ilaplace(fs1,s,t);
ft2=ilaplace(fs2,s,t);
ft3=ilaplace(fs3,s,t);
ezplot(t,ft1);
hold on;
ezplot(t,ft2);
hold on;
ezplot(t,ft3);
訊號編碼
對[1 1 0 1 1 1 0 1 0 0 1]進行編碼。
clear;
clc;
c=[1 1 0 1 1 1 0 1 0 0 1]
for i=1:length(c)
if i==1
d1(i)=0;d2(i)=0;
elseif i==2
d1(i)=c(i-1);d2(i)=c(i-1);
elseif i==3
d1(i)=mod(c(i-1)+c(i-2),2);
d2(i)=c(i-1);
else
d1(i)=mod(c(i-1)+c(i-2),2);
d2(i)=mod(c(i-1)+c(i-3),2);
endendd1
d2迭代法
使用一般迭代法求解方程的[=e^', 'altimg': '', 'w': '67', 'h': '21'}]解。
第一根:
clear;
clc;
x0=10;err=1;
while (err>10^-6)
f1=x0;
f2=log(3*x0^2);
err=abs(f1-f2);
x0=f2;
endx0
第二根:
clear;
clc;
x=1;err=1;
while (err>10^-6)
f1=x;
f2=sqrt(exp(x)/3);
err=abs(f1-f2);
x=f2;
endx
第三根:
clear;
clc;
x=10;err=1;
while (err>10^-6)
f1=x;
f2=-sqrt(exp(x)/3);
err=abs(f1-f2);
x=f2;
endx
牛頓迭代法
使用牛頓迭代法求解方程[4x^+3x6=0', 'altimg': '', 'w': '170', 'h': '21'}]的解。
syms x fx fx1;
fx=2*x^3-4*x^2+3*x-6;
fx1=diff(fx)
err=1;k=1;x0=1.5;
while (err>10^-6 && k<=1000)
x1=x0-(subs(fx,x,x0))/(subs(fx1,x,x0))
x0=x1;
err=abs(x0-x1)
k=k+1;
endk
x0牛頓迭代法求解
使用牛頓迭代法求解方程[12x=\\sin x1', 'altimg': '', 'w': '154', 'h': '21'}]的解。
第一解:
syms x fx fx1;
fx=x.^3-sin(x)-12*x+1;
fx1=diff(fx);
err=1;k=1;x0=-4;
while (err>10^-6 && k<=1000)
x1=x0-(subs(fx,x,x0))/(subs(fx1,x,x0));
err=abs(x0-x1);
x0=x1;
k=k+1;
endk
x0第二解:
syms x fx fx1;
fx=x.^3-sin(x)-12*x+1;
fx1=diff(fx);
err=1;k=1;x0=0;
while (err>10^-6) && (k<=1000)
x1=x0-(subs(fx,x,x0))/(subs(fx1,x,x0));
err=abs(x0-x1);
x0=x1;
k=k+1;
endk
x0第三解:
syms x fx fx1;
fx=x.^3-sin(x)-12*x+1;
fx1=diff(fx);
err=1;k=1;x0=4;
while (err>10^-6) && (k<=1000)
x1=x0-(subs(fx,x,x0))/(subs(fx1,x,x0));
err=abs(x0-x1);
x0=x1;
k=k+1;
endk
x0普通迭代法求解
使用普通迭代法求解方程[12x=\\sin x1', 'altimg': '', 'w': '154', 'h': '21'}]的解。
clear;
clc;
err1=1;err2=1;err3=1;x01=3;x02=0;x03=-3;k=3;
while (err1>10^-6) && (err2>10^-6) && (err3>10^-6) && (k<1000)
f11=x01;
f12=(sin(x01)+12*x01-1)^(1/3);
err1=abs(f11-f12);
x01=f12; %第一解
f21=x02;
f22=(x02^3-sin(x02)+1)/12;
err2=abs(f21-f22);
x02=f22; %第二解
f31=x03;
f32=-sqrt((sin(x03)-1)/x03+12);
err3=abs(f31-f32);
x03=f32; %第三解
k=k+1;
endx=[x01 x02 x03]
呼叫matlab函式求解方程。
clear;
clc;
syms x a b;
p=x.^2-a*x-4*b;
r=solve (p,x);
rclear; %roots求解法
clc;
syms x;
p=[1 0 0 0 -1 2 -3];
r=roots (p);
rclear; %solve求解法
clc;
syms x;
p=x.^6-x.^2+2*x-3;
r=solve (p,x);
rclear;clc;
syms x f;
f=@(x)x.*sin(x)-1;
x0=1;tol=1e-6;
z1=fzero (f,x0);
x0=2.7;tol=1e-6;
z2=fzero (f,x0);
z=[z1 z2]
clear;clc; %求非線性方程組的解
sym x;
fun='[x(1)-x(2)^2,x(2)-cos(x(1))]'
x0=[1,2]
f=fsolve(fun,x0)
clear;clc; %求解非線性方程的解
syms x;
fun='[sin(x)-log(x+0.1)]';
x0=1;
f=fsolve(fun,x0)
矩陣基本操作
clear;clc; %加、減、乘省略,只需補零即可運算
a=[1 2 3 4;2 3 1 2;1 1 1 -1;1 0 -2 -6];
b=[5;6;7;8];
c1=a\b左除
c2=det(a) %求行列式
c3=inv(a) %求逆
c4=a^3求冪
c5=rank(a) %求秩
c6=a求轉置
c7=rref(a) %行變換簡
求解方程[4x_+2x_x_=2\\\\ 3x_x_+2x_=10\\\\ 11x_+3x_=8\\end\\right.', 'altimg': '', 'w':
'162', 'h': '123'}]的解。
clear;clc;
a=[4 2 -1;3 -1 2;11 3 0];
b=[2;10;8];
rref(a)
rref([a,b])
r1=rank(a)
r2=rank([a,b])
求解方程[2x+3y+z=4\\\\ x2y+4z=5\\\\ 3x+8y2z=13\\\\ 4xy+9z=6\\end\\right.', 'altimg': '', 'w':
'151', 'h': '150'}]的解。
clear;clc;
a=[2 3 1 0;1 -2 4 0;3 8 -2 0;4 -1 9 0];
b=[4;-5;13;-6];
rref(a)
rref([a,b])
matlab的程式設計例項解析
matlab程式設計方法及若干程式例項 樊双喜 河南大學數學與資訊科學學院開封475004 摘要本文通過對 matlab 程式設計中的若干典型問題做簡要的分析和總結,並 在此基礎上著重討論了有關演算法設計 程式的除錯與測試 演算法與程式的優化以 及迴圈控制等方面的問題.還通過對一些程式例項做具體解析...
學習MATLAB感想
學習了matlab這門課程,我了解該軟體的基本功能,也知道了該軟體在我們生活中的重要地位。隨著社會的不斷發展,科技的不斷進步,計算機的普及,它也被應用在越來越多的方面。matlab的基本資料單位是矩陣,它的指令表示式與數學 工程中常用的形式十分相似,故用matlab來解算問題要比用c,fortran...
matlab學習總結
matlab總結 一 matlab常用函式 1 特殊變數與常數 2 操作符與特殊字元 3 基本數學函式 4 基本矩陣和矩陣操作 5 數值分析和傅利葉變換 6 多項式與插值 7 繪圖函式 二 matlab工作間常用命令 1 常用的視窗命令 2 有關檔案及其操作的語句 3 啟動與退出的命令 4 管理變數...